Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,451,613 C→T intergenic (+51/‑617) 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,6130CT100.0% 70.1 / NA 26intergenic (+51/‑617)unknown/unknownFIG00698722: hypothetical protein/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/14);  total (12/14)

TCGGTGCTGATGTGCCGGTGATTTCGGATTTTACAGAACGAAACGGAAAGATTATTCAATCTGTTCAATATCGTCAAAGCGGGGTGATTTTTAATATTGAGCCGGTTGTTAAAATTAATGCGATAGATTTAAAAATTAATCAACAGCTATCCAATTTTATTAAGACGGATAAAGGCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACACTTCTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATCACGTTTAACTAGTAATAGTGACCTGTGAATAACTCTGTGTTTTCATTAAATCGTAGTGTTCATTTACTTTGTCACTAATTTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGGTTTGCTCTCTTTAAAAATGA  >  1/1451381‑1451837
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                
attttactccaagtgaaacttgtctaactatttttgagcaagctcaaaaatgtagggcgaaggcttgcagcccggtaaaccatgctcgatgttattcgggcatttttattttttgcgcgtacctaaatttaatcgagaaatgaaaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTc                                                                                                                                                                                                                        <  1:426022‑M2/13‑1 (MQ=255)
attttactccaagtgaaacttgtctaactatttttgagcaagctcaaaaatgtagggcgaaggcttgcagcccggtaaaccatgctcgatgttattcgggcatttttattttttgcgcgtacctaaatttaatcgagaaatgaaaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATaaa                                                                                                                                                                                                                <  1:268366‑M2/21‑1 (MQ=255)
 ttttactccaagtgaaacttgtctaactatttttgagcaagctcaaaaatgtagggcgaaggcttgcagcccggtaaaccatgctcgatgttattcgggcatttttattttttgcgcgtacctaaatttaatcgagaaatgaaaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAAt                                                                                                                                                                                                               >  1:784907‑M2/230‑251 (MQ=255)
          aagtgaaacttgtctaactatttttgagcaagctcaaaaatgtagggcgaaggcttgcagcccggtaaaccatgctcgatgttattcgggcatttttattttttgcgcgtacctaaatttaatcgagaaatgaaaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAAtattta                                                                                                                                                                                                      >  1:844796‑M2/221‑251 (MQ=255)
                         aactatttttgggcaagctcaaaaatgtagggcgaaggcttgcagcccggtaaaccatgctcgatgttattcgggcatttttattttttgcgcgtacctaaatttaatcgagaaatgaaaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAgg                                                                                                                                                                                        >  1:346280‑M2/206‑250 (MQ=255)
                         aactatttttgagcaagctcaaaaatgtagggcgaaggcttgcagcccggtaaaccatgctcgatgttattcgggcatttttattttttgcgcgtacctaaatttaatcgagaaatgaaaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGc                                                                                                                                                                                       <  1:236940‑M2/46‑1 (MQ=255)
                                                         cgaaggcttgcagcccggtaaaccatgctcgatgttattcgggcatttttattttttgcgcgtacctaaatttaatcgagaaatgaaaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAacc                                                                                                                                                        <  1:367120‑M2/77‑2 (MQ=255)
                                                                     gcccggtaaaccatgctcgatgttattcgggcatttttattttttgcgcgtacctaaatttaatcgagaaatgaaaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACaa                                                                                                                                                                                          <  1:828074‑M2/43‑1 (MQ=255)
                                                                      cccggtaaaccatgctcgatgttattcgggcatttttattttttgcgcgtacctaaatttaatcgagaaatgaaaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAatat                                                                                                                                                                                                         <  1:530159‑M2/28‑1 (MQ=255)
                                                                                                    catttttattttttgcgcgtacctaaatttaatcgagaaatgaaaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATa                                                                                                                                                                                                                            >  1:623113‑M2/131‑139 (MQ=255)
                                                                                                      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)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                
TCGGTGCTGATGTGCCGGTGATTTCGGATTTTACAGAACGAAACGGAAAGATTATTCAATCTGTTCAATATCGTCAAAGCGGGGTGATTTTTAATATTGAGCCGGTTGTTAAAATTAATGCGATAGATTTAAAAATTAATCAACAGCTATCCAATTTTATTAAGACGGATAAAGGCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACACTTCTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATCACGTTTAACTAGTAATAGTGACCTGTGAATAACTCTGTGTTTTCATTAAATCGTAGTGTTCATTTACTTTGTCACTAATTTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGGTTTGCTCTCTTTAAAAATGA  >  1/1451381‑1451837

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