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% 156.9 / NA 39intergenic (+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 (25/14);  total (25/14)

TGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAC  >  NC_000913/257748‑257838
                                              |                                            
tGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                            >  1:1544044/1‑49 (MQ=255)
  gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                          <  1:265022/49‑1 (MQ=255)
   aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                          >  1:56203/1‑48 (MQ=255)
   aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                         <  1:1190222/49‑1 (MQ=255)
      gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaa                                       <  1:789274/48‑1 (MQ=255)
      gAAAAACGCTTGTGTTAAAAATTTGCTAAAGTTTGCCAATCTGGTaaaa                                      <  1:1491152/49‑1 (MQ=38)
       aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                      >  1:590989/1‑48 (MQ=255)
       aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                      >  1:132952/1‑48 (MQ=255)
       aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                      >  1:506619/1‑48 (MQ=255)
         aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAg                                   >  1:672511/1‑49 (MQ=255)
          aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                  >  1:1671949/1‑49 (MQ=255)
          aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                  <  1:494555/49‑1 (MQ=255)
          aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                  >  1:18419/1‑49 (MQ=255)
           aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGtt                                 >  1:219738/1‑49 (MQ=255)
            cGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                <  1:1588849/49‑1 (MQ=255)
            cGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                >  1:545904/1‑49 (MQ=255)
             gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                >  1:1270889/1‑48 (MQ=255)
              cTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                              >  1:1417198/1‑49 (MQ=255)
              cTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                              <  1:1737026/49‑1 (MQ=255)
                 gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                           >  1:861100/1‑49 (MQ=255)
                 gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                           >  1:128105/1‑49 (MQ=255)
                 gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGCAAAACAGTTGCATca                           >  1:170546/1‑49 (MQ=39)
                  tgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                           >  1:1291961/1‑48 (MQ=255)
                    ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                        <  1:1936202/49‑1 (MQ=255)
                      aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaac                       <  1:424114/48‑1 (MQ=255)
                        aaaTTTGCTAAATTTTGTCAATCTGGTAAAACAGTTGCATCACAACAgg                    <  1:449741/49‑1 (MQ=38)
                           tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                     >  1:248983/1‑45 (MQ=255)
                              gCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg              >  1:1856465/1‑49 (MQ=255)
                               cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg              >  1:1469126/1‑48 (MQ=255)
                               cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg              >  1:641671/1‑48 (MQ=255)
                                 aaaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCaa           >  1:127909/1‑49 (MQ=255)
                                 aaaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCaa           >  1:271282/1‑49 (MQ=255)
                                  aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCaa           <  1:1094704/48‑1 (MQ=255)
                                       tGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt     >  1:350394/1‑49 (MQ=255)
                                         ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGtt    <  1:1631962/48‑1 (MQ=255)
                                         ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTa   >  1:1774666/1‑49 (MQ=255)
                                         ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTa   <  1:2057920/49‑1 (MQ=255)
                                          cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc  <  1:751072/49‑1 (MQ=255)
                                          cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc  >  1:1891972/1‑49 (MQ=255)
                                              |                                            
TGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAC  >  NC_000913/257748‑257838

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: