Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 950,036 (T)5→4 intergenic (‑114/+292) focA ← / ← ycaO formate transporter FocA/30S ribosomal protein S12 methylthiotransferase accessory factor YcaO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273950,0320T.100.0% 59.5 / NA 15intergenic (‑110/+296)focA/ycaOformate transporter FocA/30S ribosomal protein S12 methylthiotransferase accessory factor YcaO
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (8/7);  total (8/7)

TCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTTGAATTTTTTCCGCA  >  NZ_CP009273/949944‑950117
                                                                                        |                                                                                     
tCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGtt                                                                                      >  1:134603/1‑90 (MQ=255)
    cgacgaGATACTAACAAAGCATTATAGATGAGAAACTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGttttat                                                                                  <  2:715172/90‑3 (MQ=255)
             aCTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACg                                                                        <  2:342432/90‑1 (MQ=255)
              cTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACgg                                                                       <  2:392613/90‑1 (MQ=255)
                  cAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTa                                                                   <  1:443084/90‑1 (MQ=255)
                  cAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTa                                                                   <  2:605033/90‑1 (MQ=255)
                       cATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTAAatat                                                              <  1:632045/90‑1 (MQ=255)
                                tGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATc                                                     >  1:410177/1‑90 (MQ=255)
                                                    tatCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTaa                                 >  1:427589/1‑90 (MQ=255)
                                                      tctcGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTaaaa                               >  2:669455/1‑90 (MQ=255)
                                                         cGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTa                                                                   >  1:638315/1‑51 (MQ=255)
                                                         cGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTa                                                                   <  2:638315/51‑1 (MQ=255)
                                                                         ccGTAACTCGCATGGTTTTATGCAAATACGGAGTAAatat                                                              >  2:716774/1‑40 (MQ=38)
                                                                         ccGTAACTCGCATGGTTTTATGCAAATACGGAGTAAatat                                                              <  1:716774/40‑1 (MQ=38)
                                                                         ccGTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTTGAAt            >  1:478042/1‑90 (MQ=255)
                                                                           gTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTTGAAttt          >  2:163870/1‑90 (MQ=255)
                                                                                   cATGGTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTTGAATTTTTTCCGCa  >  2:67946/1‑90 (MQ=255)
                                                                                        |                                                                                     
TCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTTGAATTTTTTCCGCA  >  NZ_CP009273/949944‑950117

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