Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,907,399 C→T E13K (GAG→AAG)  cspC ← stress protein, member of the CspA family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,907,3990CT98.0% 195.7 / ‑3.9 51E13K (GAG→AAG) cspCstress protein, member of the CspA family
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (22/28);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAT  >  NC_000913/1907352‑1907444
                                               |                                             
aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt                                             >  1:1977663/1‑50 (MQ=255)
 tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                            >  1:8740/1‑50 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                            <  1:2606723/46‑1 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa                                         <  1:1536179/49‑1 (MQ=255)
       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                      >  1:1440202/1‑50 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                     <  1:2971133/48‑1 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                    >  1:2703121/1‑49 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                    <  1:506797/49‑1 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACt                                    <  1:215366/49‑1 (MQ=255)
            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                 >  1:66476/1‑50 (MQ=255)
            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                 >  1:2161081/1‑50 (MQ=255)
            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                 <  1:3614868/50‑1 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                               <  1:1067705/49‑1 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                               <  1:2167537/49‑1 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                              <  1:270511/50‑1 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                              >  1:3656879/1‑50 (MQ=255)
                 ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                            <  1:2769586/50‑1 (MQ=255)
                 ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                            >  1:2852887/1‑50 (MQ=255)
                 ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                            >  1:3052800/1‑50 (MQ=255)
                 ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                            >  1:1521743/1‑50 (MQ=255)
                  cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                           >  1:25667/1‑50 (MQ=255)
                  cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                           >  1:2627735/1‑50 (MQ=255)
                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                          <  1:626654/50‑1 (MQ=255)
                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                          <  1:1274282/50‑1 (MQ=255)
                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                         <  1:3519735/50‑1 (MQ=255)
                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                        >  1:167890/1‑50 (MQ=255)
                      gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                       <  1:3371732/50‑1 (MQ=255)
                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt                    <  1:3036752/50‑1 (MQ=255)
                          tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                   <  1:3259142/50‑1 (MQ=255)
                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                 >  1:941270/1‑50 (MQ=255)
                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                 <  1:2653183/50‑1 (MQ=255)
                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                  <  1:1185959/49‑1 (MQ=255)
                             aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                >  1:2460378/1‑50 (MQ=255)
                             aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                >  1:1525338/1‑50 (MQ=255)
                              gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg               <  1:1374577/50‑1 (MQ=255)
                               ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc              <  1:3932144/50‑1 (MQ=255)
                                cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGcc             <  1:12652/50‑1 (MQ=255)
                                  aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt           <  1:299061/50‑1 (MQ=255)
                                    aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt         <  1:2558996/50‑1 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt         >  1:3161783/1‑49 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg        <  1:3022770/50‑1 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg        <  1:2286023/50‑1 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg        <  1:1985952/50‑1 (MQ=255)
                                      cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg        >  1:280625/1‑49 (MQ=255)
                                       tttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa       >  1:76057/1‑49 (MQ=255)
                                        ttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa      >  1:2803567/1‑49 (MQ=255)
                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa     <  1:2305855/49‑1 (MQ=255)
                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa     >  1:2349149/1‑49 (MQ=255)
                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa      <  1:2566297/48‑1 (MQ=255)
                                           gACTTGTTGAACCACTTACCCTGACCTTTAATCTTTGCCATTTGAAAAAt  <  1:2901923/50‑1 (MQ=38)
                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt  >  1:1209025/1‑50 (MQ=255)
                                               |                                             
ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAT  >  NC_000913/1907352‑1907444

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