Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,220,940 A→G E220G (GAA→GGA)  rluF → 23S rRNA pseudouridine(2604) synthase RluF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,220,9400AG100.0% 84.0 / NA 24E220G (GAA→GGA) rluF23S rRNA pseudouridine(2604) synthase RluF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/13);  total (11/13)

GCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAA  >  NZ_CP009273/4220871‑4221018
                                                                     |                                                                              
gCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACcgacg                                                            >  1:30363/1‑90 (MQ=255)
 cTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACcgacga                                                           <  1:156528/90‑1 (MQ=255)
  tATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGAt                                                          <  1:230062/90‑1 (MQ=255)
     gAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAg                                                       <  1:399767/90‑1 (MQ=255)
          gAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAAt                                                  >  1:396818/1‑90 (MQ=255)
           aaaaaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATc                                                 <  1:270836/90‑1 (MQ=255)
           aaaaaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATc                                                 <  2:30363/90‑1 (MQ=255)
              aaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGAc                                              <  1:353536/90‑1 (MQ=255)
                  tGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGACctct                                          >  1:410286/1‑90 (MQ=255)
                  tGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGACctct                                          >  2:222346/1‑90 (MQ=255)
                  tGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGACctct                                          >  1:449106/1‑90 (MQ=255)
                            gcgcATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCAt                                <  2:204310/90‑1 (MQ=255)
                                     gAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGAcc                                             >  2:94496/1‑68 (MQ=255)
                                     gAACGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGAcc                                             <  1:94496/68‑1 (MQ=255)
                                       aCGTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCct                     >  1:510719/1‑90 (MQ=255)
                                         gTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGAc                                              <  1:387883/63‑1 (MQ=255)
                                         gTTAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGAc                                              >  2:387883/1‑63 (MQ=255)
                                           tAGCTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTc                 >  1:448706/1‑90 (MQ=255)
                                              cTTAAGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGa              <  2:410286/90‑1 (MQ=255)
                                                 aaGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGt           >  2:178608/1‑90 (MQ=255)
                                                 aaGCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGt           >  1:51054/1‑90 (MQ=255)
                                                   gCGGCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTaa         <  2:396818/90‑1 (MQ=255)
                                                     ggCATTCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATCCCTCTCCCGAGGTaaaa       <  1:296092/90‑1 (MQ=255)
                                                          tCCGCTGGGGGGATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTaa  <  2:129709/90‑1 (MQ=255)
                                                                     |                                                                              
GCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAA  >  NZ_CP009273/4220871‑4221018

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: