Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,837,474 G→C 36.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,4740GC36.7% 66.1 / 75.2 79I315M (ATC→ATGnlpDpredicted outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (28/22);  new base C (17/12);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

TGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAA  >  minE/1837406‑1837544
                                                                    |                                                                      
tgttgtTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTtgatga                                                                       >  1:3877715/1‑70 (MQ=255)
 gttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTtgatga                                                                       <  1:5424680/69‑1 (MQ=255)
 gttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTtgatga                                                                       <  1:4754187/69‑1 (MQ=255)
 gttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat                                                                      <  1:4205122/70‑1 (MQ=255)
  ttgttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATa                                                                     <  1:343340/70‑1 (MQ=255)
     ttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATc                                                                  <  1:1129091/70‑1 (MQ=255)
     ttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATc                                                                  <  1:470958/70‑1 (MQ=255)
     ttCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCa                                                                 <  1:4905618/71‑1 (MQ=255)
         tttaCCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGAt                                                              <  1:940574/67‑1 (MQ=255)
         tttaCCAGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGAt                                                              <  1:216349/67‑1 (MQ=255)
         cGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGAtt                                                             <  1:2104168/71‑1 (MQ=255)
              cAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGt                                                        <  1:210467/71‑1 (MQ=255)
               aGCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGt                                                        <  1:2273518/70‑1 (MQ=255)
                gCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat                                                                      <  1:3767932/55‑1 (MQ=255)
                gCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat                                                                      <  1:5178212/55‑1 (MQ=255)
                gCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat                                                                      <  1:760127/55‑1 (MQ=255)
                gCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAg                                                      <  1:4288213/71‑1 (MQ=255)
                gCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAg                                                      <  1:969777/71‑1 (MQ=255)
                gCATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAg                                                      >  1:4278880/1‑71 (MQ=255)
                 cATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGc                                                     <  1:169398/71‑1 (MQ=255)
                 cATCGTGTCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGc                                                     >  1:5110646/1‑71 (MQ=255)
                  atcgtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGc                                                     >  1:2584178/4‑70 (MQ=255)
                   ttGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCg                                                         >  1:2007803/1‑65 (MQ=255)
                   ttGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAg                                                      >  1:1663468/1‑68 (MQ=255)
                    tgtgTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcgc                                                  >  1:4724543/1‑71 (MQ=255)
                    tgtgTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcg                                                   >  1:4044360/1‑70 (MQ=255)
                    tgtgTCGTTATGGGAGTAGGCACTCATGTAATCATCATTATGTTTGATGATAATCATATTAACg                                                         >  1:5221032/1‑64 (MQ=255)
                     gtgtCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcgcg                                                 >  1:3937281/1‑71 (MQ=255)
                     gtgtCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcgcg                                                 >  1:1290134/1‑71 (MQ=255)
                     gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGccac                                                  <  1:3115612/70‑3 (MQ=255)
                     gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCcacg                                                 >  1:2447017/1‑71 (MQ=255)
                     gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCcacg                                                 >  1:865233/1‑71 (MQ=255)
                     gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCcacg                                                 >  1:3550347/1‑71 (MQ=255)
                     gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACAGTAGCcacg                                                 >  1:5366931/1‑71 (MQ=255)
                      tgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATgaga                                                               >  1:744114/1‑56 (MQ=255)
                      tgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATgaga                                                               >  1:3683978/1‑56 (MQ=255)
                        tcattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAg                                                      >  1:832223/4‑63 (MQ=255)
                        tcattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCa                                               <  1:4151842/67‑1 (MQ=255)
                        tcattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCa                                               <  1:147037/67‑1 (MQ=255)
                        tcattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAg                                              >  1:4608661/4‑71 (MQ=255)
                        tcattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAg                                              >  1:4552884/4‑71 (MQ=255)
                        tcattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAg                                              >  1:435145/4‑71 (MQ=255)
                        tCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAg                                                                <  1:5070809/53‑1 (MQ=255)
                         cattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTAcc                                                          >  1:1162839/3‑58 (MQ=255)
                         cattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAg                                              >  1:5488917/3‑70 (MQ=255)
                          gTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgcg                                            >  1:279509/1‑71 (MQ=255)
                          gTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgcg                                            >  1:171463/1‑71 (MQ=255)
                          gTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAga                                             >  1:2485721/1‑69 (MQ=255)
                          attATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGAt                                                              >  1:2727593/2‑53 (MQ=255)
                          attATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTAcc                                                          >  1:3314409/2‑57 (MQ=255)
                          attATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGccac                                                  <  1:2958374/64‑3 (MQ=255)
                          attATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgct                                            <  1:4579251/70‑2 (MQ=255)
                           ttATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgctg                                           <  1:2302439/71‑3 (MQ=255)
                           ttATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgctg                                           <  1:4749363/71‑3 (MQ=255)
                           ttATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgctg                                           <  1:5610259/71‑3 (MQ=255)
                           taaTGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCa                                               <  1:1184967/65‑1 (MQ=255)
                            tATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgcgc                                           >  1:4100491/1‑70 (MQ=255)
                              tGGGCGTAGGCACTCAGGTAATCATCATTTTGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGt                                         >  1:1414191/1‑70 (MQ=255)
                              tGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGtt                                        <  1:2870711/71‑1 (MQ=255)
                               gggCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATaa                                                                    >  1:3229625/1‑42 (MQ=255)
                               gggCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAAt                                                                   >  1:4294273/1‑43 (MQ=255)
                               gggCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAAt                                                                   >  1:4626672/1‑43 (MQ=255)
                                 gCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat                                                                      <  1:1852372/38‑1 (MQ=255)
                                    tAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCa                                                                 <  1:1068833/40‑1 (MQ=255)
                                      ggCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCa                                 >  1:2287529/1‑70 (MQ=255)
                                       gCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATa                               <  1:1831448/71‑1 (MQ=255)
                                          cTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTAcc                                     <  1:4423282/62‑1 (MQ=255)
                                          cTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                            <  1:2100593/71‑1 (MQ=255)
                                            cAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaa                          >  1:4076349/1‑71 (MQ=255)
                                              ggTAATCATCATTATGTTTGATGATAATCAGATTACCGTAg                                                      <  1:1838597/41‑1 (MQ=255)
                                              ggTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAAcg                        >  1:1033989/1‑71 (MQ=255)
                                               gTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATa                               <  1:893411/63‑1 (MQ=255)
                                                 aaTCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                            >  1:775615/1‑64 (MQ=255)
                                                 aaTCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                            >  1:5385861/1‑64 (MQ=255)
                                                 aaTCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                            >  1:2230364/1‑64 (MQ=255)
                                                    catcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAg                                              >  1:3969487/1‑43 (MQ=255)
                                                    catcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAg                                              >  1:3919408/1‑43 (MQ=255)
                                                    catcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAg                                              >  1:2394416/1‑43 (MQ=255)
                                                    catcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                            <  1:1029486/61‑1 (MQ=255)
                                                     atcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaa                          >  1:2835029/1‑62 (MQ=255)
                                                     atcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCAt                 <  1:43273/71‑1 (MQ=255)
                                                        atTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCt               <  1:5115252/70‑1 (MQ=255)
                                                        atTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCt               <  1:5096701/70‑1 (MQ=255)
                                                         ttATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATaaa                             >  1:4966377/1‑55 (MQ=255)
                                                            tGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaa                          >  1:4051478/1‑55 (MQ=255)
                                                             gTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTc         >  1:2657380/1‑71 (MQ=255)
                                                              tttGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTa                                       <  1:4623498/40‑1 (MQ=255)
                                                                tgatgaTAATCAGATTACCGTAGCCGCGCAGCGCGTTAcc                                     >  1:2486644/1‑40 (MQ=255)
                                                                tgatgaTAATCAGATTACCGTAGCCGCGCAGCGCGTTAcc                                     >  1:3657265/1‑40 (MQ=255)
                                                                tgatgaTAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAg                                   >  1:1107523/1‑42 (MQ=255)
                                                                tgatgaTAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTcgc       <  1:4085857/70‑1 (MQ=255)
                                                                  atgatAATCAGATTACCGTAGCCGCGCAGCGCGTTAc                                      >  1:2184482/1‑37 (MQ=255)
                                                                  atgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCg            <  1:5142663/63‑1 (MQ=255)
                                                                  atgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGAt    <  1:4127558/71‑1 (MQ=255)
                                                                  atgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGAt    <  1:1269571/71‑1 (MQ=255)
                                                                   tgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGAt    <  1:2406421/70‑1 (MQ=255)
                                                                    gatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATaa  >  1:2650116/1‑71 (MQ=255)
                                                                    |                                                                      
TGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAA  >  minE/1837406‑1837544

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: