Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,837,444 G→A 10.7% A325A (GCC→GCT nlpD ← predicted outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,837,4440GA10.7% 220.6 / 13.5 75A325A (GCC→GCTnlpDpredicted outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (40/27);  new base A (2/6);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCA  >  minE/1837376‑1837513
                                                                    |                                                                     
gTCGATATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAgg                                                                      >  1:1301271/1‑70 (MQ=255)
   gCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAgg                                                                      >  1:1903038/1‑67 (MQ=255)
   gCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACt                                                                  <  1:1794422/71‑1 (MQ=255)
    cTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTc                                                                 <  1:3227877/71‑1 (MQ=255)
    cTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTc                                                                 <  1:212305/71‑1 (MQ=255)
      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:1694774/1‑71 (MQ=255)
                       aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                              >  1:3200491/1‑71 (MQ=255)
                       aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                              >  1:2833653/1‑71 (MQ=255)
                       aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                              >  1:2712619/1‑71 (MQ=255)
                       aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                              >  1:2198573/1‑71 (MQ=255)
                         cttcttGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACt                                                                  <  1:258571/49‑1 (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: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)
                                    tCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCa                                                                    >  1:1515143/1‑36 (MQ=255)
                                    tCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCa                                                                >  1:3789843/1‑40 (MQ=255)
                                     cccGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTaa                                                           >  1:964700/1‑44 (MQ=255)
                                     cccGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTa                                                            >  1:2299161/1‑43 (MQ=255)
                                       cGGACCAGCATTGTGTCGTTATGGGCGTAGGCACt                                                                  <  1:2448293/35‑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: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: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:1278284/1‑38 (MQ=255)
                                                        gTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt                                              >  1:874666/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:2535124/71‑1 (MQ=255)
                                                            tGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGtt         <  1:971080/71‑1 (MQ=255)
                                                              ggCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAAt                                    <  1:973303/42‑1 (MQ=255)
                                                                    ggCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCa  >  1:2378536/1‑70 (MQ=255)
                                                                    |                                                                     
GTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCA  >  minE/1837376‑1837513

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: