Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,837,474 G→C 11.7% I315M (ATC→ATG nlpD ← predicted outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,837,4740GC11.7% 185.6 / 13.0 68I315M (ATC→ATGnlpDpredicted outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (33/27);  new base C (2/6);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAA  >  minE/1837407‑1837544
                                                                   |                                                                      
gttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCATGTAATCATCATTATGTTtgatga                                                                       <  1:2130903/69‑1 (MQ=255)
gttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTtgatga                                                                       <  1:2880654/69‑1 (MQ=255)
gttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat                                                                      >  1:4082354/1‑70 (MQ=255)
gttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat                                                                      >  1:1436407/1‑70 (MQ=255)
 ttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATaa                                                                    <  1:1072264/71‑1 (MQ=255)
 ttgttCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATaa                                                                    <  1:77653/70‑1 (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:2241509/1‑71 (MQ=255)
            ccAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCg                                                         >  1:1733518/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:3613809/1‑71 (MQ=255)
                  ttGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcg                                                   >  1:690536/1‑71 (MQ=255)
                  tcgtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTAc                                                           <  1:4132440/61‑1 (MQ=255)
                    gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATgaga                                                               >  1:618454/1‑57 (MQ=255)
                    gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATgaga                                                               >  1:1107573/1‑57 (MQ=255)
                     tgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCacgc                                                <  1:447301/71‑1 (MQ=255)
                      gtCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcgcg                                                 >  1:1399836/1‑69 (MQ=255)
                       tCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAg                                              >  1:3758102/1‑71 (MQ=255)
                        cGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTAcc                                                          >  1:3222141/1‑58 (MQ=255)
                         gTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTAcc                                                          <  1:2710619/57‑1 (MQ=255)
                         attATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgct                                            <  1:3853165/70‑2 (MQ=255)
                         attATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgct                                            <  1:218434/70‑2 (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:3090907/1‑35 (MQ=255)
                                           cAGGTAATCATCATTATGTTTGATGATAATCAGAt                                                              >  1:2298775/1‑35 (MQ=255)
                                            aGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaa                          <  1:3012238/70‑1 (MQ=255)
                                            aGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaa                          <  1:2809451/70‑1 (MQ=255)
                                             ggTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAAcg                        <  1:1719454/71‑1 (MQ=255)
                                              gTAATCATCATTATGTTTGATGATAATCAGATTACCGt                                                        <  1:422719/38‑1 (MQ=255)
                                              gTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                            <  1:2403726/66‑1 (MQ=255)
                                              gTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                            <  1:1663837/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)
                                                         tATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCAt                 >  1:3851388/1‑66 (MQ=255)
                                                          aTGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaa                          >  1:2411442/1‑56 (MQ=255)
                                                          aTGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGc                    >  1:3075543/1‑62 (MQ=255)
                                                          aTGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGc                    >  1:1477141/1‑62 (MQ=255)
                                                           tGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGt          <  1:245045/71‑1 (MQ=255)
                                                             tttGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTcg        >  1:76041/1‑71 (MQ=255)
                                                             tttGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTc         >  1:2990616/1‑70 (MQ=255)
                                                               tgatgaTAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTcgc       <  1:1985776/70‑1 (MQ=255)
                                                               tgatgaTAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTcgc       >  1:2093652/1‑70 (MQ=255)
                                                               tgatgaTAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTcgc       >  1:4035365/1‑70 (MQ=255)
                                                               tgatgaTAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTcgc       >  1:979232/1‑70 (MQ=255)
                                                                gatgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGa     >  1:473133/1‑71 (MQ=255)
                                                                 atgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGAt    <  1:3954642/71‑1 (MQ=255)
                                                                 atgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGAt    <  1:3284165/71‑1 (MQ=255)
                                                                  tgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                            <  1:2175110/46‑1 (MQ=255)
                                                                  tgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATa   <  1:2834663/71‑1 (MQ=255)
                                                                   gatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGcc                   <  1:63106/54‑1 (MQ=255)
                                                                   gatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGcc                   <  1:1160973/54‑1 (MQ=255)
                                                                   gatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATaa  >  1:721381/1‑71 (MQ=255)
                                                                   |                                                                      
GTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAA  >  minE/1837407‑1837544

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: