Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,415,809 C→T 23.2% noncoding (70/77 nt) proM ← tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,415,8090CT23.2% 88.3 / 19.6 43noncoding (70/77 nt)proMtRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base C (18/15);  new base T (5/5);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

TGCGAAGGGAGGGTGACGAAATGCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCC  >  minE/2415739‑2415877
                                                                      |                                                                    
tGCGAAGGGAGGGTGACGAAATGCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGc                                                                      >  1:2154799/1‑71 (MQ=255)
 gCGAAGGGAGGGTGACGAAATGCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGc                                                                      >  1:3449086/1‑70 (MQ=255)
       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:3691230/71‑1 (MQ=255)
                        acaGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAcc                                              <  1:4330351/71‑1 (MQ=255)
                         caGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAccc                                             <  1:4732647/71‑1 (MQ=255)
                         caGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAccc                                             <  1:4726192/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:2281846/1‑59 (MQ=255)
                                                             aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                     >  1:1358257/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)
                                                                      |                                                                    
TGCGAAGGGAGGGTGACGAAATGCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCC  >  minE/2415739‑2415877

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: