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,3990CT99.1% 414.5 / ‑4.0 107E13K (GAG→AAG) cspCstress protein, member of the CspA family
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (54/52);  total (55/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

Only 100 of 109 total aligned reads displayed.

ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT  >  NC_000913/1907352‑1907445
                                               |                                              
aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt                                              >  1:756562/1‑50 (MQ=255)
aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt                                              <  1:6247446/50‑1 (MQ=255)
 tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                             >  1:162552/1‑50 (MQ=255)
  cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                             >  1:2174188/1‑49 (MQ=255)
  cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                             >  1:5702148/1‑49 (MQ=255)
   tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg                                            >  1:5338025/1‑49 (MQ=255)
    ttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa                                          >  1:6040330/1‑50 (MQ=255)
    ttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGaa                                          >  1:6253498/1‑50 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg                                            <  1:2018014/47‑1 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtg                                             <  1:4230664/46‑2 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa                                          <  1:1668825/49‑1 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGa                                           <  1:3401997/48‑1 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGa                                           <  1:6330402/48‑1 (MQ=255)
      gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                         <  1:6808688/49‑1 (MQ=255)
      gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                         >  1:2401607/1‑49 (MQ=255)
       ttgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                       >  1:2324667/2‑50 (MQ=255)
       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                        <  1:6073522/49‑1 (MQ=255)
       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                        <  1:1511584/49‑1 (MQ=255)
       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                       >  1:4528000/1‑50 (MQ=255)
       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                       <  1:3943090/50‑1 (MQ=255)
        tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      >  1:7965906/1‑50 (MQ=255)
        tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      >  1:4330316/1‑50 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      <  1:7071317/48‑1 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      <  1:2631300/48‑1 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      <  1:3843772/48‑1 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                    <  1:2106475/50‑1 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                     >  1:3090869/1‑49 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                     <  1:7426640/49‑1 (MQ=255)
           cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                    >  1:640354/1‑49 (MQ=255)
           cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa                                   <  1:6100587/50‑1 (MQ=255)
           cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa                                   <  1:89049/50‑1 (MQ=255)
            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCGCTTaa                                  <  1:1802053/50‑1 (MQ=37)
            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                  >  1:5788743/1‑50 (MQ=255)
            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                  >  1:4123196/1‑50 (MQ=255)
             tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAc                                 >  1:2340231/1‑50 (MQ=255)
              cAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                                >  1:4795291/1‑50 (MQ=255)
              cAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                                >  1:2780816/1‑50 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               <  1:7218030/50‑1 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               <  1:6346235/50‑1 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:5491776/1‑50 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:4419589/1‑50 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               <  1:8084106/50‑1 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:625144/1‑50 (MQ=255)
                gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:2145829/1‑49 (MQ=255)
                gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTg                              <  1:579003/50‑1 (MQ=255)
                 ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                             <  1:1938462/50‑1 (MQ=255)
                 ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                             >  1:113063/1‑50 (MQ=255)
                  cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                            <  1:1438295/50‑1 (MQ=255)
                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           <  1:3906135/50‑1 (MQ=255)
                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           >  1:7836013/1‑50 (MQ=255)
                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           >  1:2584567/1‑50 (MQ=255)
                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           <  1:5478160/50‑1 (MQ=255)
                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          <  1:1839434/50‑1 (MQ=255)
                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          <  1:4572283/50‑1 (MQ=255)
                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          >  1:6910534/1‑50 (MQ=255)
                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         >  1:2316449/1‑50 (MQ=255)
                       tAATGAAGCCGAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        >  1:4739766/1‑49 (MQ=25)
                       tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTa                       <  1:1221520/50‑1 (MQ=255)
                        aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      >  1:4890350/1‑50 (MQ=255)
                        aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      >  1:2569434/1‑50 (MQ=255)
                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      <  1:4863296/49‑1 (MQ=255)
                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      <  1:3673664/49‑1 (MQ=255)
                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt                     >  1:7472186/1‑50 (MQ=255)
                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt                     >  1:3811466/1‑50 (MQ=255)
                          tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                    <  1:2825903/50‑1 (MQ=255)
                           gAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                   >  1:7291584/1‑50 (MQ=255)
                           gAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                   >  1:8163518/1‑50 (MQ=255)
                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                  <  1:2732046/50‑1 (MQ=255)
                             aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                 >  1:2264497/1‑50 (MQ=255)
                             aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                 <  1:2549765/50‑1 (MQ=255)
                              gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg                <  1:8084341/50‑1 (MQ=255)
                              gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg                <  1:1167577/50‑1 (MQ=255)
                              gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg                <  1:3641612/50‑1 (MQ=255)
                               ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               <  1:559393/50‑1 (MQ=255)
                               ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               >  1:2226294/1‑50 (MQ=255)
                                cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               >  1:6226797/1‑49 (MQ=255)
                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             <  1:8533601/50‑1 (MQ=255)
                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             <  1:1080085/50‑1 (MQ=255)
                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             >  1:6732387/1‑50 (MQ=255)
                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             >  1:1025838/1‑50 (MQ=255)
                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             >  1:5302508/1‑50 (MQ=255)
                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             >  1:2118285/1‑50 (MQ=255)
                                  aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt            >  1:4218909/1‑50 (MQ=255)
                                   aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt           <  1:5395696/50‑1 (MQ=255)
                                    aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          >  1:820485/1‑50 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          >  1:6646220/1‑49 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg         <  1:4643541/50‑1 (MQ=255)
                                      cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg         >  1:4526110/1‑49 (MQ=255)
                                      cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg         >  1:5195191/1‑49 (MQ=255)
                                      cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa        <  1:3154591/50‑1 (MQ=255)
                                       tttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa        >  1:513106/1‑49 (MQ=255)
                                       tttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa        >  1:3842224/1‑49 (MQ=255)
                                       tttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa        >  1:1495988/1‑49 (MQ=255)
                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa      <  1:7692073/49‑1 (MQ=255)
                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa       <  1:4776482/48‑1 (MQ=255)
                                         tAGACTTGTTGAACCACTTAACCTGAACTTTAATCTTTGCCATTTGaaa      >  1:7444155/1‑49 (MQ=38)
                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa    <  1:3719156/49‑1 (MQ=255)
                                            aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt  >  1:7867127/1‑50 (MQ=255)
                                            aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt  <  1:6391629/50‑1 (MQ=255)
                                            aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt  >  1:6625642/1‑50 (MQ=255)
                                               |                                              
ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT  >  NC_000913/1907352‑1907445

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