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,3990CT84.1% 594.1 / 94.8 214E13K (GAG→AAG) cspCstress protein, member of the CspA family
Reads supporting (aligned to +/- strand):  ref base C (18/16);  new base T (86/94);  total (104/110)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 216 total aligned reads displayed.

CATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCC  >  NC_000913/1907351‑1907447
                                                |                                                
 aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt                                                <  2:5718543/50‑1 (MQ=255)
 aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt                                                >  2:3491803/1‑50 (MQ=255)
 aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGt                                                <  1:526153/50‑1 (MQ=255)
    tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg                                              >  1:4680830/1‑49 (MQ=255)
    tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg                                              >  2:3413513/1‑49 (MQ=255)
      tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                               <  2:4567089/46‑1 (MQ=255)
      tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa                                            <  2:389273/49‑1 (MQ=255)
      tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                           >  2:756393/1‑50 (MQ=255)
      tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                           >  1:4866345/1‑50 (MQ=255)
      tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                           >  2:3035425/1‑50 (MQ=255)
      tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                           >  1:2919723/1‑50 (MQ=255)
        ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                          <  2:1969767/49‑1 (MQ=255)
        ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                          <  2:2562585/49‑1 (MQ=255)
        ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                         >  1:2566055/1‑50 (MQ=255)
        ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                         >  1:728405/1‑50 (MQ=255)
        ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGGTGAACCa                                         <  1:2396900/50‑1 (MQ=25)
          gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                       >  1:4727654/1‑50 (MQ=255)
          gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                       >  2:114558/1‑50 (MQ=255)
           ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAGCCAc                                        <  2:4998572/48‑1 (MQ=37)
           ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                      <  2:1457705/50‑1 (MQ=255)
           ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                      <  1:5056317/50‑1 (MQ=255)
           ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                      <  1:182191/50‑1 (MQ=255)
           ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                       <  1:3714348/49‑1 (MQ=255)
           ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                       <  2:719601/49‑1 (MQ=255)
           ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACtt                                      <  2:3776372/50‑1 (MQ=255)
           ccATCAGCCAGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                       >  1:5069752/1‑49 (MQ=38)
            cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa                                     <  1:389387/50‑1 (MQ=255)
            cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa                                     <  1:4310086/50‑1 (MQ=255)
            cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa                                     <  1:3940637/50‑1 (MQ=255)
            cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa                                     <  1:3294154/50‑1 (MQ=255)
            cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACtt                                      >  2:2469410/1‑49 (MQ=255)
             aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                    <  1:1405255/50‑1 (MQ=255)
             aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                    <  1:52662/50‑1 (MQ=255)
             aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTaa                                    >  2:4239656/1‑50 (MQ=255)
                aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                                 <  1:422914/50‑1 (MQ=255)
                aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                                 >  2:3225016/1‑50 (MQ=255)
                aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                                 <  1:3573756/50‑1 (MQ=255)
                 gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                                 >  2:2576492/1‑49 (MQ=255)
                   cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                              <  1:1268865/50‑1 (MQ=255)
                   cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                              >  2:1776791/1‑50 (MQ=255)
                   cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                              <  1:2356451/50‑1 (MQ=255)
                   cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                              >  1:3499213/1‑50 (MQ=255)
                   cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                              >  1:3784949/1‑50 (MQ=255)
                    ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                             >  1:3200024/1‑50 (MQ=255)
                    ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                             >  1:5484293/1‑50 (MQ=255)
                     gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                            >  1:2327191/1‑50 (MQ=255)
                     gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                            <  1:3835778/50‑1 (MQ=255)
                      aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                           <  2:5252440/50‑1 (MQ=255)
                       gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                           >  1:1610032/1‑49 (MQ=255)
                       gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                           >  1:4373540/1‑49 (MQ=255)
                        tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                          >  2:3920614/1‑49 (MQ=255)
                        tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                          >  2:5485586/1‑49 (MQ=255)
                        tAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                         <  2:472039/50‑1 (MQ=255)
                         aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                        >  2:1399875/1‑50 (MQ=255)
                         aaTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa                        >  1:2434550/1‑50 (MQ=255)
                          aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                        <  1:609434/49‑1 (MQ=255)
                          aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt                       >  2:3638320/1‑50 (MQ=255)
                           tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                      >  2:1301571/1‑50 (MQ=255)
                           tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                      >  1:4784753/1‑50 (MQ=255)
                           tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                      <  2:631631/50‑1 (MQ=255)
                           tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                      >  1:5459331/1‑50 (MQ=255)
                            gAAGCCAAAACCTTTGGACTTGTTGAACCACTTAACCTGACCTTTAATCt                     >  2:4327760/1‑50 (MQ=25)
                             aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                    >  1:5435138/1‑50 (MQ=255)
                              aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                   <  2:1623337/50‑1 (MQ=255)
                              aGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCttt                   >  1:1077764/1‑50 (MQ=255)
                               gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg                  <  1:2206617/50‑1 (MQ=255)
                               gCCAAAACCTTTAGACTTGTTAATCCACTTAACCTGACCTTTAATCTTTg                  <  1:1903241/50‑1 (MQ=25)
                                ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc                 <  1:5798383/50‑1 (MQ=255)
                                ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc                 >  2:1727837/1‑50 (MQ=255)
                                ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc                 <  2:3207099/50‑1 (MQ=255)
                                ccAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGc                 >  1:5364842/1‑50 (MQ=255)
                                 cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGcc                <  2:2708364/50‑1 (MQ=255)
                                 cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc                 >  2:1005789/1‑49 (MQ=255)
                                  aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa               <  2:1564293/50‑1 (MQ=255)
                                  aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa               >  1:5729024/1‑50 (MQ=255)
                                   aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt              >  2:5220156/1‑50 (MQ=255)
                                   aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt              >  2:4493840/1‑50 (MQ=255)
                                   aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt              >  2:4861696/1‑50 (MQ=255)
                                   aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt              >  2:3620118/1‑50 (MQ=255)
                                    aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt             <  2:1342641/50‑1 (MQ=255)
                                    aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt             <  1:3632634/50‑1 (MQ=255)
                                     aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt            >  2:5287760/1‑50 (MQ=255)
                                      ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt            >  2:1684801/1‑49 (MQ=255)
                                      ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg           <  1:3980452/50‑1 (MQ=255)
                                       cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa          <  1:2732280/50‑1 (MQ=255)
                                       cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa          <  1:173308/50‑1 (MQ=255)
                                        tttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa          >  1:64550/1‑49 (MQ=255)
                                         ttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa         >  2:2853062/1‑49 (MQ=255)
                                         ttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa         >  2:2260868/1‑49 (MQ=255)
                                          tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa        >  2:4382785/1‑49 (MQ=255)
                                          tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa        >  1:191652/1‑49 (MQ=255)
                                          tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa        >  2:5200288/1‑49 (MQ=255)
                                          tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa         <  1:3724790/48‑1 (MQ=255)
                                          tAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa        <  2:2039338/49‑1 (MQ=255)
                                           aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa      >  2:868283/1‑50 (MQ=255)
                                           aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa      >  1:1870253/1‑50 (MQ=255)
                                            gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa      <  2:3092790/49‑1 (MQ=255)
                                            gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt     <  2:3269078/50‑1 (MQ=255)
                                            gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt     >  1:843134/1‑50 (MQ=255)
                                            gACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt     <  1:93587/50‑1 (MQ=255)
                                                |                                                
CATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCC  >  NC_000913/1907351‑1907447

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