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% 112.2 / NA 31intergenic (+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 (18/13);  total (18/13)

TGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACC  >  NC_000913/257748‑257839
                                              |                                             
tGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                            >  1:1669751/1‑50 (MQ=255)
 ggATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                            >  1:917631/1‑49 (MQ=255)
  gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                           <  1:2377664/49‑1 (MQ=255)
  gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                           <  1:2211106/49‑1 (MQ=255)
   aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                          >  1:2636291/1‑49 (MQ=255)
      gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:1433082/1‑49 (MQ=255)
      gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       <  1:2010157/49‑1 (MQ=255)
       aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACa                                     <  1:2073215/50‑1 (MQ=255)
        aaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAg                                    >  1:2602515/1‑50 (MQ=255)
           aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                 >  1:2626670/1‑50 (MQ=255)
             gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                               >  1:800878/1‑50 (MQ=255)
             gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                               <  1:1177421/50‑1 (MQ=255)
                tgtgTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                            >  1:1066681/1‑50 (MQ=255)
                  tgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATcac                           >  1:1736469/1‑49 (MQ=255)
                    ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                         >  1:1935680/1‑49 (MQ=255)
                     tAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                         >  1:700153/1‑48 (MQ=255)
                      aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:1041913/1‑50 (MQ=255)
                       aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAgg                     <  1:1204127/50‑1 (MQ=255)
                       aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      <  1:926620/49‑1 (MQ=255)
                       aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:112549/1‑49 (MQ=255)
                        aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:2459252/1‑48 (MQ=255)
                          aTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:2621771/1‑46 (MQ=255)
                           tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGAt                 <  1:1437386/50‑1 (MQ=255)
                              gCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg               >  1:1388392/1‑49 (MQ=255)
                                  aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATg          <  1:789537/50‑1 (MQ=255)
                                    ttttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc        <  1:2220686/50‑1 (MQ=255)
                                       tGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt      <  1:322234/49‑1 (MQ=255)
                                        gCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGtt     <  1:1106381/49‑1 (MQ=255)
                                        gCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTa    >  1:1226018/1‑50 (MQ=255)
                                         ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc   >  1:959284/1‑50 (MQ=255)
                                          cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc  <  1:201603/50‑1 (MQ=255)
                                              |                                             
TGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACC  >  NC_000913/257748‑257839

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