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,3990CT97.5% 618.3 / 3.3 160E13K (GAG→AAG) cspCstress protein, member of the CspA family
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (82/74);  total (85/75)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

Only 100 of 164 total aligned reads displayed.

ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT  >  NC_000913/1907352‑1907445
                                               |                                              
aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt                                              <  1:5209960/50‑1 (MQ=255)
aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt                                              <  1:4170583/50‑1 (MQ=255)
aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt                                              <  1:3421949/50‑1 (MQ=255)
  cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                             >  2:5038494/1‑49 (MQ=255)
   tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg                                            >  2:4949528/1‑49 (MQ=255)
   tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg                                            >  2:2262656/1‑49 (MQ=255)
    ttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa                                          >  1:198452/1‑50 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg                                            <  2:3096324/47‑1 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                             <  2:1614450/46‑1 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGa                                           <  2:5329250/48‑1 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGa                                           <  1:2704388/48‑1 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGa                                           <  2:1715653/48‑1 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                         >  2:4463978/1‑50 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                         >  2:4335750/1‑50 (MQ=255)
      gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                         >  1:4618913/1‑49 (MQ=255)
        tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      >  2:616892/1‑50 (MQ=255)
        tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      >  1:5357409/1‑50 (MQ=255)
        tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      >  1:3365831/1‑50 (MQ=255)
         gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                     >  1:2729549/1‑50 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      <  1:2831582/48‑1 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                    <  1:4462329/50‑1 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                     <  1:5293757/49‑1 (MQ=255)
           cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa                                   <  1:3197695/50‑1 (MQ=255)
            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                  >  1:4091427/1‑50 (MQ=255)
            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                  >  1:1691451/1‑50 (MQ=255)
               aGCCGGAGTAATGAAGCCACAACCTTTAGACTTGTTGAACCACTTAACCt                               >  2:1186650/1‑50 (MQ=37)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               <  2:2062563/50‑1 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:706196/1‑50 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:1929132/1‑50 (MQ=255)
                gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTg                              <  2:1212808/50‑1 (MQ=255)
                 ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                             >  2:2708231/1‑50 (MQ=255)
                  cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                            <  2:1760151/50‑1 (MQ=255)
                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           >  1:1024699/1‑50 (MQ=255)
                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           >  2:381071/1‑50 (MQ=255)
                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           <  1:209997/50‑1 (MQ=255)
                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          >  1:2907163/1‑50 (MQ=255)
                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          >  1:4553450/1‑50 (MQ=255)
                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         >  1:3579521/1‑50 (MQ=255)
                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         >  1:3957936/1‑50 (MQ=255)
                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         <  1:4811346/50‑1 (MQ=255)
                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         >  2:4756700/1‑50 (MQ=255)
                      gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        <  1:2720744/50‑1 (MQ=255)
                      gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        <  1:5518318/50‑1 (MQ=255)
                      gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        <  2:5229378/50‑1 (MQ=255)
                       tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        >  1:4371589/1‑49 (MQ=255)
                       tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        >  1:2356016/1‑49 (MQ=255)
                       tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTa                       <  2:4005360/50‑1 (MQ=255)
                        aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      >  1:1087657/1‑50 (MQ=255)
                        aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      >  1:5032911/1‑50 (MQ=255)
                        aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      >  1:4986272/1‑50 (MQ=255)
                        aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      >  1:1916397/1‑50 (MQ=255)
                        aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      >  1:4606742/1‑50 (MQ=255)
                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      <  2:1361302/49‑1 (MQ=255)
                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt                     >  2:4630912/1‑50 (MQ=255)
                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt                     <  1:3069632/50‑1 (MQ=255)
                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt                     >  2:3591034/1‑50 (MQ=255)
                          tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                    >  2:1583094/1‑50 (MQ=255)
                          tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                    >  1:4635416/1‑50 (MQ=255)
                           gAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                   >  1:1767113/1‑50 (MQ=255)
                           gAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                   >  2:297853/1‑50 (MQ=255)
                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                  >  2:3852957/1‑50 (MQ=255)
                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                  <  1:3638453/50‑1 (MQ=255)
                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                  >  2:1264292/1‑50 (MQ=255)
                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                  >  2:2486563/1‑50 (MQ=255)
                             aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                 <  1:991141/50‑1 (MQ=255)
                             aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                 >  2:2858014/1‑50 (MQ=255)
                              gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg                <  2:219010/50‑1 (MQ=255)
                              gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg                <  1:250520/50‑1 (MQ=255)
                               ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               >  1:76521/1‑50 (MQ=255)
                                cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               >  1:5406311/1‑49 (MQ=255)
                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             >  1:2119326/1‑50 (MQ=255)
                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             <  2:745811/50‑1 (MQ=255)
                                  aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt            >  2:195687/1‑50 (MQ=255)
                                  aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt            >  1:846691/1‑50 (MQ=255)
                                   aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt           >  1:1844087/1‑50 (MQ=255)
                                    aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          <  1:4484784/50‑1 (MQ=255)
                                    aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          <  1:1534623/50‑1 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          >  2:4263132/1‑49 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          >  2:4054704/1‑49 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg         <  1:2829395/50‑1 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTAAATCTTTGCCATTTg         <  1:3098813/50‑1 (MQ=38)
                                      cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg         >  1:5566972/1‑49 (MQ=255)
                                      cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg         >  2:1512092/1‑49 (MQ=255)
                                       tttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa        >  2:2478122/1‑49 (MQ=255)
                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGac       <  2:1047918/48‑2 (MQ=255)
                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa      >  1:2339327/1‑49 (MQ=255)
                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa       <  1:610302/48‑1 (MQ=255)
                                          aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa    >  1:2439846/1‑50 (MQ=255)
                                          aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa    >  1:1158572/1‑50 (MQ=255)
                                          aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa    >  1:1960303/1‑50 (MQ=255)
                                          aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa    >  1:1691921/1‑50 (MQ=255)
                                          aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaa     <  1:5187786/49‑1 (MQ=255)
                                          aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaa     <  1:3018589/49‑1 (MQ=255)
                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa    <  2:2570490/49‑1 (MQ=255)
                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa    <  2:1565994/49‑1 (MQ=255)
                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt   <  1:1349436/50‑1 (MQ=255)
                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt   <  2:2607275/50‑1 (MQ=255)
                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt   >  2:913525/1‑50 (MQ=255)
                                            aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt  >  2:1454880/1‑50 (MQ=255)
                                            aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt  <  1:2941081/50‑1 (MQ=255)
                                               |                                              
ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT  >  NC_000913/1907352‑1907445

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