Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 157,096 +G intergenic (‑26/+41) h16_B0136 ← / ← h16_B0137 cysteine desulfurase/conserved hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260480157,0941.G100.0% 132.9 / NA 35intergenic (‑24/+43)h16_B0136/h16_B0137cysteine desulfurase/conserved hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (19/16);  total (19/16)

TCGGGCGGCAACGAGGCCGCGGCGGTGTCGAGATAGATGGTGGAGGCCGTCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAGCCGAAGAAGCTGAACAGGGCGTCGCCGGCAATCACGCCGCCGGCGAAGGCGCTCATC  >  AM260480/156976‑157218
                                                                                                                       |                                                                                                                            
tCGGGCGGCAACGAGGCCGCGGCGGTGTCGAGATAGATGGTGGAGGCCGTCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACcg                                                                                                                >  1:379008/1‑134 (MQ=255)
  gggCGGCAACGAGGCCGCGGCGGTGTCGAGATAGATGGTGGAGGCCGTCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACcgcg                                                                                                              >  2:1009134/1‑134 (MQ=255)
  gggCGGCAACGAGGCCGCGGCGGTGTCGAGATAGATGGTGGAGGCCGTCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACcgcg                                                                                                              >  1:976179/1‑134 (MQ=255)
        cAACGAGGCCGCGGCGGTGTCGAGATAGATGGTGGAGGCCGTCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGa                                                                                                        <  2:976179/134‑1 (MQ=255)
                  gcggcgGTGTCGAGATAGATGGTGGAGGCCGTCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACt                                                                                               >  2:315316/1‑133 (MQ=255)
                                    atGGTGGAGGCCGTCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCgcgc                                                                                                             <  2:583846/101‑1 (MQ=255)
                                    atGGTGGAGGCCGTCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCgcgc                                                                                                             >  1:583846/1‑101 (MQ=255)
                                            ggCCGTCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTg                                                                    >  1:1033623/1‑134 (MQ=255)
                                              ccGTCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGcccc                                                                                                          <  1:299144/94‑1 (MQ=255)
                                              ccGTCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGcccc                                                                                                          >  2:299144/1‑94 (MQ=255)
                                                 tCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGAc                                                               >  2:795309/1‑134 (MQ=255)
                                                      ggggTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGt                                                                  >  1:591203/1‑126 (MQ=255)
                                                      ggggTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGt                                                                  <  2:591203/126‑1 (MQ=255)
                                                            tCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAGCCgaag                                                     <  2:1033623/133‑1 (MQ=255)
                                                             cGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAGCCgaag                                                     >  1:1044285/1‑132 (MQ=255)
                                                                 gcgGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAGCCGAAGAAGCt                                                >  2:223886/1‑133 (MQ=255)
                                                                                   ttGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGt                                                                                                     <  1:249266/62‑1 (MQ=255)
                                                                                   ttGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGt                                                                                                     >  2:249266/1‑62 (MQ=255)
                                                                                       atatCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAGCCGAAGAAGCTGAACAGGGCGTCGCCGGCAATCa                         <  1:223886/134‑1 (MQ=255)
                                                                                          tCCATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAGCCGAAGAAGCTGAACAGGGCGTCGCCGGCAATCAcgc                      <  1:795309/134‑1 (MQ=255)
                                                                                           ccATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCactta                                                                                 >  2:707281/1‑74 (MQ=255)
                                                                                           ccATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCactta                                                                                 <  1:707281/74‑1 (MQ=255)
                                                                                            cATGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAGCCGAAGAAGCTGAACAGGGCGTCGCCGGCAATCAcgccg                    >  2:227812/1‑134 (MQ=255)
                                                                                              tGGAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAGCCGAAGAAGCTGAACAGGGCGTCGCCGGCAATCACgccgcc                  <  1:227812/134‑1 (MQ=255)
                                                                                               ggAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAg                                                           <  1:719724/92‑1 (MQ=255)
                                                                                               ggAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAg                                                           >  2:719724/1‑92 (MQ=255)
                                                                                                gAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACcgc                                                                                                               >  1:288012/1‑39 (MQ=37)
                                                                                                gAACGTCTCGGGTAATGACAGGCGGGACGCGGCGACcgc                                                                                                               <  2:288012/39‑1 (MQ=37)
                                                                                                    gTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAGCCGAAGAAGCTGAACAGGGCGTCGCCGGCAATCACGCCGCCGGCg              >  1:66705/1‑132 (MQ=255)
                                                                                                    gTCTCGGGTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAGCCGAAGAAGCTGAACAGGGCGTCGCCGGCAATCACGCCGCCGGCGaa            <  1:315316/134‑1 (MQ=255)
                                                                                                         gggTAATGACAGGCGGGACGCGGCGACCGCGcccc                                                                                                          <  1:69344/35‑1 (MQ=25)
                                                                                                         gggTAATGACAGGCGGGACGCGGCGACCGCGcccc                                                                                                          >  2:69344/1‑35 (MQ=25)
                                                                                                          ggTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGt                                                                                           <  1:93994/49‑1 (MQ=38)
                                                                                                          ggTAATGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGt                                                                                           >  2:93994/1‑49 (MQ=38)
                                                                                                              aTGACAGGCGGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAGCCGAAGAAGCTGAACAGGGCGTCGCCGGCAATCACGCCGCCGGCGAAGGCGCTCAtc  <  1:1009134/134‑1 (MQ=255)
                                                                                                                       |                                                                                                                            
TCGGGCGGCAACGAGGCCGCGGCGGTGTCGAGATAGATGGTGGAGGCCGTCAGCGGGGTGTCGGCGCGGATCGCCTGGATGGCTTGGATATCCATGGAACGTCTCGGGTAATGACAGGCGGACGCGGCGACCGCGCCCCGAGGTAAGACTCAGTTGCGCACTTACTTGCGCAGCTTGGTGACGGAGCCGAAGAAGCTGAACAGGGCGTCGCCGGCAATCACGCCGCCGGCGAAGGCGCTCATC  >  AM260480/156976‑157218

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