Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 822,278 C→T 7.0% intergenic (‑254/+150) 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,2780CT 7.0% 222.1 / 4.8 71intergenic (‑254/+150)ldrC/chaAtoxic polypeptide, small/calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base T (3/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

GTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTAACGA  >  minE/822208‑822343
                                                                      |                                                                 
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)
            cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa                                                          >  1:184691/1‑68 (MQ=40)
            cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa                                                          >  1:2580893/1‑68 (MQ=40)
            cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAg                                                         >  1:1772258/1‑69 (MQ=40)
            cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt                                                       <  1:3704765/71‑1 (MQ=39)
            cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt                                                       <  1:3229484/71‑1 (MQ=39)
              tCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTg                                                      <  1:3484285/70‑1 (MQ=40)
              tCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGt                                                     <  1:1046257/71‑1 (MQ=39)
                    ttCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGt                                               <  1:530263/71‑1 (MQ=34)
                     tCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAg                                                         >  1:1849298/1‑60 (MQ=11)
                     tCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAg                                                         >  1:1907939/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:2167625/70‑1 (MQ=255)
                                     ccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgtg                               <  1:718169/70‑1 (MQ=255)
                                      cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGtgtg                             <  1:3909709/71‑1 (MQ=255)
                                         aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTgg                            <  1:1173134/69‑1 (MQ=255)
                                         aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt                           <  1:151232/70‑1 (MQ=255)
                                         aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt                           <  1:617587/70‑1 (MQ=255)
                                         aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTa                          <  1:1221587/71‑1 (MQ=255)
                                         aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTa                          >  1:1960187/1‑71 (MQ=255)
                                         aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTa                          >  1:3236550/1‑71 (MQ=255)
                                          aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                                    >  1:222464/1‑44 (MQ=38)
                                          aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                                    >  1:421502/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)
                                                              cATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTa      <  1:2997317/70‑1 (MQ=255)
                                                                ttCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTAc                         >  1:1845729/1‑49 (MQ=255)
                                                                ttCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGCGTGGTACAGGCTCGCACTTCTGATTaa     >  1:1687952/1‑69 (MQ=255)
                                                                 tCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTAc                         >  1:323578/1‑48 (MQ=255)
                                                                  ccGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTAACGa  >  1:3833384/1‑70 (MQ=255)
                                                                  ccGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTAACGa  >  1:620519/1‑70 (MQ=255)
                                                                     gCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGc                     <  1:3046460/48‑1 (MQ=255)
                                                                      |                                                                 
GTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTAACGA  >  minE/822208‑822343

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: