Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,415,838 G→C 27.3% noncoding (41/77 nt) proM ← tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,415,8380GC27.3% 82.0 / 27.3 44noncoding (41/77 nt)proMtRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base G (19/13);  new base C (6/6);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

ACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTC  >  minE/2415770‑2415905
                                                                    |                                                                   
aCAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTgg                                                                    >  1:548331/1‑70 (MQ=255)
    aaaaCCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTccc                                                                >  1:3533273/1‑70 (MQ=255)
    aaaaCCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCa                                                               <  1:4770576/71‑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)
                           aaaaTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCa                                                          <  1:1771869/53‑1 (MQ=255)
                              tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc                                     <  1:4180925/70‑2 (MQ=16)
                              tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc                                     <  1:1308051/70‑2 (MQ=16)
                              tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcg                                       >  1:2772012/2‑69 (MQ=16)
                              tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcg                                       >  1:1922640/2‑69 (MQ=16)
                              aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                 >  1:1358257/1‑59 (MQ=255)
                              aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                 >  1:2281846/1‑59 (MQ=255)
                              aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                 >  1:4671372/1‑59 (MQ=255)
                               ttGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAAcc                                                           <  1:4283347/48‑1 (MQ=255)
                               ttGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAAcc                                                           <  1:1723109/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)
                                        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:4269871/71‑1 (MQ=255)
                                                           cGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt        <  1:2270917/71‑1 (MQ=255)
                                                           cGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt        <  1:1273409/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)
                                                                    |                                                                   
ACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTC  >  minE/2415770‑2415905

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: