Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,837,452 G→T 12.3% L323M (CTG→ATG) ‡ nlpD ← predicted outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,837,4520GT12.3% 224.4 / 9.9 72L323M (CTG→ATG) ‡nlpDpredicted outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (38/25);  new base T (2/7);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 6.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

ATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACG  >  minE/1837382‑1837522
                                                                      |                                                                      
aTTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAg                                                                        >  1:2241295/1‑71 (MQ=255)
aTTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAg                                                                        >  1:1856763/1‑71 (MQ=255)
  tttttGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAgg                                                                       >  1:1710542/1‑70 (MQ=255)
       gCCCCGCCTTAACTTCTAGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTaa                                                                    <  1:3050616/68‑1 (MQ=255)
        ccccGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAatca                                                                 <  1:640369/70‑1 (MQ=255)
            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:2712619/1‑71 (MQ=255)
                 aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                                       >  1:3200491/1‑71 (MQ=255)
                 aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                                       >  1:2833653/1‑71 (MQ=255)
                 aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                                       >  1:1694774/1‑71 (MQ=255)
                 aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                                       >  1:2198573/1‑71 (MQ=255)
                   cttcttGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAg                                                                        >  1:3804242/1‑52 (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: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)
                            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)
                               cccGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTaa                                                                    >  1:964700/1‑44 (MQ=255)
                               cccGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTa                                                                     >  1:2299161/1‑43 (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)
                                             gtgtCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGttt                                                      >  1:2922585/1‑44 (MQ=255)
                                             gtgtCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTtgat                                                   >  1:1564132/1‑47 (MQ=255)
                                             gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATgaga                                         >  1:618454/1‑57 (MQ=255)
                                             gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATgaga                                         >  1:1107573/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:3853165/70‑2 (MQ=255)
                                                  attATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgct                      <  1:218434/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:3650222/1‑71 (MQ=255)
                                                     aTGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGt                   >  1:3501462/1‑71 (MQ=255)
                                                     aTGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGt                   >  1:1492891/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:2809451/70‑1 (MQ=255)
                                                                     aGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaa    <  1:3012238/70‑1 (MQ=255)
                                                                      ggTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAAcg  <  1:1719454/71‑1 (MQ=255)
                                                                      |                                                                      
ATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACG  >  minE/1837382‑1837522

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: