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% 67.2 / NA 17intergenic (‑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 . (10/7);  total (10/7)

TCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTTGAATTTTTTCCGC  >  NZ_CP009273/949944‑950116
                                                                                        |                                                                                    
tCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGtt                                                                                     >  2:251155/1‑90 (MQ=255)
      acgaGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGttttatgc                                                                               <  1:63621/90‑5 (MQ=255)
            tACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATAc                                                                        <  2:750/90‑1 (MQ=255)
            tACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATAc                                                                        <  1:251155/90‑1 (MQ=255)
            tACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATAc                                                                        <  1:288842/90‑1 (MQ=255)
             aCTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACg                                                                       >  1:217429/1‑90 (MQ=255)
             aCTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACg                                                                       <  2:175606/90‑1 (MQ=255)
             aCTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACg                                                                       <  1:706031/90‑1 (MQ=255)
                        aTTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAGATACgg                                                                      >  1:617300/1‑80 (MQ=255)
                        aTTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACgg                                                                      <  2:617300/80‑1 (MQ=255)
                                tGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATc                                                    >  1:21594/1‑90 (MQ=255)
                                tGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATc                                                    >  2:80727/1‑90 (MQ=255)
                                   gaAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTACCCaa                                                 >  2:744979/1‑90 (MQ=255)
                                                aTCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAAtt                                    >  2:337158/1‑90 (MQ=255)
                                                         cGCCTGGCTTATAGGCCCGTAACTCGCATGGtttt                                                                                   >  1:579521/1‑35 (MQ=255)
                                                         cGCCTGGCTTATAGGCCCGTAACTCGCATGGtttt                                                                                   <  2:579521/35‑1 (MQ=255)
                                                         cGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAAtt                           >  1:558368/1‑90 (MQ=255)
                                                         cGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAAtt                           >  2:31390/1‑90 (MQ=255)
                                                          gCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTa                          >  2:653670/1‑90 (MQ=255)
                                                                        cccGTAACTCGCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTTGaa            >  2:96846/1‑90 (MQ=255)
                                                                                  gCATGGTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTTGAATTTTTTCCGc  <  1:31390/90‑1 (MQ=255)
                                                                                        |                                                                                    
TCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTAAAAATTAACAAATGAGTTGAATTTTTTCCGC  >  NZ_CP009273/949944‑950116

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: