Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,925,766:2 +C 25.2% intergenic (‑102/+58) holC ← / ← pepA DNA polymerase III, chi subunit/aminopeptidase A, a cyteinylglycinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,925,7662.C25.2% 200.3 / 142.9 111intergenic (‑102/+58)holC/pepADNA polymerase III, chi subunit/aminopeptidase A, a cyteinylglycinase
Reads supporting (aligned to +/- strand):  ref base . (42/41);  new base C (18/10);  total (60/51)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 113 total aligned reads displayed.

AGGTATACACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCG  >  minE/2925701‑2925833
                                                                     |                                                                   
aGGTATACACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTa                                                                  >  1:1670600/1‑69 (MQ=255)
aGGTATACACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAgg                                                                >  1:3317544/1‑71 (MQ=255)
aGGTATACACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAgg                                                                >  1:5562827/1‑71 (MQ=255)
aGGTATACACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAgg                                                                >  1:3846808/1‑71 (MQ=255)
    ataCACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCg                                                             >  1:1273586/1‑70 (MQ=255)
     taCACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGa                                                           <  1:210658/71‑1 (MQ=255)
     taCACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGa                                                           <  1:767328/71‑1 (MQ=255)
       cacCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGAt                                                          <  1:1314499/70‑1 (MQ=255)
       cacCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGAt                                                          >  1:5343328/1‑70 (MQ=255)
       cacCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGAt                                                          >  1:1518991/1‑70 (MQ=255)
           gCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGc                                                     <  1:428961/71‑1 (MQ=255)
            cAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCg                                                    <  1:3880629/71‑1 (MQ=255)
            cAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCg                                                    <  1:5350910/71‑1 (MQ=255)
              gctggaattgcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCg                                                    <  1:423301‑M2/32‑1 (MQ=255)
              gcgcCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGtt                                                  <  1:1483915/71‑1 (MQ=255)
              gcgcCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGtt                                                  <  1:726348/71‑1 (MQ=255)
               cgcCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGtt                                                  >  1:945448/1‑70 (MQ=255)
                gcCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑Ctgt‑‑g                                                                    <  1:261892/51‑1 (MQ=255)
                gcCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCa                                                >  1:2894623/1‑71 (MQ=255)
                 ggaattgcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGc                                                     <  1:2998829‑M2/31‑1 (MQ=255)
                 ggaattgcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGtt                                                  >  1:4305531‑M2/37‑70 (MQ=255)
                 ggaattgcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGtt                                                  >  1:1703141‑M2/37‑70 (MQ=255)
                 ccTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGAt                                                          <  1:1373728/60‑1 (MQ=255)
                  gaattgcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGtt                                                  <  1:1040741‑M2/34‑1 (MQ=255)
                  gaattgcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGtt                                                  <  1:2241225‑M2/34‑1 (MQ=255)
                  gaattgcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCa                                                >  1:2826639‑M2/36‑71 (MQ=255)
                  cTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCAcg                                              >  1:4238028/1‑71 (MQ=255)
                   tGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCAc                                               <  1:2914178/69‑1 (MQ=255)
                   aattgcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGtt                                                  <  1:4719136‑M2/34‑1 (MQ=255)
                     ttgcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCAcgc                                             >  1:5025649‑M2/33‑71 (MQ=255)
                     ttgcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCAcgc                                             >  1:1888675‑M2/33‑71 (MQ=255)
                     aTGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCAcgccg                                           <  1:4559033/71‑1 (MQ=255)
                     aTGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCAcgccg                                           >  1:1721923/1‑71 (MQ=255)
                     aTGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCAcgccg                                           >  1:5057695/1‑71 (MQ=255)
                     aTGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCAcgc                                             >  1:491877/1‑69 (MQ=255)
                     aTGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCAcgc                                             >  1:1318144/1‑69 (MQ=255)
                      tgcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCAcgcc                                            >  1:4708630‑M2/32‑71 (MQ=255)
                      tgcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCAcgcc                                            >  1:4258136‑M2/32‑71 (MQ=255)
                      tGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGAt                                                          >  1:1858267/1‑55 (MQ=255)
                       gcctgatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCAcgccg                                           >  1:4888826‑M2/31‑71 (MQ=255)
                        cGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCa                                         <  1:4922827/70‑1 (MQ=255)
                        cGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCa                                         <  1:1819582/70‑1 (MQ=255)
                        cGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCAt                                        >  1:2017604/1‑71 (MQ=255)
                           gatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGtt                                                  <  1:5394254‑M2/34‑1 (MQ=255)
                           gatgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCACGCCGCATc                                       >  1:4182421‑M2/27‑71 (MQ=255)
                            gCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑Ctgt‑‑g                                                                    >  1:359141/1‑39 (MQ=255)
                            gCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATaa                                                        >  1:2659578/1‑51 (MQ=255)
                            gCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATaa                                                        >  1:2968987/1‑51 (MQ=255)
                            gCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATa                                                         >  1:5563422/1‑50 (MQ=255)
                            gCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATa                                                         >  1:4471487/1‑50 (MQ=255)
                            atgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                      >  1:3107430‑M2/26‑71 (MQ=255)
                            atgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                      >  1:1673904‑M2/26‑71 (MQ=255)
                            atgcgctacgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCACGCCGCATc                                       >  1:3474405‑M2/26‑70 (MQ=255)
                             cTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGc                                   <  1:4423034/71‑1 (MQ=255)
                              tgttgtGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGc                                   >  1:1607875/1‑70 (MQ=255)
                               gttgtGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCaa                                 >  1:5329858/1‑71 (MQ=255)
                                ttgtGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAt                                <  1:4786068/71‑1 (MQ=255)
                                   acgcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCACGCCGCATCCGgcata                                >  1:2012476‑M2/19‑68 (MQ=255)
                                    gtCTTATCATGCCTACGAGTTCGATGCTCTGT‑‑GTAGGCCGGATAAGGTGTTCACGCCGCATCCGGCattg                               <  1:1293923/70‑4 (MQ=9)
                                     gcttatcaggcctacgAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCACGCCGCATCCGgcata                                <  1:271830‑M2/52‑3 (MQ=255)
                                       ttatcaggcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATaa                                                        >  1:1171725‑M2/15‑40 (MQ=255)
                                       ttatcaggcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATaa                                                        >  1:5583384‑M2/15‑40 (MQ=255)
                                       ttatcaggcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATaa                                                        >  1:4212871‑M2/15‑40 (MQ=255)
                                         aTCAGGCCTACGAGTTCAGTG‑‑CTCTACGTAGGCCGGATaa                                                        >  1:3305296/1‑40 (MQ=16)
                                         aTCAGGCCTACGAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGtt                                                  <  1:1028943/46‑1 (MQ=12)
                                         aTCAGGCCTACGAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGtt                                                  <  1:4962770/46‑1 (MQ=12)
                                          tcaggcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGc                                   <  1:765018‑M2/47‑1 (MQ=255)
                                          tcaggcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGc                                   <  1:2927832‑M2/47‑1 (MQ=255)
                                          tcaggcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTg                      <  1:1496778‑M2/60‑1 (MQ=255)
                                          tcaggcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTg                      <  1:3019297‑M2/60‑1 (MQ=255)
                                           caggcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGc                     >  1:1660041‑M2/11‑71 (MQ=255)
                                           caggcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGc                     >  1:5031944‑M2/11‑71 (MQ=255)
                                           cAGGCCTACGAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGtt                                                  <  1:349833/44‑1 (MQ=12)
                                           cAGGCCTACGAGTTCAGTG‑‑CTCTACGTAGGCCGGATAAGGCGTTCa                                                >  1:2941399/1‑46 (MQ=12)
                                            aggcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCAc                                               <  1:2827852‑M2/35‑1 (MQ=255)
                                             ggcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCt                   <  1:3217045‑M2/63‑1 (MQ=255)
                                              gcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTg                  <  1:3302510‑M2/64‑1 (MQ=255)
                                              gcctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTg                  <  1:270384‑M2/64‑1 (MQ=255)
                                                ctacgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGc                                   <  1:2745454‑M2/47‑1 (MQ=255)
                                                  acgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATaa                              >  1:4640862‑M2/4‑55 (MQ=255)
                                                  acgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATa                               >  1:4634188‑M2/4‑54 (MQ=255)
                                                  acgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACg               >  1:1725065‑M2/4‑70 (MQ=255)
                                                  acgAGTTCAGTG‑‑CTGT‑‑GTACGTCGGATAAGGCGTTCACGCCGCATCCGGCAATaa                              >  1:1313807‑M2/4‑55 (MQ=255)
                                                   cgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAt                                <  1:3222217‑M2/50‑1 (MQ=255)
                                                   cgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAt                                <  1:4216311‑M2/50‑1 (MQ=255)
                                                   cgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATaa                              >  1:2448353‑M2/3‑54 (MQ=255)
                                                   cgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATaa                              >  1:5196982‑M2/3‑54 (MQ=255)
                                                   cgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGc              <  1:4029648‑M2/68‑1 (MQ=255)
                                                   cgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGc              <  1:2146912‑M2/68‑1 (MQ=255)
                                                   cgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCa             <  1:3559362‑M2/69‑1 (MQ=255)
                                                   cgAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCCTGCCTGACGc              <  1:137436‑M2/68‑1 (MQ=255)
                                                    gAGTTCAGTG‑‑CTGT‑‑GAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAAt           >  1:2675517‑M2/2‑71 (MQ=255)
                                                       ttCAGTG‑‑CTGT‑‑CTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGc                          <  1:5680109/54‑1 (MQ=255)
                                                         cAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATa                               >  1:3135213/1‑47 (MQ=255)
                                                         cAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTAct       >  1:5107347/1‑71 (MQ=255)
                                                         cAGTG‑‑CTGT‑‑GTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTAc        <  1:2231166/70‑1 (MQ=255)
                                                          aGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTActc      <  1:5446806/71‑1 (MQ=255)
                                                          aGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTAct       <  1:2988778/70‑1 (MQ=255)
                                                           gTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTActct     <  1:5580577/71‑1 (MQ=255)
                                                                cTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCg  >  1:1452985/1‑71 (MQ=255)
                                                                     |                                                                   
AGGTATACACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTG‑‑CTGT‑‑GTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCG  >  minE/2925701‑2925833

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: