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,7940TC98.9% 352.4 / ‑6.4 87intergenic (+23/‑35)frsA/crlfermentation‑respiration switch protein/RNA polymerase holoenzyme assembly factor Crl
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (59/27);  minor base A (1/0);  total (60/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACC  >  NC_000913/257747‑257839
                                               |                                             
cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                             >  1:5042763/1‑50 (MQ=255)
cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                             >  1:2760875/1‑50 (MQ=255)
cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                             >  1:2697414/1‑50 (MQ=255)
 tGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                            >  1:616405/1‑50 (MQ=255)
 tGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                            >  1:4426587/1‑50 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                            <  1:4414032/48‑1 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                           <  1:1481585/49‑1 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                           >  1:3251464/1‑49 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                          <  1:3310662/50‑1 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAACATGGt                                           >  1:5177642/1‑49 (MQ=38)
    aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaa                                         <  1:5724987/50‑1 (MQ=255)
    aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaa                                         <  1:6384230/50‑1 (MQ=255)
    aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                          >  1:2745018/1‑49 (MQ=255)
    aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                          >  1:5099904/1‑49 (MQ=255)
     tCGAAAAACGCTTGTGTTAAAAATTTGCTAACTTTTGCCAATCTGGTaa                                         >  1:4770250/1‑49 (MQ=38)
       gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:5334901/1‑49 (MQ=255)
        aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:4241256/1‑48 (MQ=255)
        aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:103443/1‑48 (MQ=255)
        aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:907027/1‑48 (MQ=255)
        aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:922780/1‑48 (MQ=255)
        aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:1700049/1‑48 (MQ=255)
          aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                   >  1:1265731/1‑50 (MQ=255)
          aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                   <  1:4823642/50‑1 (MQ=255)
          aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                   >  1:1526278/1‑50 (MQ=255)
           aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGtt                                  >  1:1022761/1‑50 (MQ=255)
           aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGtt                                  <  1:504379/50‑1 (MQ=255)
            aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                 <  1:5184739/50‑1 (MQ=255)
             cGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                 >  1:761157/1‑49 (MQ=255)
              gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                               >  1:935820/1‑50 (MQ=255)
              gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                               >  1:3698941/1‑50 (MQ=255)
              gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                               <  1:87163/50‑1 (MQ=255)
              gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                               >  1:3248769/1‑50 (MQ=255)
              gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                               >  1:5209826/1‑50 (MQ=255)
              gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                               >  1:3095466/1‑50 (MQ=255)
               cTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:5691315/1‑50 (MQ=255)
               cTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              <  1:6332250/50‑1 (MQ=255)
               cTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:2179714/1‑50 (MQ=255)
                ttGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:1130348/1‑49 (MQ=255)
                ttGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:872013/1‑49 (MQ=255)
                ttGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:454642/1‑49 (MQ=255)
                ttGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:373172/1‑49 (MQ=255)
                  gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                            >  1:769530/1‑49 (MQ=255)
                  gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                            <  1:4317622/49‑1 (MQ=255)
                  gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                            >  1:7211794/1‑49 (MQ=255)
                  gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                            <  1:5743951/49‑1 (MQ=255)
                  gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                            <  1:4863234/49‑1 (MQ=255)
                     ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                         <  1:389757/49‑1 (MQ=255)
                     ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                         >  1:5111065/1‑49 (MQ=255)
                      tAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                         >  1:5443261/1‑48 (MQ=255)
                      tAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                         >  1:7377454/1‑48 (MQ=255)
                       aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaaca                       <  1:4787566/49‑1 (MQ=255)
                       aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaaca                       <  1:6963032/49‑1 (MQ=255)
                       aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:3597501/1‑50 (MQ=255)
                        aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAgg                     <  1:1994300/50‑1 (MQ=255)
                         aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:3371866/1‑48 (MQ=255)
                         aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:322185/1‑48 (MQ=255)
                         aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGga                    <  1:5978628/50‑1 (MQ=255)
                          aaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:5414663/1‑47 (MQ=255)
                           aTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:1609244/1‑46 (MQ=255)
                           aTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:4108043/1‑46 (MQ=255)
                            tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                  >  1:3426609/1‑49 (MQ=255)
                            tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGAt                 <  1:2782490/50‑1 (MQ=255)
                             ttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                  >  1:3268754/1‑48 (MQ=255)
                             ttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                  >  1:683869/1‑48 (MQ=255)
                              tGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATa                <  1:1099502/49‑1 (MQ=255)
                               gCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg               >  1:7446681/1‑49 (MQ=255)
                               gCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg               >  1:508963/1‑49 (MQ=255)
                                cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGc              <  1:4908347/49‑1 (MQ=255)
                                cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGc              >  1:308964/1‑49 (MQ=255)
                                cTAAATCTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGc              >  1:160105/1‑49 (MQ=38)
                                   aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCaa            <  1:3897460/48‑1 (MQ=255)
                                   aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt           >  1:6063720/1‑49 (MQ=255)
                                    aTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGa         <  1:2652836/50‑1 (MQ=255)
                                     ttttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc        <  1:1726082/50‑1 (MQ=255)
                                     ttttGCCAAGCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc        <  1:2900636/50‑1 (MQ=255)
                                      tttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACg       >  1:6256306/1‑50 (MQ=255)
                                      tttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACg       >  1:1688424/1‑50 (MQ=255)
                                        tGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt      >  1:4002599/1‑49 (MQ=255)
                                        tGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt      >  1:3307281/1‑49 (MQ=255)
                                         gCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTa    >  1:6535055/1‑50 (MQ=255)
                                          ccAATCTGGTAAGACAGTTGCATCACAACAGGAGATAGCAATGACGTTa    <  1:5445708/49‑1 (MQ=255)
                                          ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc   <  1:264627/50‑1 (MQ=255)
                                          ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc   >  1:3968778/1‑50 (MQ=255)
                                          ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc   >  1:2359282/1‑50 (MQ=255)
                                           cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc  >  1:5605799/1‑50 (MQ=255)
                                            aaTCTGGTAAAACAGTTGCATCACAACAGGAGATATCAATGACGTTAcc  >  1:2780192/1‑49 (MQ=38)
                                            aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc  >  1:911384/1‑49 (MQ=255)
                                               |                                             
CTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACC  >  NC_000913/257747‑257839

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