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% 94.6 / NA 26intergenic (+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 (13/13);  total (13/13)

GGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGT  >  NZ_CP009273/1306006‑1306173
                                                                                 |                                                                                      
ggCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCt                                                                               >  1:336505/1‑91 (MQ=255)
ggCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCt                                                                               >  1:442458/1‑91 (MQ=255)
  cAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCttt                                                                             <  1:659971/91‑1 (MQ=255)
           tgtgGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTa                                                                    >  2:196528/1‑91 (MQ=255)
              ggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTg                                                                           >  1:216860/1‑81 (MQ=255)
              ggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTg                                                                           <  2:216860/81‑1 (MQ=255)
                 gAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAAc                                                              <  1:223620/91‑1 (MQ=255)
                      tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTa                                                         >  2:335198/1‑91 (MQ=255)
                       ccTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTaa                                                        <  1:538165/91‑1 (MQ=255)
                         tGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAg                                                                                  >  1:102732/1‑63 (MQ=255)
                         tGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAg                                                                                  <  2:102732/63‑1 (MQ=255)
                                                        gTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGAc                       <  1:335198/91‑1 (MQ=255)
                                                        gTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGAc                       <  1:719911/91‑1 (MQ=255)
                                                          aaGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGgcgc                                                   <  2:524800/61‑1 (MQ=255)
                                                          aaGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGgcgc                                                   >  1:524800/1‑61 (MQ=255)
                                                                  aGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaat             >  2:585481/1‑91 (MQ=255)
                                                                       aaaGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGc        >  2:440204/1‑91 (MQ=255)
                                                                          gCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTc     >  2:266185/1‑91 (MQ=255)
                                                                             tCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtttt                                           >  1:559512/1‑50 (MQ=255)
                                                                             tCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtttt                                           >  2:258968/1‑50 (MQ=255)
                                                                             tCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtttt                                           >  2:287298/1‑50 (MQ=255)
                                                                             tCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtttt                                           <  1:287298/50‑1 (MQ=255)
                                                                             tCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtttt                                           <  2:559512/50‑1 (MQ=255)
                                                                             tCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtttt                                           <  1:258968/50‑1 (MQ=255)
                                                                             tCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGt  <  1:585481/91‑1 (MQ=255)
                                                                             tCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGt  <  1:517999/91‑1 (MQ=255)
                                                                             tCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGt  <  1:277983/91‑1 (MQ=255)
                                                                                 |                                                                                      
GGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGT  >  NZ_CP009273/1306006‑1306173

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