Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 257,794 T→C intergenic (+23/‑35) frsA → / → crl fermentation‑respiration switch protein/RNA polymerase holoenzyme assembly factor Crl

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913257,7940TC100.0% 122.7 / NA 32intergenic (+23/‑35)frsA/crlfermentation‑respiration switch protein/RNA polymerase holoenzyme assembly factor Crl
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/11);  total (21/11)

CTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACC  >  NC_000913/257747‑257839
                                               |                                             
cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                             <  1:924949/50‑1 (MQ=255)
 tGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                            >  1:822561/1‑50 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                            <  1:1252202/48‑1 (MQ=255)
      cGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:1749663/1‑50 (MQ=255)
      cGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:1464486/1‑50 (MQ=255)
          aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                   <  1:1065401/50‑1 (MQ=255)
             cGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGc                                <  1:1635944/50‑1 (MQ=255)
               cTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:982904/1‑50 (MQ=255)
                  gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                            >  1:971179/1‑49 (MQ=255)
                  gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATc                             <  1:1552562/48‑1 (MQ=255)
                   tgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATcac                           >  1:2184158/1‑49 (MQ=255)
                     ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                         >  1:2326438/1‑49 (MQ=255)
                       aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:1356877/1‑50 (MQ=255)
                        aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:593531/1‑49 (MQ=255)
                        aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      <  1:2097512/49‑1 (MQ=255)
                         aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:1115839/1‑48 (MQ=255)
                         aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:448948/1‑48 (MQ=255)
                         aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGga                    <  1:1134041/50‑1 (MQ=255)
                            tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                  >  1:1358398/1‑49 (MQ=255)
                             ttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                  >  1:1609240/1‑48 (MQ=255)
                                cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGc              >  1:1953361/1‑49 (MQ=255)
                                  aaaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt           >  1:2321464/1‑50 (MQ=255)
                                   aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt           <  1:620318/49‑1 (MQ=255)
                                   aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATg          <  1:2306862/50‑1 (MQ=255)
                                    aTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGa         >  1:930690/1‑50 (MQ=255)
                                       ttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt      >  1:382239/1‑50 (MQ=255)
                                       ttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt      >  1:1314063/1‑50 (MQ=255)
                                        tGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt      >  1:2281357/1‑49 (MQ=255)
                                          ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTa    <  1:914377/49‑1 (MQ=255)
                                          ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc   >  1:25622/1‑50 (MQ=255)
                                           cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc  <  1:2199115/50‑1 (MQ=255)
                                            aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc  >  1:967948/1‑49 (MQ=255)
                                               |                                             
CTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACC  >  NC_000913/257747‑257839

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