Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,837,463 T→C 10.6% N319S (AAT→AGT)  nlpD ← predicted outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,837,4630TC10.6% 203.8 / 13.8 75N319S (AAT→AGT) nlpDpredicted outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (35/32);  new base C (2/6);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCAT  >  minE/1837394‑1837529
                                                                     |                                                                  
gCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATcatcat                                                                    >  1:2630090/1‑70 (MQ=255)
   ttAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTa                                                                  >  1:977391/1‑69 (MQ=255)
    tAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATcatcat                                                                    >  1:1026135/1‑66 (MQ=255)
    tAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGt                                                               >  1:1922686/1‑71 (MQ=255)
     aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                                              >  1:2198573/1‑71 (MQ=255)
     aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                                              >  1:2712619/1‑71 (MQ=255)
     aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                                              >  1:2833653/1‑71 (MQ=255)
     aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                                              >  1:1694774/1‑71 (MQ=255)
     aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                                              >  1:3200491/1‑71 (MQ=255)
            tgttgtTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTtg                                                            >  1:2977308/1‑66 (MQ=255)
             gttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCATGTAATCATCATTATGTTtgatga                                                        <  1:2130903/69‑1 (MQ=255)
             gttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTtgatga                                                        <  1:2880654/69‑1 (MQ=255)
             gttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat                                                       >  1:1436407/1‑70 (MQ=255)
             gttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat                                                       >  1:4082354/1‑70 (MQ=255)
              ttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATaa                                                     <  1:1072264/71‑1 (MQ=255)
              ttgttCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATaa                                                     <  1:77653/70‑1 (MQ=255)
                gttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATcatcat                                                                    >  1:2521679/1‑54 (MQ=255)
                gttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATa                                                      <  1:3941138/68‑1 (MQ=255)
                gttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAAt                                                    <  1:3212871/70‑1 (MQ=255)
                 ttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCa                                                  <  1:3954538/71‑1 (MQ=255)
                         ccAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCa                                                  <  1:1624757/63‑1 (MQ=255)
                         ccAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCg                                          >  1:1733518/1‑71 (MQ=255)
                         ccAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCg                                          >  1:2241509/1‑71 (MQ=255)
                          cAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCg                                          >  1:2681483/1‑70 (MQ=255)
                             cATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTAcc                                           >  1:2130067/1‑66 (MQ=255)
                             cATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGc                                      <  1:4123061/71‑1 (MQ=255)
                              atcgtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGc                                      <  1:102012/67‑1 (MQ=255)
                               ttGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTtgatg                                                         <  1:2810053/50‑1 (MQ=255)
                               ttGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAg                                       >  1:223673/1‑68 (MQ=255)
                               ttGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcg                                    >  1:690536/1‑71 (MQ=255)
                               ttGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcg                                    >  1:3613809/1‑71 (MQ=255)
                               tcgtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTAc                                            <  1:4132440/61‑1 (MQ=255)
                                 gtgtCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGttt                                                             >  1:2922585/1‑44 (MQ=255)
                                 gtgtCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTtgat                                                          >  1:1564132/1‑47 (MQ=255)
                                 gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATgaga                                                >  1:1107573/1‑57 (MQ=255)
                                 gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATgaga                                                >  1:618454/1‑57 (MQ=255)
                                  tgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCacgc                                 <  1:447301/71‑1 (MQ=255)
                                   gtCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                                              >  1:2530778/1‑41 (MQ=255)
                                   gtCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcgcg                                  >  1:1399836/1‑69 (MQ=255)
                                    tCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATg                                                                <  1:3982061/38‑1 (MQ=255)
                                    tCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGt                                                               <  1:256968/39‑1 (MQ=255)
                                    tCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAg                               >  1:3758102/1‑71 (MQ=255)
                                     cGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGt                                                               <  1:3245449/38‑1 (MQ=255)
                                     cGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGt                                                               <  1:1470513/38‑1 (MQ=255)
                                     cGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTAcc                                           >  1:3222141/1‑58 (MQ=255)
                                      gTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                                              >  1:874666/1‑38 (MQ=255)
                                      gTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                                              >  1:1278284/1‑38 (MQ=255)
                                      gTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTAcc                                           <  1:2710619/57‑1 (MQ=255)
                                      attATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgct                             <  1:218434/70‑2 (MQ=255)
                                      attATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgct                             <  1:3853165/70‑2 (MQ=255)
                                       ttATGGGCGTAGGCACTCAGGTAATCATCATTATGttt                                                             <  1:172997/38‑1 (MQ=255)
                                       ttATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgcgc                            <  1:375140/71‑1 (MQ=255)
                                         aTGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAAt                                                    >  1:1332449/1‑45 (MQ=255)
                                         aTGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGt                          >  1:1492891/1‑71 (MQ=255)
                                         aTGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGt                          >  1:3650222/1‑71 (MQ=255)
                                         aTGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGt                          >  1:3501462/1‑71 (MQ=255)
                                          tGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgcg                             <  1:2810475/67‑1 (MQ=255)
                                          tGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGtt                         <  1:971080/71‑1 (MQ=255)
                                          tGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGtt                         <  1:2535124/71‑1 (MQ=255)
                                            ggCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAAt                                                    <  1:973303/42‑1 (MQ=255)
                                                  ggCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCa                  >  1:2378536/1‑70 (MQ=255)
                                                   gCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATa                <  1:2449639/71‑1 (MQ=255)
                                                     aCTCAGGTAATCATCATTATGTTTGATGATAATCAg                                                 <  1:3845218/36‑1 (MQ=255)
                                                       tCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac             >  1:2795219/1‑70 (MQ=255)
                                                        cAGGTAATCATCATTATGTTTGATGATAATCAGAt                                               >  1:2298775/1‑35 (MQ=255)
                                                        cAGGTAATCATCATTATGTTTGATGATAATCAGAt                                               >  1:3090907/1‑35 (MQ=255)
                                                         aGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaa           <  1:2809451/70‑1 (MQ=255)
                                                         aGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaa           <  1:3012238/70‑1 (MQ=255)
                                                          ggTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAAcg         <  1:1719454/71‑1 (MQ=255)
                                                           gTAATCATCATTATGTTTGATGATAATCAGATTACCGt                                         <  1:422719/38‑1 (MQ=255)
                                                           gTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac             <  1:1663837/66‑1 (MQ=255)
                                                           gTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac             <  1:2403726/66‑1 (MQ=255)
                                                            tAATCATCATTATGTTTGATGATAATCAGATTACCGTAg                                       <  1:1779845/39‑1 (MQ=255)
                                                             aaTCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgc                              >  1:2243037/1‑47 (MQ=255)
                                                              atcatcATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCa                                >  1:2517615/1‑44 (MQ=255)
                                                               tcatcaTTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgcg                             <  1:2282723/46‑1 (MQ=255)
                                                               tcatcaTTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgcg                             <  1:607334/46‑1 (MQ=255)
                                                               tcatcaTTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgcg                             <  1:1166691/46‑1 (MQ=255)
                                                                catcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGt                          >  1:1875175/1‑48 (MQ=255)
                                                                 atcatTATGTTTGATGATAATCAGATTACCGTAGCCgcgc                                 <  1:3116356/40‑1 (MQ=255)
                                                                 atcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCa                                >  1:1988318/1‑41 (MQ=255)
                                                                 atcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCAt  >  1:1062104/1‑71 (MQ=255)
                                                                  tcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgcg                             >  1:2617998/1‑43 (MQ=255)
                                                                  tcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgc                              <  1:4062406/42‑1 (MQ=255)
                                                                    atTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTa                        >  1:3024430/1‑46 (MQ=255)
                                                                    atTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCa                     >  1:1320598/1‑49 (MQ=255)
                                                                     ttATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaac          >  1:835186/1‑59 (MQ=255)
                                                                     |                                                                  
GCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCAT  >  minE/1837394‑1837529

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: