Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,140,377 T→A 11.1% D244V (GAT→GTT)  infB ← fused protein chain initiation factor 2, IF2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,140,3770TA11.1% 238.7 / 20.0 90D244V (GAT→GTT) infBfused protein chain initiation factor 2, IF2
Reads supporting (aligned to +/- strand):  ref base T (31/49);  new base A (5/5);  total (36/54)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TTTGTTGCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAGTCGG  >  minE/2140307‑2140447
                                                                      |                                                                      
tttgttGCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAt                                                                        >  1:230167/1‑71 (MQ=255)
  tgttgCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATc                                                                       >  1:1704543/1‑70 (MQ=255)
    ttgCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCt                                                                    <  1:1372049/71‑1 (MQ=255)
     tgCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGc                                                                     >  1:460325/1‑69 (MQ=255)
       ccTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTc                                                                 <  1:2070082/71‑1 (MQ=255)
        cTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCttt                                                                  >  1:3677323/1‑69 (MQ=255)
        cTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCg                                                                >  1:3842567/1‑71 (MQ=255)
          ttcttcGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGtc                                                              >  1:160476/1‑71 (MQ=255)
           tcttcGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAt                                                                        >  1:3146660/1‑60 (MQ=255)
           tcttcGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGtct                                                             >  1:1090174/1‑71 (MQ=255)
           tcttcGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGtc                                                              <  1:276013/70‑1 (MQ=255)
             ttcGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGtcttc                                                           >  1:493716/1‑71 (MQ=255)
               cGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGtctt                                                            <  1:788259/68‑1 (MQ=255)
               cGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGtcttct                                                          <  1:2489949/70‑1 (MQ=255)
                 gACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGcc                                                       <  1:1205661/71‑1 (MQ=255)
                   cgcgCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTg                                                     <  1:415165/71‑1 (MQ=255)
                   cgcgCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTg                                                     <  1:2686128/71‑1 (MQ=255)
                    gcgcTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCTCTTTCGTCTTCTGCCTg                                                     <  1:1214311/70‑1 (MQ=255)
                    gcgcTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTg                                                     <  1:2447316/70‑1 (MQ=255)
                    gcgcTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTgg                                                    >  1:3897163/1‑71 (MQ=255)
                     cgcTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTc                                                                 <  1:322260/57‑1 (MQ=255)
                         gctTTCGCGTCACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCtt                                                                   <  1:2624612/51‑1 (MQ=255)
                          ctTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAg                                                <  1:2837696/69‑1 (MQ=255)
                          ctTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGc                                               <  1:2954257/70‑1 (MQ=255)
                          ctTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCa                                              >  1:1538349/1‑71 (MQ=255)
                          ctTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCa                                              >  1:3188144/1‑71 (MQ=255)
                            ttCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCAt                                             >  1:3554874/1‑70 (MQ=255)
                             tCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGc                                                   >  1:2017923/1‑63 (MQ=255)
                             tCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATg                                            <  1:422260/70‑1 (MQ=255)
                                 gTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCt                                                                    <  1:1819560/42‑1 (MQ=255)
                                 gTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCt                                                                    <  1:2761721/42‑1 (MQ=255)
                                 gTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAg                                                <  1:3187810/62‑1 (MQ=255)
                                 gTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                       <  1:1464170/71‑1 (MQ=255)
                                 gTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                       <  1:718016/71‑1 (MQ=255)
                                 gTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCAAGt                                           >  1:1928051/1‑67 (MQ=255)
                                  ttACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCa                                              >  1:997239/1‑63 (MQ=255)
                                  ttACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                       <  1:2367835/70‑1 (MQ=255)
                                   tACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGt                                                               >  1:915950/1‑45 (MQ=255)
                                   tACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGTCTGGCGAGCATGTTgag                                       >  1:652772/1‑69 (MQ=255)
                                   tACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                       >  1:1768729/1‑69 (MQ=255)
                                    aCGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGtt                                          >  1:485012/1‑65 (MQ=255)
                                    aCGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAag                                    >  1:2330502/1‑71 (MQ=255)
                                    aCGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAag                                    >  1:2703072/1‑71 (MQ=255)
                                     cGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAagt                                   <  1:2681875/71‑1 (MQ=255)
                                     cGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAag                                    <  1:3485994/70‑1 (MQ=255)
                                       aCCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTg                                                         <  1:451639/47‑1 (MQ=255)
                                         cacggccacggcCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGcc                                                       >  1:3968191/1‑47 (MQ=255)
                                         cacggccacggcCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTg                               <  1:1904078/71‑1 (MQ=255)
                                         cacggccacggcCGCCTTCGACTTCACGAACGCTTTCGTCTTCCGCCTGGCGAGCATGTTGAGAagtagt                                >  1:1472123/1‑70 (MQ=255)
                                            ggccacggccGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACg                            <  1:3653677/71‑1 (MQ=255)
                                             gccacggccGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGt                           <  1:2552257/71‑1 (MQ=255)
                                             gccacggccGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGt                           <  1:4139089/71‑1 (MQ=255)
                                             gccacggccGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGt                           <  1:1141604/71‑1 (MQ=255)
                                             gccacggccGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGt                           <  1:2132710/71‑1 (MQ=255)
                                              ccacggccGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGt                           <  1:916714/70‑1 (MQ=255)
                                                acggccGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGa                         >  1:2782174/1‑70 (MQ=255)
                                                 cggccGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATa                       <  1:2262967/71‑1 (MQ=255)
                                                 cggccGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATa                       <  1:1620568/71‑1 (MQ=255)
                                                    ccgccTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAagta                                  >  1:1418210/1‑57 (MQ=255)
                                                     cgccTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                       <  1:2220014/51‑1 (MQ=255)
                                                      gccTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCa                                              <  1:2198922/43‑1 (MQ=255)
                                                         ttCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                       >  1:2845967/1‑47 (MQ=255)
                                                         ttCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTg                          <  1:3160695/60‑1 (MQ=255)
                                                         ttCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTg                          <  1:2424450/60‑1 (MQ=255)
                                                         ttCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTg                          <  1:199974/60‑1 (MQ=255)
                                                         ttCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTg                          <  1:2774448/60‑1 (MQ=255)
                                                         ttCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCg                   <  1:4018164/67‑1 (MQ=255)
                                                         ttCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAag                                    <  1:2271277/50‑1 (MQ=255)
                                                         ttCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAag                                    <  1:2276058/50‑1 (MQ=255)
                                                          tCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTgg               >  1:3264173/1‑70 (MQ=255)
                                                          tCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTgg               >  1:3888524/1‑70 (MQ=255)
                                                           cGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGa                              >  1:928205/1‑54 (MQ=255)
                                                           cGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGaa             >  1:3749343/1‑71 (MQ=255)
                                                           agACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGa              >  1:2959490/2‑70 (MQ=255)
                                                             aCTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTgg               <  1:851556/67‑1 (MQ=255)
                                                              cTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCt          <  1:1757000/71‑1 (MQ=255)
                                                                 cACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                       <  1:1343898/39‑1 (MQ=255)
                                                                 cACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAAt            >  1:424384/1‑66 (MQ=255)
                                                                 cACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCtt         <  1:4091066/69‑1 (MQ=255)
                                                                 cACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCtt         <  1:846296/69‑1 (MQ=255)
                                                                 cACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCtt         >  1:1601054/1‑69 (MQ=255)
                                                                 cACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCtt         >  1:35549/1‑69 (MQ=255)
                                                                  aCGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                       <  1:1718053/38‑1 (MQ=255)
                                                                  aCGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAAt                     >  1:3543186/1‑56 (MQ=255)
                                                                  aCGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAAt                     >  1:2884161/1‑56 (MQ=255)
                                                                  aCGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCt          <  1:854779/67‑1 (MQ=255)
                                                                  aCGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAg      >  1:3659386/1‑71 (MQ=255)
                                                                     aTCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                       <  1:3662326/35‑1 (MQ=255)
                                                                     aTCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                       <  1:769128/35‑1 (MQ=255)
                                                                     aTCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                       <  1:3663259/35‑1 (MQ=255)
                                                                     aTCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                       <  1:2568588/35‑1 (MQ=255)
                                                                     aTCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTg                               <  1:126231/43‑1 (MQ=255)
                                                                      tCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAGTc    >  1:1095714/1‑69 (MQ=255)
                                                                      tCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAGTCgg  >  1:2442217/1‑71 (MQ=255)
                                                                      |                                                                      
TTTGTTGCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAGTCGG  >  minE/2140307‑2140447

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: