Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,306,087 A→G intergenic (+22/+18) tonB → / ← yciA TonB system transport protein TonB/acyl‑CoA thioester hydrolase YciA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,306,0870AG100.0% 58.3 / NA 18intergenic (+22/+18)tonB/yciATonB system transport protein TonB/acyl‑CoA thioester hydrolase YciA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/10);  total (8/10)

AGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCA  >  NZ_CP009273/1305959‑1306221
                                                                                                                                |                                                                                                                                      
aGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCtt                                                                                                                              <  2:271718/139‑1 (MQ=255)
         atgcgatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTAtt                                                                                                                     <  2:124844/139‑1 (MQ=255)
                                      ggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGt                                                                                                          >  1:47990/1‑121 (MQ=255)
                                      ggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGt                                                                                                          <  2:47990/121‑1 (MQ=255)
                                                gCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCaa                                                                              >  2:314219/1‑139 (MQ=255)
                                                        aTTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACa                                                                      <  1:116146/139‑1 (MQ=255)
                                                            tggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACt                                                                  >  2:15417/1‑139 (MQ=255)
                                                                         gTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTc                                                     >  1:212124/1‑139 (MQ=255)
                                                                          tttAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTc                                                     <  1:425391/138‑1 (MQ=255)
                                                                               aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGAc                                                                                                                         <  1:81236/65‑1 (MQ=255)
                                                                               aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGAc                                                                                                                         >  2:81236/1‑65 (MQ=255)
                                                                                        caccacCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCg                                      <  2:280900/139‑1 (MQ=255)
                                                                                           caccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTg                                   >  1:345899/1‑139 (MQ=255)
                                                                                                 aaTTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAAt                             <  2:345899/139‑1 (MQ=255)
                                                                                                        tAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGtt                        >  2:49483/1‑137 (MQ=255)
                                                                                                          aGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGa                    <  2:82478/139‑1 (MQ=255)
                                                                                                                    cAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACtttt          >  2:16052/1‑139 (MQ=255)
                                                                                                                          ccTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCAcc    <  1:15417/139‑1 (MQ=255)
                                                                                                                             ccggccggAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCca  >  2:423281/1‑138 (MQ=255)
                                                                                                                             ccggccggAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCca  >  2:400047/1‑138 (MQ=255)
                                                                                                                                |                                                                                                                                      
AGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCA  >  NZ_CP009273/1305959‑1306221

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