Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 980,239 C→A 12.0% T22T (ACG→ACT)  cspI ← cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE980,2390CA12.0% 229.5 / 12.1 75T22T (ACG→ACT) cspIcold shock protein
Reads supporting (aligned to +/- strand):  ref base C (29/37);  new base A (2/7);  total (31/44)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATT  >  minE/980169‑980308
                                                                      |                                                                     
ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                       >  1:2636296/1‑71 (MQ=255)
ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                       <  1:1732959/71‑1 (MQ=255)
ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                       <  1:3820552/71‑1 (MQ=255)
   aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtg                                                                    <  1:1122383/71‑1 (MQ=255)
   aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtg                                                                    <  1:3160479/71‑1 (MQ=255)
        ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                <  1:2133592/70‑1 (MQ=255)
        ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                <  1:266283/70‑1 (MQ=255)
        ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                <  1:621809/70‑1 (MQ=255)
           aaaTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                            <  1:3570037/71‑1 (MQ=255)
           aaaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                            <  1:3829845/71‑1 (MQ=255)
            aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                       >  1:107367/1‑59 (MQ=255)
            aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                       >  1:3214165/1‑59 (MQ=255)
            aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                            <  1:79450/70‑1 (MQ=255)
             aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                          <  1:74153/71‑1 (MQ=255)
             aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                          <  1:1164842/71‑1 (MQ=255)
              tCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                          >  1:2807331/1‑70 (MQ=255)
                gTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                         <  1:1274742/69‑1 (MQ=255)
                   gCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtt                                                    >  1:3253011/1‑71 (MQ=255)
                      ctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGaa                                                               >  1:2839880/1‑57 (MQ=255)
                       tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                                                   >  1:2051136/1‑52 (MQ=255)
                       tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                            >  1:1881350/1‑59 (MQ=255)
                         aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                              >  1:2917685/1‑71 (MQ=255)
                         aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                              >  1:1779617/1‑71 (MQ=255)
                          aTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                             >  1:1237321/1‑71 (MQ=255)
                          aTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                             >  1:1634379/1‑71 (MQ=255)
                           ttGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                     <  1:1176248/46‑1 (MQ=255)
                              cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa                                         <  1:2521743/71‑1 (MQ=255)
                              cTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa                                         <  1:3571031/71‑1 (MQ=255)
                                gagaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                          <  1:2008069/52‑1 (MQ=255)
                                 agaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGaaa                                                              >  1:1054258/1‑47 (MQ=255)
                                  gaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtg                                                                    >  1:3158922/1‑40 (MQ=255)
                                  gaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaaa                                                        >  1:2162911/1‑52 (MQ=255)
                                  gaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                      >  1:3664975/1‑70 (MQ=255)
                                   aaaTGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCa                                    >  1:1007396/1‑71 (MQ=255)
                                    aaTGGACAAACACATCTTTGCTGCCATCTTTCGGAGTGATGAAACCAAAACCTTTTTCAg                                              <  1:494509/60‑1 (MQ=25)
                                     atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                  <  1:987032/68‑1 (MQ=255)
                                      tggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAggg                                            >  1:232045/3‑60 (MQ=255)
                                      tggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAggg                                            >  1:1273661/3‑60 (MQ=255)
                                      tAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                                 <  1:2211814/71‑1 (MQ=255)
                                       gcaCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCAttt                                 >  1:2548954/3‑70 (MQ=255)
                                        gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                          >  1:2141863/1‑44 (MQ=255)
                                        gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                          >  1:1701962/1‑44 (MQ=255)
                                        caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCAtttta                               <  1:658298/70‑3 (MQ=255)
                                          cAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                            <  1:2518770/40‑1 (MQ=255)
                                          cAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCAc                              >  1:1534400/1‑70 (MQ=255)
                                          cAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACt                             <  1:1450160/71‑1 (MQ=255)
                                          cAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACt                             <  1:2728011/71‑1 (MQ=255)
                                           aaaCACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACTa                            <  1:3525091/71‑1 (MQ=255)
                                            aaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTaa                           <  1:3226106/71‑1 (MQ=255)
                                            aaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTaa                           <  1:1430417/71‑1 (MQ=255)
                                             acacATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAggg                                            >  1:2235155/1‑53 (MQ=255)
                                             acacATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTaaa                          >  1:3760251/1‑71 (MQ=255)
                                               acgtCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAAcc                        <  1:2092435/68‑1 (MQ=9)
                                               acgtCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAAc                         >  1:3505348/4‑70 (MQ=9)
                                               acaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAcc                        <  1:3379269/71‑1 (MQ=255)
                                                       gctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAttt                >  1:1066236/1‑71 (MQ=255)
                                                       gctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtt                 <  1:2909648/70‑1 (MQ=255)
                                                        gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtttt               <  1:2536224/69‑1 (MQ=255)
                                                        gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtttt               <  1:3804793/69‑1 (MQ=255)
                                                         tgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                                 <  1:2182225/52‑1 (MQ=255)
                                                         tgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAt                                   <  1:531124/50‑1 (MQ=255)
                                                         agcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtttt               <  1:2864534/69‑1 (MQ=255)
                                                           ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCAc                              <  1:833754/53‑1 (MQ=255)
                                                           ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt             <  1:2761567/70‑1 (MQ=255)
                                                            cATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                               <  1:3158695/35‑1 (MQ=255)
                                                            cATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaa                                        >  1:2866869/1‑42 (MQ=255)
                                                            cATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCa                               >  1:400185/1‑51 (MQ=255)
                                                            cATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa           <  1:1984849/71‑1 (MQ=255)
                                                             aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                              <  1:2663073/35‑1 (MQ=255)
                                                             aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                         <  1:2536476/56‑1 (MQ=255)
                                                              tCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                             <  1:1487346/35‑1 (MQ=255)
                                                              tCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                             <  1:1740353/35‑1 (MQ=255)
                                                               cTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg              >  1:1881533/1‑65 (MQ=255)
                                                                tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtt                 >  1:3310296/1‑61 (MQ=255)
                                                                tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa       >  1:1521667/1‑71 (MQ=255)
                                                                   ctgCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATa     <  1:2200645/68‑1 (MQ=255)
                                                                   cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATaa    >  1:1203435/1‑71 (MQ=255)
                                                                   cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATaa    >  1:3118615/1‑71 (MQ=255)
                                                                     gCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtt                 <  1:3418694/56‑1 (MQ=255)
                                                                     gCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt  <  1:3671569/71‑1 (MQ=255)
                                                                      cGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat   <  1:3395416/69‑1 (MQ=255)
                                                                      |                                                                     
TTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATT  >  minE/980169‑980308

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: ‑