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% 133.9 / NA 33intergenic (+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 (19/14);  total (19/14)

GGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACC  >  NC_000913/257749‑257839
                                             |                                             
ggATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                            >  1:2028141/1‑49 (MQ=255)
  aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                          <  1:2033171/49‑1 (MQ=255)
   tCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaa                                         <  1:532592/49‑1 (MQ=255)
   tCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                          >  1:1226647/1‑48 (MQ=255)
     gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       <  1:1896002/49‑1 (MQ=255)
     gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:1660344/1‑49 (MQ=255)
       aaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:1916809/1‑47 (MQ=255)
       aaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:386086/1‑47 (MQ=255)
       aaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACa                                     <  1:1373860/49‑1 (MQ=255)
         aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                   >  1:2080074/1‑49 (MQ=255)
         aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                   >  1:892180/1‑49 (MQ=255)
          aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGtt                                  >  1:68404/1‑49 (MQ=255)
            gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                 >  1:1538725/1‑48 (MQ=255)
            gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                 >  1:1357635/1‑48 (MQ=255)
              ttGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:1808599/1‑49 (MQ=255)
                gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATc                             <  1:26208/48‑1 (MQ=255)
                    tAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                         >  1:1391564/1‑48 (MQ=255)
                     aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                         >  1:1542948/1‑47 (MQ=255)
                      aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      <  1:1408830/49‑1 (MQ=255)
                      aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:55099/1‑49 (MQ=255)
                         aTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgag                   <  1:1937451/49‑1 (MQ=255)
                         aTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgag                   <  1:1634029/49‑1 (MQ=255)
                          tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgag                   <  1:226131/48‑1 (MQ=255)
                           ttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                  >  1:822878/1‑48 (MQ=255)
                            tGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                  >  1:659576/1‑47 (MQ=255)
                            tGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                  >  1:495829/1‑47 (MQ=255)
                             gCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg               >  1:1365689/1‑49 (MQ=255)
                              cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGc              <  1:1246853/49‑1 (MQ=255)
                                 aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCa             <  1:1404317/47‑1 (MQ=255)
                                   ttttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGa         <  1:395094/49‑1 (MQ=255)
                                     ttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc        <  1:105178/48‑1 (MQ=255)
                                         cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc   >  1:450398/1‑49 (MQ=255)
                                          aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc  <  1:1786149/49‑1 (MQ=255)
                                             |                                             
GGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACC  >  NC_000913/257749‑257839

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