Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,952,505 (G)8→7 noncoding (72/87 nt) leuP ← tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,952,4980G.98.4% 239.4 / ‑1.6 61noncoding (79/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (50/10);  total (51/10)
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.57e-01

GCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAA  >  minE/2952430‑2952561
                                                                    |                                                               
gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                               <  1:2306/71‑1 (MQ=255)
gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                               <  1:2434880/71‑1 (MQ=255)
gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                               <  1:842780/71‑1 (MQ=255)
gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                               <  1:2205044/71‑1 (MQ=255)
gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                               <  1:2667585/71‑1 (MQ=255)
gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                               <  1:2176193/71‑1 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggac                                                         <  1:2309633/70‑3 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggac                                                         <  1:1944586/70‑3 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggac                                                         <  1:3819633/70‑3 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggac                                                         <  1:2906824/70‑3 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggac                                                         <  1:3406745/70‑3 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggac                                                         <  1:3420873/70‑3 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggac                                                         <  1:2205936/70‑3 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                           <  1:3855494/68‑1 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                           <  1:3661374/68‑1 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                           <  1:3324751/68‑1 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                           <  1:3360664/68‑1 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                           <  1:1919381/68‑1 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                           <  1:3579776/68‑1 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                           <  1:3563804/68‑1 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                           <  1:1361593/68‑1 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                           <  1:1481646/68‑1 (MQ=255)
       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                           <  1:189260/68‑1 (MQ=255)
        cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggac                                                         <  1:1699326/69‑3 (MQ=255)
              cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                           <  1:2242640/61‑1 (MQ=255)
                   ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACg                                           <  1:767533/71‑1 (MQ=255)
                   ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCGCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACg                                           <  1:46546/71‑1 (MQ=38)
                    cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACg                                           <  1:1736674/70‑1 (MQ=255)
                    cTTCCCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACg                                           <  1:1856297/70‑1 (MQ=38)
                              tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACg                                           <  1:2655267/60‑1 (MQ=255)
                                aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                           <  1:1365328/43‑1 (MQ=255)
                                aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                               <  1:259273/70‑1 (MQ=255)
                                aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                               <  1:4269048/70‑1 (MQ=255)
                                   ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAAccc                                                >  1:3605941/1‑50 (MQ=37)
                                   ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTaa                                    >  1:2700732/1‑62 (MQ=255)
                                   ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTaa                                    >  1:1950531/1‑62 (MQ=255)
                                   ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTaa                                    >  1:2007549/1‑62 (MQ=255)
                                   ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTaa                                    >  1:3189101/1‑62 (MQ=255)
                                   ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTaa                                    >  1:3706343/1‑62 (MQ=255)
                                   ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTaa                                    >  1:3393707/1‑62 (MQ=255)
                                   ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTaa                                    >  1:2196388/1‑62 (MQ=255)
                                   ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTaa                                    >  1:1040151/1‑62 (MQ=255)
                                                 atCGCCTCGTTTGGGGCGAGGGGGGGACTTGAACCCCCACg                                           <  1:2499725/41‑1 (MQ=18)
                                                    gCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                           >  1:3956265/1‑54 (MQ=255)
                                                    gCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                           >  1:219417/1‑54 (MQ=255)
                                                    gCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                           >  1:2187530/1‑54 (MQ=255)
                                                    gCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                           >  1:2154505/1‑54 (MQ=255)
                                                    gCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                           >  1:1719420/1‑54 (MQ=255)
                                                    gCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                           >  1:4117853/1‑54 (MQ=255)
                                                    gCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                           >  1:3375821/1‑54 (MQ=255)
                                                    gCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                           >  1:4194054/1‑54 (MQ=255)
                                                    gCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                           >  1:454321/1‑54 (MQ=255)
                                                    gCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                           >  1:2477755/1‑54 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTa                            >  1:1360788/1‑52 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGt                                          >  1:223972/1‑37 (MQ=14)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTa                                     >  1:2831362/1‑42 (MQ=14)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                               <  1:3533177/48‑1 (MQ=17)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                               <  1:3994454/48‑1 (MQ=17)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGt         >  1:2825123/1‑70 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:3829833/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:826982/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:365165/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:4005133/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:362543/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:976937/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:3551182/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:1233962/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:1464279/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:1813717/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:2029007/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:2059666/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:2091374/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:23725/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:2640742/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:2645802/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:2648615/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:2922022/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:3101579/1‑71 (MQ=255)
                                                      cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc        >  1:3214176/1‑71 (MQ=255)
                                                            tGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                         >  1:1840496/1‑48 (MQ=11)
                                                            tGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaa                   >  1:360019/1‑54 (MQ=255)
                                                            tGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa  >  1:4034284/1‑71 (MQ=255)
                                                            tGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa  >  1:4039486/1‑71 (MQ=255)
                                                            tGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa  >  1:3045665/1‑71 (MQ=255)
                                                            tGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa  >  1:3715943/1‑71 (MQ=255)
                                                            tGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa  >  1:1689516/1‑71 (MQ=255)
                                                            tGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACATGAAGCTAGCGCGTCTACCaa  >  1:2151412/1‑71 (MQ=11)
                                                                    |                                                               
GCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAA  >  minE/2952430‑2952561

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: