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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,415,8120GT23.8% 99.2 / 16.8 42noncoding (67/77 nt)proMtRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base G (18/14);  new base T (5/5);  total (23/19)
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.42e-01

GGAGGGTGACGAAATGCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGA  >  minE/2415746‑2415879
                                                                  |                                                                   
ggagggTGACGAAATGCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGagagag                                                                 <  1:694645/71‑1 (MQ=255)
     gTGACGAAATGCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGa                                                               <  1:3141884/68‑1 (MQ=255)
               gCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTCGTCGGCGAGAGAGGATTCGAAc                                                        >  1:657077/1‑65 (MQ=255)
                cacaGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAc                                                 >  1:5593864/1‑71 (MQ=255)
                cacaGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAc                                                 >  1:1024987/1‑71 (MQ=255)
                 acaGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAcc                                                <  1:4330351/71‑1 (MQ=255)
                 acaGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAcc                                                <  1:3691230/71‑1 (MQ=255)
                  caGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAccc                                               <  1:4726192/71‑1 (MQ=255)
                  caGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAccc                                               <  1:4732647/71‑1 (MQ=255)
                  caGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAcc                                                >  1:3209176/1‑70 (MQ=255)
                    gAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCAc                                             >  1:2324492/1‑71 (MQ=255)
                        aCAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCa                                              >  1:4755699/1‑66 (MQ=255)
                        aCAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTgg                                          >  1:548331/1‑70 (MQ=255)
                            aaaaCCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTccc                                      >  1:3533273/1‑70 (MQ=255)
                            aaaaCCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCa                                     <  1:4770576/71‑1 (MQ=255)
                               aCCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAcc                                                <  1:5360046/57‑1 (MQ=255)
                                  ccgccgAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGa                                                          <  1:4475856/44‑1 (MQ=255)
                                       gAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                          >  1:3122587/1‑71 (MQ=255)
                                       gAAGCGGGGTTAAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                          >  1:3793516/1‑71 (MQ=255)
                                          gCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                          >  1:577056/1‑68 (MQ=255)
                                             gggTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAAcc                                 <  1:5097850/58‑1 (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)
                                                   aaaaTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCa                                              >  1:1731257/1‑39 (MQ=255)
                                                   aaaaTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCa                                <  1:1771869/53‑1 (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)
                                                       ttGGTCGGTGATAGAGGATTCGAACCTCCGAcccc                                              <  1:1540883/35‑2 (MQ=18)
                                                       ttGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAAcc                                 <  1:1723109/48‑1 (MQ=255)
                                                       ttGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAAcc                                 <  1:4283347/48‑1 (MQ=255)
                                                         ggTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTAc                      <  1:3823222/57‑1 (MQ=12)
                                                             ggCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGcc    <  1:1147009/71‑1 (MQ=255)
                                                              gCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGcc    >  1:1599202/1‑70 (MQ=255)
                                                                gagagaGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGa  >  1:3291001/1‑70 (MQ=255)
                                                                  |                                                                   
GGAGGGTGACGAAATGCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGA  >  minE/2415746‑2415879

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: