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,7940TC88.1% 424.2 / 40.5 134intergenic (+23/‑35)frsA/crlfermentation‑respiration switch protein/RNA polymerase holoenzyme assembly factor Crl
Reads supporting (aligned to +/- strand):  ref base T (12/4);  new base C (67/51);  total (79/55)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 134 total aligned reads displayed.

CTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCG  >  NC_000913/257747‑257840
                                               |                                              
cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                              <  2:2412239/50‑1 (MQ=255)
cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                              <  1:2240914/50‑1 (MQ=255)
 tGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                             >  2:2664295/1‑50 (MQ=255)
 tGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                             >  2:1934966/1‑50 (MQ=255)
  ggATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                             >  2:4083897/1‑49 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTgg                                             <  2:3118456/48‑1 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGt                                            >  2:4621530/1‑49 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                             <  1:4181815/48‑1 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                             <  2:5040618/48‑1 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                            <  2:1270868/49‑1 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                           <  2:4092134/50‑1 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                           <  1:4271088/50‑1 (MQ=255)
   gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                           <  2:3541211/50‑1 (MQ=255)
    aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTa                                           >  1:4883510/1‑49 (MQ=255)
    aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTa                                           >  1:4952827/1‑49 (MQ=255)
     tCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaa                                          >  1:974352/1‑49 (MQ=255)
      cGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  2:2748192/1‑50 (MQ=255)
       gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAAc                                       <  1:718507/50‑1 (MQ=255)
        aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  2:1365618/1‑48 (MQ=255)
        aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  2:2355571/1‑48 (MQ=255)
        aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACa                                      <  2:567216/50‑1 (MQ=255)
         aaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAg                                     <  1:1895536/50‑1 (MQ=255)
          aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                    <  2:3327969/50‑1 (MQ=255)
           aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGtt                                   <  1:996297/50‑1 (MQ=255)
           aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGtt                                   >  2:1422620/1‑50 (MQ=255)
           aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGtt                                   >  1:4740812/1‑50 (MQ=255)
              gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                                >  2:1938910/1‑50 (MQ=255)
              gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                                >  1:637790/1‑50 (MQ=255)
               cTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                               >  2:580462/1‑50 (MQ=255)
               cTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                               >  1:519063/1‑50 (MQ=255)
               cTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                               <  1:1301060/50‑1 (MQ=255)
                ttGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                               >  1:4645126/1‑49 (MQ=255)
                 tgtgTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                             >  2:2844826/1‑50 (MQ=255)
                  gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                             <  1:896644/49‑1 (MQ=255)
                  gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                             >  1:4688996/1‑49 (MQ=255)
                  gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                             <  2:4879472/49‑1 (MQ=255)
                  gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATc                              <  2:242394/48‑1 (MQ=255)
                   tgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATcaca                           <  1:2861042/50‑1 (MQ=255)
                   tgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATcac                            >  2:2420643/1‑49 (MQ=255)
                   tgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATcac                            >  2:5172952/1‑49 (MQ=255)
                   tgtTAAAAATTGGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATcaca                           <  2:4736221/50‑1 (MQ=255)
                   tgtTAAAAAGTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATcac                            <  1:200242/49‑1 (MQ=38)
                    gtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                          >  1:1103041/1‑50 (MQ=255)
                    gtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                          >  1:3179719/1‑50 (MQ=255)
                     ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                          >  1:1545373/1‑49 (MQ=255)
                      tAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                          >  2:2436780/1‑48 (MQ=255)
                       aaaaaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAg                       >  1:5054486/1‑50 (MQ=255)
                       aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                       >  2:890402/1‑50 (MQ=255)
                       aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                       >  1:2509462/1‑50 (MQ=255)
                       aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                       >  1:3336475/1‑50 (MQ=255)
                        aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAgg                      <  1:4576986/50‑1 (MQ=255)
                        aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAgg                      <  1:10620/50‑1 (MQ=255)
                        aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                       >  2:1937037/1‑49 (MQ=255)
                         aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGga                     <  2:1125659/50‑1 (MQ=255)
                          aaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgag                    <  2:3854884/50‑1 (MQ=255)
                           aTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                       >  2:440805/1‑46 (MQ=255)
                           aTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                       >  1:3881634/1‑46 (MQ=255)
                            tttGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGgaga                   >  1:4993986/1‑49 (MQ=255)
                            tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGtaga                   >  1:2453679/1‑49 (MQ=38)
                            tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGAt                  <  1:283207/50‑1 (MQ=255)
                             ttGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGgaga                   >  2:3638597/1‑48 (MQ=255)
                             ttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                   >  1:1618376/1‑48 (MQ=255)
                              tGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATa                 <  2:3277788/49‑1 (MQ=255)
                              tGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATa                 <  2:5127428/49‑1 (MQ=255)
                              tGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg                >  1:724072/1‑50 (MQ=255)
                               gCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg                <  1:3693494/49‑1 (MQ=255)
                               gCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg                <  1:1117162/49‑1 (MQ=255)
                               gCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg                <  2:5327642/49‑1 (MQ=255)
                                cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGc               >  2:2302680/1‑49 (MQ=255)
                                cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGc               <  2:3096655/49‑1 (MQ=255)
                                cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCa              <  1:1903271/50‑1 (MQ=255)
                                 tAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCaa             >  1:5263914/1‑50 (MQ=255)
                                 tAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCaa             >  2:1392490/1‑50 (MQ=255)
                                   aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCaa             <  2:3060323/48‑1 (MQ=255)
                                   aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt            >  1:3963846/1‑49 (MQ=255)
                                   aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt            >  2:2596812/1‑49 (MQ=255)
                                   aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt            <  1:4185013/49‑1 (MQ=255)
                                   aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt            <  1:3741904/49‑1 (MQ=255)
                                   aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt            <  1:3440026/49‑1 (MQ=255)
                                   aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt            <  2:2239988/49‑1 (MQ=255)
                                    aTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGa          >  1:1962744/1‑50 (MQ=255)
                                     ttttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc         >  1:3026066/1‑50 (MQ=255)
                                     ttttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc         >  2:1467870/1‑50 (MQ=255)
                                     ttttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc         <  2:4874524/50‑1 (MQ=255)
                                     ttttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc         <  2:4866554/50‑1 (MQ=255)
                                      tttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACg        >  1:1372279/1‑50 (MQ=255)
                                      tttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACg        >  1:169835/1‑50 (MQ=255)
                                       ttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACg        <  1:3025281/49‑1 (MQ=255)
                                        tGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt       >  1:2458088/1‑49 (MQ=255)
                                        tGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt       >  2:80535/1‑49 (MQ=255)
                                        tGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt       >  1:3926861/1‑49 (MQ=255)
                                        tGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt       >  1:3644075/1‑49 (MQ=255)
                                         gCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTa     >  1:643542/1‑50 (MQ=255)
                                         gCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTa     >  2:4711486/1‑50 (MQ=255)
                                          ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc    >  1:3090427/1‑50 (MQ=255)
                                           cAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   <  1:3667621/50‑1 (MQ=255)
                                           cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   <  2:3832741/50‑1 (MQ=255)
                                           cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   >  2:2683615/1‑50 (MQ=255)
                                           cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   <  1:3049325/50‑1 (MQ=255)
                                            aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCg  <  1:4067406/50‑1 (MQ=255)
                                               |                                              
CTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCG  >  NC_000913/257747‑257840

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