Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,137,310 C→G 57.7% L187L (CTG→CTC truB ← tRNA pseudouridine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,137,3100CG57.7% 23.9 / 57.8 59L187L (CTG→CTCtruBtRNA pseudouridine synthase
Reads supporting (aligned to +/- strand):  ref base C (5/20);  new base G (18/16);  total (23/36)
Fisher's exact test for biased strand distribution p-value = 1.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCAT  >  minE/2137243‑2137374
                                                                   |                                                                
ccGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCAccc                                                                  <  1:1918714/68‑1 (MQ=255)
ccGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCAccc                                                                  <  1:2911770/68‑1 (MQ=255)
ccGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCAccc                                                                  <  1:2800029/68‑1 (MQ=255)
       tttACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGtcatc                                                          >  1:1122222/1‑69 (MQ=255)
       tttACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGtcat                                                           >  1:727160/1‑68 (MQ=255)
          acttacCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGtcatca                                                         <  1:1596268/67‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:1852467/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:1031119/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:940054/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:2561199/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:839282/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:3135169/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:3233151/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:746069/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:3381719/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:650198/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:709390/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:1215347/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:1195247/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:1573370/63‑1 (MQ=255)
                  ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:1660680/63‑1 (MQ=255)
                     gACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGa                                            <  1:968483/69‑1 (MQ=255)
                     gACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGa                                            <  1:2423797/69‑1 (MQ=255)
                        ggCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGatat                                         <  1:1082734/69‑1 (MQ=255)
                         gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGatat                                         <  1:1746818/68‑1 (MQ=255)
                         gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGatat                                         <  1:2261311/68‑1 (MQ=255)
                         gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGAtata                                        <  1:2173170/69‑1 (MQ=255)
                         gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGAtata                                        <  1:2297205/69‑1 (MQ=255)
                         gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGAtata                                        <  1:166751/69‑1 (MQ=255)
                          cgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGAtata                                        <  1:1760907/68‑1 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:1099491/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:3308389/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:355264/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:582229/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:3160136/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:658112/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:749617/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:1292105/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:3134899/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:1338693/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:1363783/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:2788442/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:1453502/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:194312/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:1886243/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                             >  1:1876516/1‑69 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTCCGGATATAAGTCCCtttt                              >  1:611783/1‑68 (MQ=255)
                                    tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCCTCAATGATAGTGCGGATAAAAGTGCCTTTTg                             >  1:3351491/1‑69 (MQ=255)
                                        aTGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                                     <  1:3383302/41‑1 (MQ=19)
                                           cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                      >  1:1065793/1‑69 (MQ=255)
                                           cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                      >  1:991624/1‑69 (MQ=255)
                                           cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                      <  1:793578/69‑1 (MQ=255)
                                           cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                      <  1:882267/69‑1 (MQ=255)
                                           cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                      <  1:315293/69‑1 (MQ=255)
                                           cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                      >  1:3088524/1‑69 (MQ=255)
                                           cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                      <  1:2206041/69‑1 (MQ=255)
                                           cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                      <  1:1095463/69‑1 (MQ=255)
                                                  caGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTtcc               <  1:2040328/69‑1 (MQ=255)
                                                  caGCCGAGTTTTTCACCCAGGGCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTtcc               <  1:2767693/69‑1 (MQ=255)
                                                     ccGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTtccagc            <  1:3264520/69‑1 (MQ=255)
                                                     ccGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTtccagc            <  1:890786/69‑1 (MQ=255)
                                                                aCCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCAt  <  1:188888/68‑1 (MQ=255)
                                                                   |                                                                
CCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCAT  >  minE/2137243‑2137374

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: