Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 158,183 A→G 75.8% N41S (AAC→AGC)  dgt → deoxyguanosine triphosphate triphosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE158,1830AG75.8% 120.1 / 36.8 62N41S (AAC→AGC) dgtdeoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (7/8);  new base G (27/20);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGC  >  minE/158118‑158251
                                                                 |                                                                    
gCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGctct                                                                  >  1:838170/1‑70 (MQ=255)
gCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGctct                                                                  >  1:409481/1‑70 (MQ=255)
gCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGAGATCGCGGGCGTATCATCAGctct                                                                  >  1:1774046/1‑70 (MQ=255)
       gTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAAt                                                          <  1:552569/71‑1 (MQ=255)
       gTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAAt                                                          <  1:495277/71‑1 (MQ=255)
       gTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAAt                                                          <  1:2558184/71‑1 (MQ=255)
              ccGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCaa                                                           >  1:1899319/1‑63 (MQ=255)
              ccGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCa                                                            >  1:2600426/1‑62 (MQ=255)
                gAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATtcgt                                                      <  1:1545908/66‑1 (MQ=255)
                gAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATtcgt                                                      <  1:1644034/66‑1 (MQ=255)
                 aaCATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTg                                                <  1:3566001/71‑1 (MQ=255)
                 aaCATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTg                                                <  1:3380443/71‑1 (MQ=255)
                     tGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaa                                             >  1:671470/1‑70 (MQ=255)
                       agaTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGcaacaa                                          <  1:932396/71‑1 (MQ=255)
                       agaTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGcaaca                                           <  1:3178224/70‑1 (MQ=255)
                       agaTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGcaaca                                           <  1:3427724/70‑1 (MQ=255)
                       agaTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaacaa                                          <  1:507292/71‑1 (MQ=255)
                         aTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa                                        <  1:1449693/71‑1 (MQ=255)
                         aTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa                                        <  1:2999894/71‑1 (MQ=255)
                         aTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa                                        <  1:2534584/71‑1 (MQ=255)
                            cTGCGGATCTTCGAGAGCGATCGCGTGCGTATCATCAGCTCTCCTGCAATTCGTCGTCTGCAACAAAAGAc                                     >  1:802025/1‑71 (MQ=255)
                            cTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTcgtcgt                                                   <  1:3336338/57‑1 (MQ=255)
                            cTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTcgtcgt                                                   <  1:1462751/57‑1 (MQ=255)
                            cTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTcgtcgt                                                   <  1:3927918/57‑1 (MQ=255)
                            cTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGa                                      >  1:3646674/1‑70 (MQ=255)
                            cTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGAc                                     >  1:4197889/1‑71 (MQ=255)
                            cTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGAc                                     >  1:288026/1‑71 (MQ=255)
                            cTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGAc                                     >  1:2132722/1‑71 (MQ=255)
                            cTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATAAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGAc                                     >  1:2742967/1‑71 (MQ=255)
                             tGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGAcc                                    >  1:427697/1‑71 (MQ=255)
                             tGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGAcc                                    >  1:2185059/1‑71 (MQ=255)
                                       cGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGtttt                            <  1:3671814/69‑1 (MQ=255)
                                       cGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGtttt                            <  1:1325572/69‑1 (MQ=255)
                                       cGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttt                             <  1:4020885/68‑1 (MQ=255)
                                        gagaGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGtttt                            <  1:4197196/68‑1 (MQ=255)
                                        gagaGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                          <  1:1200265/70‑1 (MQ=255)
                                          gagCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                       >  1:1590121/1‑71 (MQ=255)
                                          gagCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                       >  1:1350903/1‑71 (MQ=255)
                                           agCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGAc                                     >  1:1720485/1‑56 (MQ=255)
                                           agCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGAc                                     >  1:4178368/1‑56 (MQ=255)
                                           agCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGAc                                     >  1:2777435/1‑56 (MQ=255)
                                                 cgcgGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCa                        >  1:1207871/1‑63 (MQ=255)
                                                 cgcgGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCa                        >  1:1230690/1‑63 (MQ=255)
                                                 cgcgGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCa                        >  1:1622970/1‑63 (MQ=255)
                                                 cgcgGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCa                        >  1:2964535/1‑63 (MQ=255)
                                                  gcgGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACaaaa                                        >  1:1879181/1‑46 (MQ=255)
                                                  gcgGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACaaa                                         >  1:417801/1‑45 (MQ=255)
                                                  gcgGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc               <  1:3488188/71‑1 (MQ=255)
                                                  gcgGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcg                >  1:2549638/1‑70 (MQ=255)
                                                    gggCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCaa             <  1:3979954/71‑1 (MQ=255)
                                                    gggCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCa              <  1:3897997/70‑1 (MQ=255)
                                                     ggCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAAt            <  1:3743393/71‑1 (MQ=255)
                                                     ggCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAAt            <  1:1203037/71‑1 (MQ=255)
                                                         tATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                          <  1:3957245/53‑1 (MQ=255)
                                                          atcatcAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgc       >  1:398898/1‑71 (MQ=255)
                                                           tcatcaGCTCTCCGGCAATTCGTCGTCTGCAAAAAAAGACACAGGttt                             >  1:1045964/1‑48 (MQ=38)
                                                             atcaGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                       >  1:1462151/1‑52 (MQ=255)
                                                             atcaGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                       >  1:4210303/1‑52 (MQ=255)
                                                             atcaGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                       >  1:3181922/1‑52 (MQ=255)
                                                             atcaGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                       >  1:897566/1‑52 (MQ=255)
                                                             atcaGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                       >  1:3186892/1‑52 (MQ=255)
                                                             atcaGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTg                     >  1:1959697/1‑54 (MQ=255)
                                                                aaCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttttt                           >  1:2807009/1‑45 (MQ=255)
                                                                aaCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                          >  1:2350676/1‑46 (MQ=255)
                                                                aaCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                          >  1:4198639/1‑46 (MQ=255)
                                                                aaCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                          >  1:2843960/1‑46 (MQ=255)
                                                                aaCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                          >  1:3822887/1‑46 (MQ=255)
                                                                aaCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTgc  <  1:3080277/70‑1 (MQ=255)
                                                                 |                                                                    
GCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGC  >  minE/158118‑158251

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: