Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,264,627 T→C 12.5% K6K (AAA→AAG nusG ← transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,264,6270TC12.5% 185.2 / 16.2 72K6K (AAA→AAGnusGtranscription termination factor
Reads supporting (aligned to +/- strand):  ref base T (37/26);  new base C (5/4);  total (42/30)
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.65e-01

GCTCACGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGAATACCATC  >  minE/2264560‑2264696
                                                                   |                                                                     
gCTCACGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCtttt                                                                      >  1:1377569/1‑69 (MQ=255)
 cTCACGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCttt                                                                       <  1:2040349/67‑1 (MQ=255)
 cTCACGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTa                                                                    <  1:2464972/70‑1 (MQ=255)
        aGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCtt                                                             <  1:1986871/70‑1 (MQ=255)
        aGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTc                                                            <  1:740811/71‑1 (MQ=255)
          cGACGCTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTc                                                            >  1:1782609/1‑69 (MQ=255)
           gACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGa                                                         >  1:246593/1‑71 (MQ=255)
             cGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACa                                                       >  1:1089922/1‑71 (MQ=255)
                tGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTa                                                                    >  1:667016/1‑55 (MQ=255)
                tGCTACGAGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGa                                                                 <  1:446268/58‑1 (MQ=255)
                  cTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACAtct                                                    >  1:2522051/1‑69 (MQ=255)
                     cgcgGCCTTCAAAACCGGAAAACGCCTGAACGACTTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAg                                                 <  1:3827349/69‑1 (MQ=255)
                      gcgGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCtt                                                             >  1:3532733/1‑56 (MQ=255)
                      gcgGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCtt                                                             >  1:336062/1‑56 (MQ=255)
                        ggCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCt                                            >  1:1268572/1‑71 (MQ=255)
                        ggCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCt                                            <  1:1403653/71‑1 (MQ=255)
                         gCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTc                                           <  1:2056295/71‑1 (MQ=255)
                           cTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAg                                         >  1:3310809/1‑71 (MQ=255)
                           cTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAg                                         >  1:949684/1‑71 (MQ=255)
                            ttCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAg                                         >  1:3314492/1‑70 (MQ=255)
                              cAAAACCGGAAAACGCCTGAACGACGTACCAGCGcttctt                                                                     <  1:1230274/40‑1 (MQ=255)
                              cAAAACCGGAAAACGCCTGAACGACGTACCAGCGcttctt                                                                     <  1:1536873/40‑1 (MQ=255)
                              cAAAACCGGAAAACGCCTGAACGACGTACCAGCGCttt                                                                       >  1:1521743/1‑38 (MQ=255)
                              cAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGc                                       >  1:1222245/1‑70 (MQ=255)
                              cAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGc                                       >  1:2830570/1‑70 (MQ=255)
                              cAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGc                                       >  1:3365598/1‑70 (MQ=255)
                              cAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGc                                       >  1:3237980/1‑70 (MQ=255)
                                aaaCCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGa                                                                 >  1:3463974/1‑42 (MQ=255)
                                aaaCCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGa                                                                 >  1:1574735/1‑42 (MQ=255)
                                aaaCCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGa                                                                 >  1:3286380/1‑42 (MQ=255)
                                 caCCGGAAAACGCCTGAACGACGTACCAGCGCttt                                                                       <  1:2941321/34‑1 (MQ=255)
                                 aaCCGGAAAACGCCTGAACGACGTACCAGCGCttt                                                                       <  1:3872057/35‑1 (MQ=255)
                                 aaCCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGAc                                                        >  1:1901869/1‑50 (MQ=255)
                                  aCCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACAt                                                      >  1:2143928/1‑51 (MQ=255)
                                  aCCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGt                                   >  1:2378916/1‑70 (MQ=255)
                                  aCCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAg                                    >  1:612577/1‑69 (MQ=255)
                                  aCCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGt                                   >  1:281783/1‑70 (MQ=255)
                                  aCCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGt                                   >  1:2779133/1‑70 (MQ=255)
                                  aCCGGAAAACGCCTGAACGACGTACCAGCGATTTTTAGGAGCTTCAGACAt                                                      >  1:2662997/1‑51 (MQ=255)
                                    cGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGa                                 >  1:2407527/1‑70 (MQ=255)
                                     ggAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATa                               <  1:682528/71‑1 (MQ=255)
                                      gAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATaa                              <  1:630997/71‑1 (MQ=255)
                                       aaaaCGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGc                                       >  1:1266886/1‑61 (MQ=255)
                                       aaaaCGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATaa                              <  1:789594/70‑1 (MQ=255)
                                        aaaCGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACAtctc                                                   >  1:2207786/1‑48 (MQ=255)
                                           cGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAAc                                              >  1:1960002/1‑50 (MQ=255)
                                            gCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATaa                              <  1:144841/65‑1 (MQ=255)
                                             ccTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATaa                              <  1:468272/64‑1 (MQ=255)
                                             ccTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATAc                       <  1:3180449/71‑1 (MQ=255)
                                             ccTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATAc                       >  1:2343259/1‑71 (MQ=255)
                                              cTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATAcc                      <  1:2493191/71‑1 (MQ=255)
                                              cTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATAcc                      <  1:3416765/71‑1 (MQ=255)
                                              cTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCa                                          <  1:3856333/51‑1 (MQ=255)
                                              cTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCa                                          <  1:3612330/51‑1 (MQ=255)
                                               tGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACAt                                                      <  1:740231/38‑1 (MQ=255)
                                               tGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGa                          <  1:1649564/66‑1 (MQ=255)
                                               tGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATAc                       >  1:2755170/1‑69 (MQ=255)
                                                gAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTc                                           >  1:769047/1‑48 (MQ=255)
                                                gAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTc                                           >  1:2526206/1‑48 (MQ=255)
                                                gAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATa                               <  1:2466128/60‑1 (MQ=255)
                                                gAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATAc                       >  1:579203/1‑68 (MQ=255)
                                                gAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAg                    >  1:3836647/1‑71 (MQ=255)
                                                       gTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAAcc             <  1:1186225/71‑1 (MQ=255)
                                                        tACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAAc              >  1:18667/1‑69 (MQ=255)
                                                        tACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCa            <  1:403235/71‑1 (MQ=255)
                                                         aCCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAg                                    >  1:2087437/1‑46 (MQ=255)
                                                          ccAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGa                          >  1:1817577/1‑55 (MQ=255)
                                                          caaGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCa                                     >  1:2583910/3‑44 (MQ=255)
                                                           cAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGaa         >  1:1324486/1‑71 (MQ=255)
                                                           cAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGaa         <  1:1627940/71‑1 (MQ=255)
                                                           cAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGaa         <  1:173228/71‑1 (MQ=255)
                                                            aGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGAt                                <  1:3304165/47‑1 (MQ=255)
                                                            aGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGAt                                <  1:1053311/47‑1 (MQ=255)
                                                            aGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATaa                              <  1:149044/49‑1 (MQ=255)
                                                            aGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGaa         >  1:2136107/1‑70 (MQ=255)
                                                             gcgcTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGAATa       <  1:1328827/71‑1 (MQ=255)
                                                              cgcTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATaa                              <  1:3884031/47‑1 (MQ=255)
                                                                cTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGAt                                <  1:2942051/43‑1 (MQ=255)
                                                                cTTTTTAGGAGCTTCAGACAACTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGa                >  1:2788911/1‑59 (MQ=255)
                                                                   tttAGGCGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAgg                   >  1:509746/1‑53 (MQ=255)
                                                                   tttAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGAATACCATc  >  1:3283468/1‑70 (MQ=255)
                                                                   |                                                                     
GCTCACGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGAATACCATC  >  minE/2264560‑2264696

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: