Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 980,248 A→G 16.5% G19G (GGT→GGC)  cspI ← cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE980,2480AG16.5% 236.0 / 34.7 103G19G (GGT→GGC) cspIcold shock protein
Reads supporting (aligned to +/- strand):  ref base A (42/44);  new base G (11/6);  total (53/50)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

Only 100 of 106 total aligned reads displayed.

AATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTT  >  minE/980181‑980317
                                                                   |                                                                     
aaTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                                                                    <  1:5691882/71‑1 (MQ=255)
 aTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                                                                    >  1:4960457/1‑70 (MQ=255)
 aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                                   <  1:4837128/71‑1 (MQ=255)
       gCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCt                                                              <  1:3456392/70‑1 (MQ=255)
          ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                          >  1:104868/1‑71 (MQ=255)
           tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                          >  1:1620991/1‑70 (MQ=255)
           tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                          >  1:128643/1‑70 (MQ=255)
           tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCATGATGAAACCaaa                                                                  <  1:5561120/62‑1 (MQ=255)
            gAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                                  >  1:4012266/1‑61 (MQ=255)
             aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                          <  1:3675894/68‑1 (MQ=255)
             aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                          <  1:1508089/68‑1 (MQ=255)
             aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                         <  1:2949205/69‑1 (MQ=255)
             aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                       <  1:4446119/71‑1 (MQ=255)
               ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                        <  1:2919994/68‑1 (MQ=255)
               ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                        <  1:4550533/68‑1 (MQ=255)
               ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                        <  1:4578964/68‑1 (MQ=255)
               ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                        <  1:5582488/68‑1 (MQ=255)
                 gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                                                                >  1:5680592/1‑58 (MQ=255)
                 gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                                    >  1:3491141/1‑70 (MQ=255)
                 gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                                    >  1:3022718/1‑70 (MQ=255)
                  cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                          <  1:5457685/63‑1 (MQ=255)
                  cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                          <  1:380875/63‑1 (MQ=255)
                  cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa                                                  <  1:3922659/71‑1 (MQ=255)
                  cTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                        >  1:2389092/1‑65 (MQ=255)
                   tGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                          <  1:3410184/62‑1 (MQ=255)
                    gagaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaaa                                                >  1:1216871/1‑71 (MQ=255)
                      gaAATGGACAAACTCATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                                                                >  1:2788336/1‑53 (MQ=37)
                       aaaTGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCa                                             <  1:2167906/71‑1 (MQ=255)
                       aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                               >  1:220450/1‑69 (MQ=255)
                         atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                           >  1:1490823/4‑71 (MQ=255)
                         atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                           >  1:1039078/4‑71 (MQ=255)
                          tggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                                          <  1:1379447/69‑1 (MQ=255)
                            gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                       >  1:2381801/1‑56 (MQ=255)
                            caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCAtttta                                        >  1:1453682/2‑69 (MQ=255)
                            caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACAAtttta                                        >  1:1298922/2‑69 (MQ=255)
                             aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                                  >  1:3443239/1‑44 (MQ=255)
                             aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                                  >  1:4190745/1‑44 (MQ=255)
                             aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                         <  1:4048922/53‑1 (MQ=255)
                             aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                         <  1:4785105/53‑1 (MQ=255)
                             aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                        <  1:395464/54‑1 (MQ=255)
                             aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                        <  1:4158835/54‑1 (MQ=255)
                               aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                                      >  1:4306035/1‑54 (MQ=255)
                               aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                     <  1:720516/71‑1 (MQ=255)
                               aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                     <  1:2315796/71‑1 (MQ=255)
                               aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                     <  1:5510559/71‑1 (MQ=255)
                               aaaCACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACt                                      >  1:5153911/1‑70 (MQ=255)
                               aaaCACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACt                                      >  1:5127617/1‑70 (MQ=255)
                                 acacATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTAcga                                     >  1:1052620/1‑67 (MQ=9)
                                 acacATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTAcga                                     >  1:4621677/1‑67 (MQ=9)
                                  cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                                          >  1:1567424/1‑63 (MQ=255)
                                  cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTc                                         >  1:1156252/1‑64 (MQ=255)
                                  cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                                  <  1:1513145/71‑1 (MQ=255)
                                    caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                       >  1:1659681/1‑48 (MQ=255)
                                    caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                       >  1:631642/1‑48 (MQ=255)
                                    caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                       >  1:5185732/1‑48 (MQ=255)
                                       cTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCa                                >  1:5511648/1‑68 (MQ=255)
                                        tttGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                               >  1:1263479/1‑36 (MQ=255)
                                        tttGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                               >  1:1587214/1‑36 (MQ=255)
                                        tttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                               <  1:5616151/52‑1 (MQ=255)
                                         ttGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                               >  1:4633623/1‑35 (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:82126/68‑1 (MQ=255)
                                            gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt                         <  1:2003423/68‑1 (MQ=255)
                                            gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt                         <  1:3892167/68‑1 (MQ=255)
                                              gcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                               <  1:3448645/46‑1 (MQ=255)
                                               ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                      <  1:4351269/70‑1 (MQ=255)
                                                 aTCTTTCGGTGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                     <  1:3328034/53‑1 (MQ=255)
                                                 aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                     <  1:3092487/53‑1 (MQ=255)
                                                 aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                     <  1:160521/53‑1 (MQ=255)
                                                    tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa                >  1:1425652/1‑71 (MQ=255)
                                                     ttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa                >  1:2488971/1‑70 (MQ=255)
                                                     ttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa                >  1:4028751/1‑70 (MQ=255)
                                                      tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt               >  1:2333624/1‑70 (MQ=255)
                                                       cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATaa             >  1:151545/1‑71 (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)
                                                             gatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatta          <  1:3978299/68‑1 (MQ=255)
                                                              atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAttt                         <  1:3203337/52‑1 (MQ=255)
                                                              atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                      <  1:4986770/55‑1 (MQ=255)
                                                              atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                      <  1:4753505/55‑1 (MQ=255)
                                                              atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                      <  1:5633836/55‑1 (MQ=255)
                                                              atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                      <  1:210327/55‑1 (MQ=255)
                                                              atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                      <  1:8492/55‑1 (MQ=255)
                                                               tgaAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtt                          >  1:899380/1‑50 (MQ=255)
                                                               tgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat            >  1:2128480/1‑64 (MQ=255)
                                                                gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg                       >  1:2865127/1‑52 (MQ=255)
                                                                gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg                       >  1:4645415/1‑52 (MQ=255)
                                                                gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt           <  1:1845520/64‑1 (MQ=255)
                                                                gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt           <  1:4098546/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)
                                                                   |                                                                     
AATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTT  >  minE/980181‑980317

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