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,3990CT94.0% 296.1 / 7.0 83E13K (GAG→AAG) cspCstress protein, member of the CspA family
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (30/48);  total (32/51)
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.98e-01

ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAT  >  NC_000913/1907352‑1907444
                                               |                                             
aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGt                                             >  1:1339143/1‑50 (MQ=255)
 tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                            >  1:2890803/1‑50 (MQ=255)
  cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                            >  1:2677379/1‑49 (MQ=255)
   tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg                                           >  1:2741747/1‑49 (MQ=255)
    ttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa                                         >  1:907295/1‑50 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa                                         <  1:1212751/49‑1 (MQ=255)
     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGa                                          <  1:2259196/48‑1 (MQ=255)
      gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                        >  1:3251355/1‑49 (MQ=255)
      gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                        >  1:2577490/1‑49 (MQ=255)
       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                       <  1:1484912/49‑1 (MQ=255)
       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                       <  1:4393350/49‑1 (MQ=255)
       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                      >  1:287370/1‑50 (MQ=255)
       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCa                                      <  1:200963/50‑1 (MQ=255)
         gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                    >  1:3876705/1‑50 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                     <  1:1761507/48‑1 (MQ=255)
          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                    <  1:792710/49‑1 (MQ=255)
          ccATCAGCCGGAGTAATGAAACCAAAACCTTTAGACTTGTTGAACCACtt                                   <  1:4299107/50‑1 (MQ=25)
           cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa                                  <  1:1075399/50‑1 (MQ=255)
            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                 <  1:1447526/50‑1 (MQ=255)
            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                 <  1:3010420/50‑1 (MQ=255)
             tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAc                                <  1:1061383/50‑1 (MQ=255)
             tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAAc                                <  1:360547/50‑1 (MQ=255)
              cAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                               >  1:2312437/1‑50 (MQ=255)
              cAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                               >  1:1566352/1‑50 (MQ=255)
               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                              <  1:3338402/50‑1 (MQ=255)
                gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                              >  1:1340953/1‑49 (MQ=255)
                gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTg                             <  1:828389/50‑1 (MQ=255)
                 ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                            >  1:3242608/1‑50 (MQ=255)
                  cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                           >  1:3341131/1‑50 (MQ=255)
                  cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                           >  1:3507779/1‑50 (MQ=255)
                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                          >  1:2072664/1‑50 (MQ=255)
                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                          >  1:3031417/1‑50 (MQ=255)
                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                          >  1:961561/1‑50 (MQ=255)
                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                         >  1:1760363/1‑50 (MQ=255)
                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                        <  1:4237838/50‑1 (MQ=255)
                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                        <  1:3942406/50‑1 (MQ=255)
                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                        <  1:3878337/50‑1 (MQ=255)
                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                        <  1:3140107/50‑1 (MQ=255)
                      gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                       <  1:450368/50‑1 (MQ=255)
                      gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                        >  1:2911666/1‑49 (MQ=255)
                       tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                       >  1:2366306/1‑49 (MQ=255)
                       tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTa                      <  1:3964633/50‑1 (MQ=255)
                       tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTa                      <  1:1390306/50‑1 (MQ=255)
                       tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTa                      <  1:2494616/50‑1 (MQ=255)
                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                     <  1:971441/49‑1 (MQ=255)
                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                     <  1:4324172/49‑1 (MQ=255)
                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt                    >  1:439508/1‑50 (MQ=255)
                          tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                   >  1:3993033/1‑50 (MQ=255)
                          tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                   <  1:1505479/50‑1 (MQ=255)
                          tGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATc                   <  1:43347/50‑1 (MQ=255)
                          tGAAACCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                   <  1:1746535/50‑1 (MQ=25)
                           gAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                  >  1:3121498/1‑50 (MQ=255)
                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                 <  1:4367012/50‑1 (MQ=255)
                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                  <  1:756871/49‑1 (MQ=255)
                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                  <  1:3701069/49‑1 (MQ=255)
                             aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                <  1:669328/50‑1 (MQ=255)
                             aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                <  1:302967/50‑1 (MQ=255)
                             aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                <  1:2713667/50‑1 (MQ=255)
                               ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc              <  1:2611294/50‑1 (MQ=255)
                                cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGcc             <  1:4184220/50‑1 (MQ=255)
                                cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc              >  1:950778/1‑49 (MQ=255)
                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa            <  1:642130/50‑1 (MQ=255)
                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa            >  1:20096/1‑50 (MQ=255)
                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa            <  1:3759523/50‑1 (MQ=255)
                                  aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt           >  1:3499097/1‑50 (MQ=255)
                                  aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt           <  1:1081943/50‑1 (MQ=255)
                                  aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt           >  1:3086892/1‑50 (MQ=255)
                                  aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt           >  1:1444179/1‑50 (MQ=255)
                                  aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt           <  1:2704058/50‑1 (MQ=255)
                                   aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt          <  1:2615977/50‑1 (MQ=255)
                                   aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt          <  1:497527/50‑1 (MQ=255)
                                   aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt          <  1:4033996/50‑1 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg        <  1:4016154/50‑1 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg        <  1:596191/50‑1 (MQ=255)
                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg        <  1:653257/50‑1 (MQ=255)
                                      cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg        >  1:4142620/1‑49 (MQ=255)
                                       tttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa       >  1:1073520/1‑49 (MQ=255)
                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa     <  1:4021133/49‑1 (MQ=255)
                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa     >  1:300970/1‑49 (MQ=255)
                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa      <  1:3184869/48‑1 (MQ=255)
                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa   <  1:4405525/49‑1 (MQ=255)
                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt  <  1:918675/50‑1 (MQ=255)
                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt  <  1:3670126/50‑1 (MQ=255)
                                           gACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt  >  1:3625107/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: