Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,368,852 T→C 40.7% D188G (GAT→GGT)  dusC ← tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,368,8520TC40.7% 31.9 / 70.0 59D188G (GAT→GGT) dusCtRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base T (20/15);  new base C (8/16);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

GTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCAC  >  minE/1368784‑1368919
                                                                    |                                                                   
gttgCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATc                                                                    <  1:786755/70‑1 (MQ=255)
gttgCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCg                                                                   >  1:2192236/1‑71 (MQ=255)
   gcgcgcTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGTAATATTCAGCCGCTGGCGAATATCGcc                                                                 <  1:1957775/70‑1 (MQ=255)
   gcgcgcTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGcc                                                                 >  1:2062778/1‑70 (MQ=255)
             ccAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAatat                                                                     >  1:1804086/1‑56 (MQ=255)
             ccAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGc                                                       >  1:3764979/1‑70 (MQ=255)
              cAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGc                                                           >  1:1467795/1‑65 (MQ=255)
                   ccAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCAAGTCa                                                 >  1:2526650/1‑70 (MQ=255)
                     aGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAata                                              >  1:914188/1‑71 (MQ=255)
                      gATTTCACCGTTGGCAATCACCGGAATATTCAGCCGATGGCGAATATCGCCAATCGCCTGCCAGTCAatat                                             >  1:1247154/1‑71 (MQ=255)
                        tttCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGc                                           >  1:369568/1‑71 (MQ=255)
                        tttCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGc                                           <  1:1243800/71‑1 (MQ=255)
                         ttCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTc                                                  >  1:3963688/1‑63 (MQ=255)
                               gTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCcgcgc                                    <  1:2898966/71‑1 (MQ=255)
                                ttGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCg                                                                   >  1:3413215/1‑39 (MQ=255)
                                   gCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGt                                 <  1:1762085/70‑1 (MQ=255)
                                   gCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGt                                 <  1:537431/70‑1 (MQ=255)
                                       tCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAAccc                            >  1:99243/1‑71 (MQ=255)
                                        cACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCt                           >  1:2748662/1‑71 (MQ=255)
                                         aCCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCt                                          <  1:1138926/55‑1 (MQ=255)
                                             gAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTc                                                  <  1:3129602/43‑1 (MQ=255)
                                             gAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAata                                              >  1:1299838/1‑47 (MQ=255)
                                             gAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctct                      <  1:4114878/71‑1 (MQ=255)
                                             gAATATTCAGCCGCTGGCGAAAATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctct                      <  1:4125291/71‑1 (MQ=255)
                                             gAATATTCAGCCGCTGGCGAAAATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctct                      <  1:4129908/71‑1 (MQ=255)
                                               atatTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCaa                                                >  1:2601446/1‑43 (MQ=255)
                                                tatTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTc                   <  1:311052/71‑1 (MQ=255)
                                                  ttCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGt                 >  1:646995/1‑71 (MQ=255)
                                                  ttCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGt                 >  1:611395/1‑71 (MQ=255)
                                                   tCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAata                                              >  1:3760441/1‑41 (MQ=255)
                                                    cAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTcc               >  1:851985/1‑71 (MQ=255)
                                                     caccgcTGGCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgc              >  1:1725937‑M1/3‑19 (MQ=255)
                                                     caccgcTGGCGAATACCGCattcaggaaggaaatgcgagtgatttcacc                                    <  1:2766843‑M1/47‑31 (MQ=255)
                                                     aGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATg                                            >  1:3111066/1‑41 (MQ=255)
                                                     aGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATg                                            >  1:3130215/1‑41 (MQ=255)
                                                      gCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTcctc             <  1:2393672/71‑3 (MQ=255)
                                                      accgCTGGCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgc              >  1:3846587‑M1/2‑18 (MQ=255)
                                                      accgCTGGCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaa                       >  1:1609139‑M1/2‑18 (MQ=255)
                                                      accgCTGGCGAATACAGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaa                       >  1:2722785‑M1/2‑18 (MQ=255)
                                                       ccGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGcc            <  1:3561759/71‑1 (MQ=255)
                                                       ccGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGcc            <  1:2820872/71‑1 (MQ=255)
                                                        cGCTGGCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgcca            >  1:2174734‑M1/1‑16 (MQ=255)
                                                        cGCTGGCGAATACCGCattcaggaaggaaatgcgagtga                                           >  1:1131353‑M1/1‑16 (MQ=255)
                                                        cGCTGGCGAATACCGCattcaggaaggaaatgcgagt                                             <  1:3350481‑M1/37‑22 (MQ=255)
                                                         gCTGGCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaa          <  1:818966‑M1/71‑57 (MQ=255)
                                                         gCTGGCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaa          <  1:1901858‑M1/71‑57 (MQ=255)
                                                         gCTGGCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaa          <  1:1522625‑M1/71‑57 (MQ=255)
                                                         gCTGGCGAATACCGCattcaggaaggaaatgcgagtgatttc                                       <  1:1867202‑M1/42‑28 (MQ=255)
                                                          cTGGCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaa          >  1:3638509‑M1/1‑14 (MQ=255)
                                                          cTGGCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaa          >  1:754419‑M1/1‑14 (MQ=255)
                                                           tGGCGAATACCGCattcaggaaggaaatgcgagtga                                           <  1:1624554‑M1/36‑24 (MQ=255)
                                                           tGGCGAATACCGCattcaggaaggaaatgcgagtga                                           <  1:656680‑M1/36‑24 (MQ=255)
                                                             gCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaaga      <  1:3254188‑M1/71‑61 (MQ=255)
                                                             gCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaa        <  1:2312531‑M1/69‑59 (MQ=255)
                                                             gCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaa        <  1:3709041‑M1/69‑59 (MQ=255)
                                                             gCGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccgg                            <  1:659449‑M1/49‑39 (MQ=255)
                                                              cGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaaga      <  1:1995104‑M1/70‑61 (MQ=255)
                                                              cGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaa        <  1:696912‑M1/68‑59 (MQ=255)
                                                              cGAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaa        <  1:79003‑M1/68‑59 (MQ=255)
                                                               gAATACCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccga                   <  1:2697180‑M1/56‑48 (MQ=255)
                                                                 atatCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTg                          >  1:2998969/1‑47 (MQ=255)
                                                                 atatCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGc                         >  1:3329087/1‑48 (MQ=255)
                                                                 atacCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcg  <  1:937855‑M1/71‑65 (MQ=255)
                                                                 atacCGCattcaggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatc                 >  1:696642‑M1/1‑7 (MQ=255)
                                                                  tatCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCttt                    >  1:2872765/1‑52 (MQ=255)
                                                                  tatCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTc                   <  1:2781221/53‑1 (MQ=255)
                                                                  tatCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTc                   <  1:180762/53‑1 (MQ=255)
                                                                  tatCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGcaccac  >  1:1372096/1‑70 (MQ=255)
                                                                    |                                                                   
GTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCAC  >  minE/1368784‑1368919

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: