Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,415,851 T→G 20.7% noncoding (28/77 nt) proM ← tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,415,8510TG20.7% 115.8 / 30.5 58noncoding (28/77 nt)proMtRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base T (30/16);  new base G (6/6);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCG  >  minE/2415785‑2415920
                                                                  |                                                                     
gAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                                                                   >  1:3122587/1‑71 (MQ=255)
gAAGCGGGGTTAAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                                                                   >  1:3793516/1‑71 (MQ=255)
   gCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                                                                   >  1:577056/1‑68 (MQ=255)
       ttttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCag                                                            >  1:3108423/3‑70 (MQ=12)
         ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCag                                                            <  1:2367851/69‑2 (MQ=12)
         ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGc                                                          >  1:5014009/1‑71 (MQ=12)
         ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGc                                                          >  1:468949/1‑71 (MQ=12)
         ttCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgc                                                                    <  1:4029939/61‑1 (MQ=255)
           cAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                                >  1:3907592/1‑63 (MQ=255)
           cAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTg                                                        >  1:4132185/1‑71 (MQ=255)
               tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc                                                    <  1:1308051/70‑2 (MQ=16)
               tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc                                                    <  1:4180925/70‑2 (MQ=16)
               tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcg                                                      >  1:2772012/2‑69 (MQ=16)
               tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcg                                                      >  1:1922640/2‑69 (MQ=16)
               aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                                >  1:4671372/1‑59 (MQ=255)
               aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                                >  1:1358257/1‑59 (MQ=255)
               aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                                >  1:2281846/1‑59 (MQ=255)
                  ggTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTAc                                                               <  1:3823222/57‑1 (MQ=12)
                      ggCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGcc                                             <  1:1147009/71‑1 (MQ=255)
                       gCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGcc                                             >  1:1599202/1‑70 (MQ=255)
                         gatagagGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGaa                                          >  1:1185964/4‑71 (MQ=12)
                         gagagaGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGa                                           >  1:3291001/1‑70 (MQ=255)
                          atagagGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTAc                                                               <  1:3747175/47‑1 (MQ=9)
                             gagGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc                                                    <  1:4244743/57‑2 (MQ=12)
                                    cGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCtttt                               >  1:2900712/1‑71 (MQ=255)
                                           ccGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCt                                                    >  1:561022/1‑43 (MQ=255)
                                           ccGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTa                                                   >  1:2476629/1‑44 (MQ=255)
                                           ccGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCg                                            <  1:3706236/51‑1 (MQ=255)
                                           ccGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAttt                        >  1:2233676/1‑71 (MQ=255)
                                            cGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                       <  1:1273409/71‑1 (MQ=255)
                                            cGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                       <  1:4269871/71‑1 (MQ=255)
                                            cGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                       <  1:2270917/71‑1 (MQ=255)
                                             gACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGaa                                          >  1:5520604/1‑51 (MQ=255)
                                             gACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                       >  1:291555/1‑70 (MQ=255)
                                                  aCTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATAc                                      >  1:3173002/1‑50 (MQ=255)
                                                   cTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGaaa                                         <  1:4971219/46‑1 (MQ=255)
                                                   cTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGaaa                                         <  1:2936474/46‑1 (MQ=255)
                                                   cTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTc                 >  1:1032025/1‑70 (MQ=255)
                                                     ggTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATa                                       >  1:1049354/1‑46 (MQ=255)
                                                     ggTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAt                          >  1:2381708/1‑59 (MQ=255)
                                                      gTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCaa               <  1:3630416/69‑1 (MQ=255)
                                                         ccAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCtt          >  1:3243/1‑71 (MQ=255)
                                                         ccAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCt           <  1:2149250/70‑1 (MQ=255)
                                                         ccAAACCAGTTGCGCTACCAAGCTGAGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCt           <  1:2672216/70‑1 (MQ=255)
                                                           aaaCCAGTTGCGCTACGAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCttt         >  1:2451534/1‑70 (MQ=255)
                                                           aaaCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCttt         >  1:2421618/1‑70 (MQ=255)
                                                           aaaCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCttt         >  1:5201225/1‑70 (MQ=255)
                                                           aaaCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCtt          >  1:2426201/1‑69 (MQ=255)
                                                             aCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGc                                   >  1:3526204/1‑42 (MQ=255)
                                                             aCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTa                     >  1:270388/1‑56 (MQ=255)
                                                             aCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTaaa      <  1:1735515/71‑1 (MQ=255)
                                                                aGTTTCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAttt                        >  1:5364138/1‑50 (MQ=255)
                                                                aGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAttt                        >  1:255122/1‑50 (MQ=255)
                                                                aGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAttt                        >  1:5595757/1‑50 (MQ=255)
                                                                 gTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCg  <  1:1975606/71‑1 (MQ=255)
                                                                 gTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCg  <  1:289974/71‑1 (MQ=255)
                                                                 gTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCg  <  1:5412913/71‑1 (MQ=255)
                                                                 gTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCg  <  1:1000378/71‑1 (MQ=255)
                                                                  |                                                                     
GAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCG  >  minE/2415785‑2415920

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