Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 980,263 A→T 5.8% P14P (CCT→CCA) ‡ cspI ← cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE980,2630AT 5.8% 337.4 / 8.8 119P14P (CCT→CCA) ‡cspIcold shock protein
Reads supporting (aligned to +/- strand):  ref base A (61/51);  new base T (4/3);  total (65/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

AATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCT  >  minE/980194‑980333
                                                                     |                                                                      
aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                                       <  1:4446119/71‑1 (MQ=255)
  ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                                        <  1:5582488/68‑1 (MQ=255)
  ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                                        <  1:4550533/68‑1 (MQ=255)
  ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                                        <  1:4578964/68‑1 (MQ=255)
    gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                                                    >  1:3491141/1‑70 (MQ=255)
    gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                                                    >  1:3022718/1‑70 (MQ=255)
     cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa                                                                  <  1:3922659/71‑1 (MQ=255)
     cTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                                        >  1:2389092/1‑65 (MQ=255)
       gagaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaaa                                                                >  1:1216871/1‑71 (MQ=255)
          aaaTGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCa                                                             <  1:2167906/71‑1 (MQ=255)
          aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                                               >  1:220450/1‑69 (MQ=255)
            atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                                           >  1:1039078/4‑71 (MQ=255)
            atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                                           >  1:1490823/4‑71 (MQ=255)
               gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                                       >  1:2381801/1‑56 (MQ=255)
               caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCAtttta                                                        >  1:1453682/2‑69 (MQ=255)
               caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACAAtttta                                                        >  1:1298922/2‑69 (MQ=255)
                aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                                        <  1:395464/54‑1 (MQ=255)
                aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                                        <  1:4158835/54‑1 (MQ=255)
                  aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                                     <  1:720516/71‑1 (MQ=255)
                  aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                                     <  1:5510559/71‑1 (MQ=255)
                  aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                                     <  1:2315796/71‑1 (MQ=255)
                  aaaCACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACt                                                      >  1:5153911/1‑70 (MQ=255)
                    acacATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTAcga                                                     >  1:1052620/1‑67 (MQ=9)
                     cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                                                          >  1:1567424/1‑63 (MQ=255)
                     cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                                                  <  1:1513145/71‑1 (MQ=255)
                       caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                                       >  1:5185732/1‑48 (MQ=255)
                       caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                                       >  1:1659681/1‑48 (MQ=255)
                          cTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCa                                                >  1:5511648/1‑68 (MQ=255)
                           tttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                                               <  1:5616151/52‑1 (MQ=255)
                            ttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGt                                              >  1:3644293/1‑68 (MQ=255)
                            ttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAt                                           >  1:5221626/1‑71 (MQ=255)
                             tgctgcCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAAc                                                  >  1:4022803/1‑63 (MQ=9)
                              gctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtt                                          <  1:4479250/70‑1 (MQ=255)
                              gctgcCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCAtt                                          <  1:2100745/70‑1 (MQ=12)
                               gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt                                         <  1:3892167/68‑1 (MQ=255)
                               gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt                                         <  1:2003423/68‑1 (MQ=255)
                               ctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAttt                                         <  1:1591862/70‑1 (MQ=255)
                                  ccATCTTTGGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                      >  1:2737445/1‑70 (MQ=255)
                                  ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                      <  1:4351269/70‑1 (MQ=255)
                                    aTCTTTCGGTGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                                     <  1:3328034/53‑1 (MQ=255)
                                    aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                                     <  1:160521/53‑1 (MQ=255)
                                    aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                                     <  1:3092487/53‑1 (MQ=255)
                                       tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa                                >  1:1425652/1‑71 (MQ=255)
                                        ttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa                                >  1:4028751/1‑70 (MQ=255)
                                            gagTGATGAAGCCGAAGCCTTTGTCAGCGTTGAGCCATTTTACGATACCAGTCAtt                                          >  1:3514997/3‑56 (MQ=255)
                                            gagTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtt                                          >  1:4568935/3‑56 (MQ=255)
                                            gagTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtt                                          >  1:1857598/3‑56 (MQ=255)
                                              gTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                          >  1:2103233/1‑70 (MQ=255)
                                               tgatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCTCTAAACCAGTCATTTTGTTAGACATAatta                          <  1:4387484/69‑1 (MQ=255)
                                               tgataaAGCCAAAACCTTTATCTGCGTTAAACCAt                                                            <  1:5295674/35‑1 (MQ=12)
                                                 atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAttt                                         <  1:3203337/52‑1 (MQ=255)
                                                 atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                      <  1:5633836/55‑1 (MQ=255)
                                                 atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                      <  1:4986770/55‑1 (MQ=255)
                                                 atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                      <  1:4753505/55‑1 (MQ=255)
                                                 atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                      <  1:8492/55‑1 (MQ=255)
                                                 atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                      <  1:210327/55‑1 (MQ=255)
                                                  tgaAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtt                                          >  1:899380/1‑50 (MQ=255)
                                                   gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg                                       >  1:2865127/1‑52 (MQ=255)
                                                   gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg                                       >  1:4645415/1‑52 (MQ=255)
                                                   gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt                           <  1:1845520/64‑1 (MQ=255)
                                                    aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAATCCAGTCATTTTGTTAGACATAATTATTAc                      >  1:3841039/1‑68 (MQ=255)
                                                    aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtt                   >  1:923052/1‑71 (MQ=255)
                                                     aaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCttt                  >  1:4390881/1‑71 (MQ=255)
                                                       ccAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                <  1:2256864/71‑1 (MQ=255)
                                                       ccAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTc                                             >  1:1237644/1‑42 (MQ=25)
                                                       ccAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCa                                            >  1:3179028/1‑43 (MQ=25)
                                                       ccAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCa                                            >  1:1923790/1‑43 (MQ=25)
                                                        cAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtttt                                        >  1:5024372/1‑46 (MQ=255)
                                                        cAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                >  1:2594071/1‑70 (MQ=255)
                                                        cAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                >  1:5017903/1‑70 (MQ=255)
                                                         aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                      >  1:5375720/1‑47 (MQ=255)
                                                         aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                >  1:5498046/1‑69 (MQ=255)
                                                         aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                >  1:5448399/1‑69 (MQ=255)
                                                         aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                >  1:2074788/1‑69 (MQ=255)
                                                         aaaaCCTTTATCTGCGTTAAACCATTTTACTAAACCAGt                                              <  1:4281685/39‑1 (MQ=25)
                                                          aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa                                    >  1:2629414/1‑48 (MQ=255)
                                                          aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                >  1:3407416/1‑68 (MQ=255)
                                                          aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                >  1:3166722/1‑68 (MQ=255)
                                                           aaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaaga            <  1:405913/71‑1 (MQ=255)
                                                           aaCCTTTATCAGCGTTAAACCATTTTACTAAACCa                                                >  1:1136215/1‑35 (MQ=25)
                                                           aaCCTTTATCAGCGTTAAACCATTTTACTAAACCa                                                >  1:1012348/1‑35 (MQ=25)
                                                            aCCTTTTTCAGGGTTAAACCATTTCACTATACCAGTCATTTTGt                                      <  1:5152605/44‑1 (MQ=255)
                                                            aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat                         >  1:4752905/1‑57 (MQ=255)
                                                            aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat                         >  1:2418489/1‑57 (MQ=255)
                                                            aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat                         >  1:4194263/1‑57 (MQ=255)
                                                            aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaagaa           <  1:2837519/71‑1 (MQ=255)
                                                             ccTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGaaa          >  1:691397/1‑71 (MQ=255)
                                                               tttgtCAGCGTTGAACCATTTTACGATACCAGTCAttt                                         <  1:5566213/38‑1 (MQ=12)
                                                               tttgtCAGCGTTGAACCATTTTACGATACCAGTCAttt                                         <  1:2054396/38‑1 (MQ=12)
                                                               tttgtCAGCGTTGAACCATTTTACGATACCAGTCAttt                                         <  1:4874612/38‑1 (MQ=12)
                                                               tttatCAGCGTTAAACCATTTTACTAAACCAGTCAttt                                         <  1:2365632/38‑1 (MQ=25)
                                                               tttatCAGCGTTAAACCATTTTACTAAACCAGTCATTTTATTTGa                                  <  1:1796779/45‑1 (MQ=255)
                                                                 tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAtt        <  1:2627508/69‑1 (MQ=255)
                                                                 tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAtt        <  1:4040715/69‑1 (MQ=255)
                                                                 tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTa       <  1:5659966/70‑1 (MQ=255)
                                                                 tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAg      >  1:697425/1‑71 (MQ=255)
                                                                   tCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGc     <  1:2869713/70‑1 (MQ=255)
                                                                   tCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGc     <  1:1376651/70‑1 (MQ=255)
                                                                   tCAGCGTTAAACCATTTTACTAAACCAGTCATTTTAttt                                    <  1:3098439/39‑2 (MQ=25)
                                                                    cAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGc     <  1:778949/69‑1 (MQ=255)
                                                                     |                                                                      
AATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCT  >  minE/980194‑980333

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