Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 822,269 A→G 8.1% intergenic (‑245/+159) ldrC ← / ← chaA toxic polypeptide, small/calcium/sodium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE822,2690AG 8.1% 205.8 / 7.8 74intergenic (‑245/+159)ldrC/chaAtoxic polypeptide, small/calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (32/36);  new base G (3/3);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

CGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATT  >  minE/822201‑822338
                                                                    |                                                                     
cGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACa                                                                     >  1:1808010/1‑71 (MQ=39)
 ggCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACAt                                                                    <  1:168253/71‑1 (MQ=39)
   cTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTc                                                                  <  1:1373289/71‑1 (MQ=39)
    ttGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTcc                                                                 >  1:819932/1‑71 (MQ=39)
      ggTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTcc                                                                 <  1:458069/69‑1 (MQ=40)
      ggTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCgg                                                               >  1:2280375/1‑71 (MQ=39)
       gtgtGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGc                                                              <  1:1834663/71‑1 (MQ=39)
            taTACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCtttt                                                          <  1:1831327/69‑1 (MQ=255)
            gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc                                                         <  1:2201998/71‑1 (MQ=39)
             aTACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCa                                                        <  1:3503504/71‑1 (MQ=39)
              tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa                                                       <  1:582289/71‑1 (MQ=39)
                ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa                                                     <  1:3813511/71‑1 (MQ=39)
                  acGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                           >  1:3703888/1‑63 (MQ=40)
                  acGTCGCAATTCATCTTAACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                           >  1:153315/1‑63 (MQ=255)
                   cGTCGCAATTCATCTTCACTTCATCCGTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt                                                  >  1:3787728/1‑71 (MQ=255)
                   cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCgg                                                               >  1:366524/1‑58 (MQ=255)
                   cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCgg                                                               >  1:2636052/1‑58 (MQ=255)
                   cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa                                                     >  1:2580893/1‑68 (MQ=40)
                   cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa                                                     >  1:184691/1‑68 (MQ=40)
                   cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAg                                                    >  1:1772258/1‑69 (MQ=40)
                   cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt                                                  <  1:3229484/71‑1 (MQ=39)
                   cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt                                                  <  1:3704765/71‑1 (MQ=39)
                     tCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTg                                                 <  1:3484285/70‑1 (MQ=40)
                     tCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGt                                                <  1:1046257/71‑1 (MQ=39)
                           ttCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTGCATTCCggtt                                                             <  1:2703103/52‑3 (MQ=255)
                           ttCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGt                                          <  1:530263/71‑1 (MQ=34)
                            tCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAg                                                    >  1:1907939/1‑60 (MQ=11)
                            tCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAg                                                    >  1:1849298/1‑60 (MQ=11)
                                 ttcacttcaTCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGc                                                   <  1:193888/56‑1 (MQ=11)
                                   cacttcaTCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCaa                                                       <  1:1492259/50‑1 (MQ=11)
                                   cacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGttgc                                    <  1:3203528/69‑3 (MQ=255)
                                   cacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGtt                                      <  1:1881329/67‑1 (MQ=40)
                                   cacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTc                                    <  1:1308077/69‑1 (MQ=40)
                                      ttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTg                                              <  1:3249432/56‑1 (MQ=35)
                                       tcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGtt                                      <  1:2846250/63‑1 (MQ=40)
                                        caTCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGtta                                              >  1:2533098/1‑53 (MQ=11)
                                        caTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa                                                       >  1:2377031/1‑45 (MQ=38)
                                        caTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAg                                                    >  1:1580054/1‑48 (MQ=38)
                                        caTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                             >  1:1188287/1‑71 (MQ=255)
                                            ccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgtgt                         <  1:3540193/71‑1 (MQ=255)
                                            ccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgtgt                         >  1:3490891/1‑71 (MQ=255)
                                            ccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgtg                          <  1:718169/70‑1 (MQ=255)
                                            ccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgtg                          <  1:2167625/70‑1 (MQ=255)
                                             cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGtgtg                        <  1:3909709/71‑1 (MQ=255)
                                                aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTgg                       <  1:1173134/69‑1 (MQ=255)
                                                aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt                      <  1:617587/70‑1 (MQ=255)
                                                aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt                      <  1:151232/70‑1 (MQ=255)
                                                aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTa                     >  1:3236550/1‑71 (MQ=255)
                                                aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTa                     >  1:1960187/1‑71 (MQ=255)
                                                aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTa                     <  1:1221587/71‑1 (MQ=255)
                                                 aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                               >  1:421502/1‑44 (MQ=38)
                                                 aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                               >  1:222464/1‑44 (MQ=38)
                                                  aTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACa                   <  1:3769323/71‑1 (MQ=255)
                                                     tttGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGc                >  1:1706967/1‑71 (MQ=255)
                                                      ttGCAATGAAGAGTACATTCCGGCTTTTTAACAGCTGTTGTAGTGGTTTCGTGCAGGGTGTTGTACAGGCt               >  1:2992040/1‑71 (MQ=21)
                                                       tGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                               >  1:1011061/1‑38 (MQ=40)
                                                       tGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTc              >  1:3688529/1‑71 (MQ=255)
                                                       tGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTc              >  1:2267544/1‑71 (MQ=255)
                                                         cAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAg                                           <  1:797970/40‑1 (MQ=40)
                                                          aaTGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                           >  1:1393725/1‑55 (MQ=255)
                                                          aaTGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTgg                       >  1:3261759/1‑59 (MQ=255)
                                                          aaTGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGc            >  1:3749856/1‑70 (MQ=255)
                                                          aaTGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGc            >  1:3857699/1‑70 (MQ=255)
                                                          aaTGAAGAGTACATTCCGCCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTgg                       >  1:1631771/1‑59 (MQ=255)
                                                           aTGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTAc                    <  1:2618971/61‑1 (MQ=255)
                                                            tGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAAt         <  1:263028/71‑1 (MQ=255)
                                                              aaGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTg                                 >  1:3129859/1‑45 (MQ=255)
                                                              aaGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgtg                          <  1:525217/52‑1 (MQ=255)
                                                                 agTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGa    <  1:2097700/71‑1 (MQ=255)
                                                                  gTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                           >  1:1829139/1‑47 (MQ=255)
                                                                   tACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAg                  <  1:1185521/55‑1 (MQ=255)
                                                                   tACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAg                  <  1:6262/55‑1 (MQ=255)
                                                                   tACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGAtt  <  1:3868239/71‑1 (MQ=255)
                                                                   tACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTATCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGAt   >  1:500221/1‑70 (MQ=255)
                                                                    |                                                                     
CGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATT  >  minE/822201‑822338

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: