Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,451,693 T→G intergenic (+131/‑537) unknown → / → unknown FIG00698722: hypothetical protein/Bacterial non‑heme ferritin (EC 1.16.3.2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,451,6930TG100.0% 72.5 / NA 25intergenic (+131/‑537)unknown/unknownFIG00698722: hypothetical protein/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/12);  total (13/12)

CGGTTGTTAAAATTAATGCGATAGATTTAAAAATTAATCAACAGCTATCCAATTTTATTAAGACGGATAAAGGCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACACTTCTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATCACGTTTAACTAGTAATAGTGACCTGTGAATAACTCTGTGTTTTCATTAAATCGTAGTGTTCATTTACTTTGTCACTAATTTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGGTTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTC  >  1/1451483‑1451932
                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
tttttattttttgcgcgtacctaaatttaatcgagaaatgaaaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAAc                                                                                                                                                                                                                                >  1:192171‑M2/129‑228 (MQ=255)
                 tacctaaatttaatcgagaaatgaaaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAAt                                                                                                                                                                                                                                   >  1:749381‑M2/112‑208 (MQ=255)
                                         aaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgactagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatta                                                                                                                                                                <  1:358946‑M2/164‑1 (MQ=255)
                                            tttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGt                                                                                                                                                             <  1:641107‑M2/167‑1 (MQ=255)
                                                            atacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATgg                                                                                                                                             <  1:334688‑M2/183‑1 (MQ=255)
                                                               cccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAAc                                                                                                                                                                                                                                >  1:520655‑M2/66‑165 (MQ=255)
                                                               cccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatt                                                                                                                                                                 <  1:332040‑M2/163‑1 (MQ=255)
                                                               cccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTaa                                                                                                                                          >  1:814500‑M2/66‑251 (MQ=255)
                                                                                     gctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatt                                                                                                                                                                 >  1:591573‑M2/44‑206 (MQ=255)
                                                                                       tgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAActc                                                                                                                                                                                                                              <  1:562196‑M2/102‑1 (MQ=255)
                                                                                        gattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACatatat                                                                                                                                                                  <  1:873395‑M2/162‑1 (MQ=255)
                                                                                        gattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatta                                                                                                                                                                <  1:882711‑M2/164‑1 (MQ=255)
                                                                                            agataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACatatat                                                                                                                                                                  >  1:413914‑M2/37‑198 (MQ=255)
                                                                                                 atgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTgtgt                                                                                                                                                                                                                         >  1:182044‑M2/32‑138 (MQ=255)
                                                                                                      atacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTaaaaa                                                                                                    <  1:732398‑M2/224‑1 (MQ=255)
                                                                                                         cgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGa                                                                                                 >  1:684168‑M2/24‑250 (MQ=255)
                                                                                                                                          gTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAAc                                                                                                                                                                                                                                <  1:324039/90‑1 (MQ=255)
                                                                                                                                                        atatTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGt                                                                                                                                                             >  1:796235/1‑143 (MQ=255)
                                                                                                                                                                  aTGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTaatccgtagcg                                                                                                                                                                                                       >  1:502915/1‑82 (MQ=255)
                                                                                                                                                                     ttACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatt                                                                                                                                                                 <  1:659191/126‑1 (MQ=255)
                                                                                                                                                                     ttACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTaaaaa                                                                                                    >  1:661888/1‑187 (MQ=255)
                                                                                                                                                                            gCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACat                                                                                                                                                                      >  1:676356/1‑114 (MQ=255)
                                                                                                                                                                                    tcttcAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTag                     <  1:352855/251‑1 (MQ=255)
                                                                                                                                                                                     cttcAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTaga                    >  1:313791/1‑251 (MQ=255)
                                                                                                                                                                                                       ttAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTtc  <  1:807161/251‑1 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
CGGTTGTTAAAATTAATGCGATAGATTTAAAAATTAATCAACAGCTATCCAATTTTATTAAGACGGATAAAGGCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACACTTCTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATCACGTTTAACTAGTAATAGTGACCTGTGAATAACTCTGTGTTTTCATTAAATCGTAGTGTTCATTTACTTTGTCACTAATTTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGGTTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTC  >  1/1451483‑1451932

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑