Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,581 T→G 6.1% K288Q (AAA→CAA) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,5810TG 6.1% 267.2 / 12.2 99K288Q (AAA→CAA) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (54/39);  new base G (4/2);  total (58/41)
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.54e-01

Only 100 of 105 total aligned reads displayed.

TCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAC  >  minE/1933514‑1933646
                                                                   |                                                                 
tCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                                                  >  1:959493/1‑69 (MQ=255)
      ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGc                                                          >  1:2941054/1‑71 (MQ=255)
      ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGAGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGc                                                          >  1:1184279/1‑71 (MQ=255)
        aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGc                                                          <  1:436360/69‑1 (MQ=255)
        aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCg                                                         >  1:1048979/1‑70 (MQ=255)
        aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCg                                                         >  1:3987714/1‑70 (MQ=255)
         ccGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                                >  1:3244948/1‑62 (MQ=255)
         ccGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGt                                                        >  1:3504792/1‑70 (MQ=255)
          cGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                                <  1:3081496/61‑1 (MQ=255)
           gagaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAg                                                     >  1:3251239/1‑71 (MQ=255)
           gagaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAg                                                     >  1:787454/1‑71 (MQ=255)
              aaTATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                                                  >  1:3586161/1‑71 (MQ=255)
                tatCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                                                >  1:493825/1‑71 (MQ=255)
                 atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                                    >  1:803082/1‑66 (MQ=255)
                 atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                                                >  1:1540856/1‑70 (MQ=255)
                 atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                                                >  1:710732/1‑70 (MQ=255)
                       tgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                                <  1:2126579/48‑1 (MQ=255)
                       tgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                         >  1:3752500/1‑71 (MQ=255)
                       tgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                         >  1:2024923/1‑71 (MQ=255)
                       tgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                         >  1:3873091/1‑71 (MQ=255)
                        gCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                                    >  1:311818/1‑59 (MQ=255)
                        gCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                                              <  1:743968/65‑1 (MQ=255)
                        gCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTc                                        >  1:808657/1‑71 (MQ=255)
                         ccAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTc                                        >  1:3769239/1‑70 (MQ=255)
                          cAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTaga                                            >  1:2596937/1‑65 (MQ=255)
                           agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                     <  1:2654301/71‑1 (MQ=255)
                           agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                     <  1:405187/71‑1 (MQ=255)
                            gagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTc                                        >  1:2552437/1‑67 (MQ=255)
                              gtgcGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCcatcat                                  <  1:4062946/68‑1 (MQ=255)
                              gTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                         <  1:4030477/64‑1 (MQ=255)
                               tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                                <  1:3224893/40‑1 (MQ=255)
                               tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                                    <  1:2990111/52‑1 (MQ=255)
                                gcGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCcatcattg                                >  1:805300/2‑69 (MQ=255)
                                aCGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg                                >  1:777398/1‑71 (MQ=255)
                                 cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATc                                                             >  1:1096470/1‑41 (MQ=255)
                                 cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCt                                                            >  1:533369/1‑42 (MQ=255)
                                 cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGcc                                                   >  1:2928884/1‑51 (MQ=255)
                                 cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                               >  1:582045/1‑71 (MQ=255)
                                 cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                               <  1:3369666/71‑1 (MQ=255)
                                   gAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcatc                                    >  1:4052373/1‑64 (MQ=255)
                                     aCAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCcatc                                    >  1:1070820/1‑62 (MQ=255)
                                     aCAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAc                           >  1:1482309/1‑71 (MQ=255)
                                      cAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcatc                                    <  1:272688/61‑1 (MQ=255)
                                      cAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAct                          >  1:3825783/1‑71 (MQ=255)
                                        gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                                           >  1:2408506/1‑52 (MQ=255)
                                        gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                        <  1:1261435/71‑1 (MQ=255)
                                          gATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTgg                      >  1:3818034/1‑71 (MQ=255)
                                           aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                                              <  1:1948008/46‑1 (MQ=255)
                                            ttCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGt                    >  1:2920777/1‑71 (MQ=255)
                                             tCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                              >  1:3181064/1‑60 (MQ=255)
                                             tCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGt                    >  1:831731/1‑70 (MQ=255)
                                              cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                         >  1:2669447/1‑48 (MQ=255)
                                              cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTgg                      >  1:909813/1‑67 (MQ=255)
                                                gTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATg                <  1:2209943/71‑1 (MQ=255)
                                                gTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATg                <  1:852328/71‑1 (MQ=255)
                                                   aGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTgg                      <  1:2262911/62‑1 (MQ=255)
                                                    gCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACg             <  1:3572921/70‑1 (MQ=255)
                                                    gCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGt            <  1:89657/71‑1 (MQ=255)
                                                    gCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGt            <  1:1667434/71‑1 (MQ=255)
                                                     cTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                                           >  1:400692/1‑39 (MQ=255)
                                                     cTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                         >  1:406979/1‑41 (MQ=255)
                                                     cTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTAt                 <  1:1580681/65‑1 (MQ=255)
                                                     cTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTa           >  1:18510/1‑71 (MQ=255)
                                                      tCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                              <  1:1292786/51‑1 (MQ=255)
                                                      tCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTAc          >  1:3990410/1‑71 (MQ=255)
                                                        gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTATGACGTACGGa       <  1:1870443/71‑1 (MQ=255)
                                                         aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACTTACGGa       <  1:533436/71‑1 (MQ=255)
                                                         aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa       <  1:968724/71‑1 (MQ=255)
                                                         aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa       <  1:1979392/71‑1 (MQ=255)
                                                         aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa       <  1:3216173/71‑1 (MQ=255)
                                                         aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa       <  1:3899779/71‑1 (MQ=255)
                                                         aTCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAAcag        >  1:3822874/1‑68 (MQ=255)
                                                         aTCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGa       <  1:2550477/71‑1 (MQ=255)
                                                          tCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTc                                        >  1:1651916/1‑37 (MQ=12)
                                                            aGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa     >  1:590360/1‑70 (MQ=255)
                                                            aGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa       <  1:3124829/68‑1 (MQ=255)
                                                            aGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAAt    <  1:2070404/71‑1 (MQ=255)
                                                            aGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAAt    <  1:310532/71‑1 (MQ=255)
                                                            aGATCTTTGTAAACTGCGTAAGCCATGTAGAGTTCCATCATGGt                               >  1:2915336/1‑44 (MQ=255)
                                                             gATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa       <  1:2337203/67‑1 (MQ=255)
                                                              aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcatcat                                  <  1:2705904/39‑1 (MQ=255)
                                                              aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                     <  1:1104750/36‑1 (MQ=255)
                                                              aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                     <  1:2762994/36‑1 (MQ=255)
                                                              aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                              <  1:2020246/43‑1 (MQ=255)
                                                              aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa       <  1:2238646/66‑1 (MQ=255)
                                                              aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAc  >  1:3417673/1‑71 (MQ=255)
                                                              aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAc  >  1:1151704/1‑71 (MQ=255)
                                                              aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAc  >  1:870379/1‑71 (MQ=255)
                                                              aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAc  >  1:3634999/1‑71 (MQ=255)
                                                              aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAc  >  1:3863377/1‑71 (MQ=255)
                                                              aTCTTTGTAATATGCGTAAGCCATGTAGAGTTCCATAATGGTGAACTCTGGGTTATg                <  1:1832478/57‑1 (MQ=255)
                                                              aTCATTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                     <  1:1360294/36‑1 (MQ=255)
                                                                cTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGAc              >  1:617691/1‑57 (MQ=255)
                                                                cTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa     >  1:356757/1‑66 (MQ=255)
                                                                  ttGTAATCTGCGTAAGCCATGTAGAGTTCcatcat                                  >  1:2554648/1‑35 (MQ=255)
                                                                  ttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATg                                 >  1:3621195/1‑36 (MQ=255)
                                                                  ttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATg                                 <  1:1406816/36‑1 (MQ=255)
                                                                  ttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATg                                 >  1:2032029/1‑36 (MQ=255)
                                                                   tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                        <  1:2881580/44‑1 (MQ=255)
                                                                   tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAc  <  1:1380633/66‑1 (MQ=255)
                                                                   |                                                                 
TCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAC  >  minE/1933514‑1933646

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: