Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,883,630 G→C 41.1% Q285E (CAA→GAA)  rumA ← 23S rRNA (uracil‑5)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,883,6300GC41.1% 36.8 / 72.6 73Q285E (CAA→GAA) rumA23S rRNA (uracil‑5)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (20/23);  new base C (17/13);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

GTCGCCAATGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGA  >  minE/1883566‑1883691
                                                                |                                                             
gtCGCCAATGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGca                                                             <  1:1623564/67‑1 (MQ=255)
      aaTGGCAGTGTAAAGTTGCCCATACCTCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:10626/67‑1 (MQ=255)
      aaTGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:1009206/67‑1 (MQ=255)
      aaTGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:3029565/67‑1 (MQ=255)
      aaTGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:371509/67‑1 (MQ=255)
      aaTGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:2984740/67‑1 (MQ=255)
      aaTGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:451456/67‑1 (MQ=255)
      aaTGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:535502/67‑1 (MQ=255)
      aaTGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:1359582/67‑1 (MQ=255)
      aaTGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:2548265/67‑1 (MQ=255)
      aaTGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:2383769/67‑1 (MQ=255)
      aaTGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:151202/67‑1 (MQ=255)
      aaTGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:2309066/67‑1 (MQ=255)
         ggCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCAt                                                  >  1:680898/1‑69 (MQ=255)
         ggCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCAt                                                  >  1:3258151/1‑69 (MQ=255)
         ggCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCAt                                                  >  1:2114615/1‑69 (MQ=255)
         ggCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCAt                                                  >  1:2680706/1‑69 (MQ=255)
         ggCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCAt                                                  >  1:1280748/1‑69 (MQ=255)
         ggCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCAt                                                  >  1:902324/1‑69 (MQ=255)
            aGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTcc                                               <  1:1835266/69‑1 (MQ=255)
            aGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTcc                                               <  1:1965688/69‑1 (MQ=255)
            aGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTcc                                               <  1:1108737/69‑1 (MQ=255)
                          cATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCAcgcg                                     >  1:1715452/1‑65 (MQ=255)
                          cATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTAc                                 >  1:796875/1‑69 (MQ=255)
                          cATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTAc                                 >  1:1254836/1‑69 (MQ=255)
                          cATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTAc                                 >  1:1218070/1‑69 (MQ=255)
                          cATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTAc                                 >  1:2970540/1‑69 (MQ=255)
                          cATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTAc                                 >  1:3433805/1‑69 (MQ=255)
                          cATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTAc                                 >  1:2171088/1‑69 (MQ=255)
                          cATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTAc                                 >  1:330057/1‑69 (MQ=255)
                          cATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTAc                                 >  1:3187586/1‑69 (MQ=255)
                          cATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTAc                                 >  1:1882575/1‑69 (MQ=255)
                          cATACCGCAGAACAGATCCAGAACGCGATCTTCAGCTTGCACAt                                                          >  1:1171626/1‑44 (MQ=38)
                               cGCAGAACAGATCCAGTACGCGATCTTCAGGTTCCACATCCa                                                       <  1:2035521/42‑1 (MQ=255)
                                      cAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt                     <  1:2226747/69‑1 (MQ=255)
                                      cAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt                     <  1:3208038/69‑1 (MQ=255)
                                      cAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt                     <  1:328782/69‑1 (MQ=255)
                                      cAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt                     <  1:2602291/69‑1 (MQ=255)
                                      cAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt                     <  1:1827694/69‑1 (MQ=255)
                                      cAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt                     <  1:268530/69‑1 (MQ=255)
                                       aGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt                     <  1:1804816/68‑1 (MQ=255)
                                       aGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt                     <  1:2695644/68‑1 (MQ=255)
                                       aGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt                     <  1:2826666/68‑1 (MQ=255)
                                          tCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTcaca                 >  1:2602469/1‑69 (MQ=255)
                                          tCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTcaca                 >  1:1624337/1‑69 (MQ=255)
                                              gTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCg             <  1:2121798/69‑1 (MQ=255)
                                              gTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCg             <  1:2042004/69‑1 (MQ=255)
                                                 cgcgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt                     <  1:300710/58‑1 (MQ=255)
                                                   cgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCTCTACCATTTTTTGGTTCACACCCGCATTg        <  1:2519852/69‑1 (MQ=255)
                                                   cgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCTCTACCATTTTTTGGTTCACACCCGCATTg        <  1:2515780/69‑1 (MQ=255)
                                                   cgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTg        <  1:1506049/69‑1 (MQ=255)
                                                   cgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTg        <  1:2263245/69‑1 (MQ=255)
                                                   cgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTg        <  1:3124318/69‑1 (MQ=255)
                                                   cgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTg        <  1:2072153/69‑1 (MQ=255)
                                                      tCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt                     <  1:1574612/53‑1 (MQ=255)
                                                       cTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGAcc     >  1:2771824/1‑68 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTAcc                                >  1:1133627/1‑38 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTg   >  1:662308/1‑67 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTg   >  1:1031326/1‑67 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTg   >  1:3152075/1‑67 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTg   >  1:2880100/1‑67 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTg   >  1:2758002/1‑67 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa  >  1:3220960/1‑68 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa  >  1:3212404/1‑68 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa  >  1:3079580/1‑68 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa  >  1:3434285/1‑68 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa  >  1:3008571/1‑68 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa  >  1:2742140/1‑68 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa  >  1:1720517/1‑68 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa  >  1:697053/1‑68 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa  >  1:1640269/1‑68 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa  >  1:1424200/1‑68 (MQ=255)
                                                          cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa  >  1:908470/1‑68 (MQ=255)
                                                                gCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt                     <  1:3378823/43‑1 (MQ=255)
                                                                |                                                             
GTCGCCAATGGCAGTGTAAAGTTGCCCATACCGCAGAACAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGA  >  minE/1883566‑1883691

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