Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 822,269 A→G 18.7% 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,2690AG18.7% 205.8 / 35.1 91intergenic (‑245/+159)ldrC/chaAtoxic polypeptide, small/calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (37/37);  new base G (12/5);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATT  >  minE/822200‑822338
                                                                     |                                                                     
ccGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTa                                                                       >  1:4998389/1‑70 (MQ=40)
    cTTGGTGTGGATACTACGTCGCAATTCATCTGCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTc                                                                  <  1:295500/71‑1 (MQ=255)
      tGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCg                                                                >  1:755966/1‑71 (MQ=39)
      tGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCg                                                                >  1:4335476/1‑71 (MQ=39)
      tGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCg                                                                >  1:4329109/1‑71 (MQ=39)
        gtgtGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCgg                                                               <  1:5155005/70‑1 (MQ=40)
        gtgtGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCgg                                                               <  1:2021548/70‑1 (MQ=40)
           tgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCt                                                             <  1:4948895/69‑1 (MQ=40)
           tgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCt                                                             <  1:1886170/69‑1 (MQ=40)
           tgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCt                                                             <  1:3742314/69‑1 (MQ=40)
             gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGTCTTTTc                                                         <  1:3874844/71‑1 (MQ=255)
             gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCtttt                                                          <  1:1141021/70‑1 (MQ=40)
             gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc                                                         <  1:1493152/71‑1 (MQ=39)
             gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc                                                         <  1:3260012/71‑1 (MQ=39)
             gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc                                                         <  1:2363604/71‑1 (MQ=39)
               tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa                                                       <  1:2968349/71‑1 (MQ=39)
                actacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACa                                                                     >  1:2034085/1‑56 (MQ=255)
                actacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACa                                                                     >  1:360394/1‑56 (MQ=255)
                actacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACa                                                                     >  1:2261381/1‑56 (MQ=255)
                actacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACa                                                                     >  1:2521769/1‑56 (MQ=255)
                actacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa                                                       >  1:3563197/1‑70 (MQ=40)
                  tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa                                                       >  1:2006096/1‑68 (MQ=40)
                      actcAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGt                                                >  1:4185514/4‑71 (MQ=11)
                       cGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                               >  1:2606835/1‑71 (MQ=38)
                         cAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTg                                              <  1:1033075/70‑1 (MQ=35)
                         cAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTg                                              >  1:3071989/1‑70 (MQ=35)
                         cAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGc                                             <  1:2307039/71‑1 (MQ=34)
                            ttCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCt                                                             >  1:5149153/1‑52 (MQ=35)
                            ttCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc                                                         <  1:544924/56‑1 (MQ=35)
                            ttCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc                                                         <  1:454445/56‑1 (MQ=35)
                             tCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTg                                         >  1:1470259/1‑71 (MQ=35)
                              cATCTTCACTTCATTCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTaaa                                                       >  1:493176/1‑56 (MQ=11)
                              cATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTGCa                                                                     >  1:4921806/1‑42 (MQ=255)
                              cATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTGCa                                                                     >  1:4338898/1‑42 (MQ=255)
                              cATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTGCa                                                                     >  1:629988/1‑42 (MQ=255)
                              cATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTAc                                                                      >  1:1247510/1‑41 (MQ=38)
                              cATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTAc                                                                      >  1:3738006/1‑41 (MQ=38)
                              cATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTAc                                                                      >  1:3378411/1‑41 (MQ=38)
                              cATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAAc                                                      >  1:2537834/1‑57 (MQ=35)
                               aTCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTGCATTc                                                                  >  1:4798178/1‑44 (MQ=255)
                               aTCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCt                                                  >  1:4791674/1‑60 (MQ=11)
                               aTCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGtca                                        <  1:659996/70‑3 (MQ=11)
                                  ttcacttcaTCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCaa                                                       >  1:4626209/1‑52 (MQ=11)
                                  ttcacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa                                                       >  1:3095258/1‑52 (MQ=35)
                                    cacttcaTCCCTGAAATGTTTGCAATGAAGAGTACa                                                                     >  1:4065603/1‑36 (MQ=40)
                                    cacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGttt                                     <  1:2063695/68‑1 (MQ=40)
                                    cacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGttt                                     <  1:2890281/68‑1 (MQ=40)
                                    cacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTc                                    <  1:1827659/69‑1 (MQ=40)
                                    cacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGt                                  >  1:5077407/1‑71 (MQ=255)
                                     acttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTg                                 >  1:707559/1‑71 (MQ=255)
                                     acttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTg                                 >  1:4011666/1‑71 (MQ=255)
                                     acttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTg                                 >  1:4098301/1‑71 (MQ=255)
                                      cttcaTCCCTGAAATGTTTGCAATGAAGAGTGCATTCCggt                                                              <  1:4814080/41‑2 (MQ=255)
                                      cttcaTCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGttt                                                            <  1:4099667/43‑1 (MQ=14)
                                       ttcaTCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTc                                                         <  1:40448/45‑1 (MQ=11)
                                       ttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc                                                         <  1:977655/45‑1 (MQ=38)
                                       ttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGc                                                   <  1:5588532/51‑1 (MQ=35)
                                       ttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGc                                                   <  1:2827122/51‑1 (MQ=35)
                                         caTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGTCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                             >  1:4729365/1‑71 (MQ=255)
                                         caTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCtttt                                                          >  1:2050910/1‑42 (MQ=38)
                                         caTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                             >  1:64072/1‑71 (MQ=255)
                                         caTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                             >  1:1634326/1‑71 (MQ=255)
                                         caTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                             >  1:2970997/1‑71 (MQ=255)
                                         caTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAg                              <  1:4052352/70‑1 (MQ=255)
                                          aTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTg                                              <  1:4285390/53‑1 (MQ=35)
                                             ccTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCaa                                                       >  1:2570016/1‑41 (MQ=14)
                                               tGAAATGTTTGCAATGAAGAGTGCATTCCGGtttt                                                           <  1:2813744/35‑1 (MQ=14)
                                                 aaaTGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTg                                                 >  1:2999634/1‑43 (MQ=14)
                                                 aaaTGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTg                                                 >  1:5420033/1‑43 (MQ=14)
                                                 aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAg                              >  1:2517851/1‑62 (MQ=255)
                                                 aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTa                     >  1:2602757/1‑71 (MQ=255)
                                                 aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTa                     >  1:2646107/1‑71 (MQ=255)
                                                 aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTa                     >  1:3545259/1‑71 (MQ=255)
                                                  aaTGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAg                                                    >  1:303373/1‑39 (MQ=14)
                                                         gCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGc                <  1:4691385/68‑1 (MQ=255)
                                                         gCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGc                <  1:1179779/68‑1 (MQ=255)
                                                         gCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGc                <  1:560414/68‑1 (MQ=255)
                                                          cAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCg             <  1:3172830/70‑1 (MQ=255)
                                                          cAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCg             >  1:3446764/1‑70 (MQ=255)
                                                          cAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGc            <  1:514652/71‑1 (MQ=255)
                                                          cAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGc            <  1:2235729/71‑1 (MQ=255)
                                                            aTGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTgg                       <  1:2068127/58‑1 (MQ=255)
                                                             tGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGtt                                      >  1:1563116/1‑42 (MQ=255)
                                                             tGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGc                                >  1:1322526/1‑48 (MQ=255)
                                                             tGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAAt         >  1:3207426/1‑71 (MQ=255)
                                                                agagTAGATTCCGGCTTTTCAACAGCTGTTGCAGTGGttt                                     <  1:2549572/40‑1 (MQ=255)
                                                                 gagTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCg                                   >  1:4536017/1‑41 (MQ=40)
                                                                 gagTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACCGGCTCGCAATTCTg     <  1:5105947/71‑1 (MQ=255)
                                                                 gagTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTg     <  1:5352516/71‑1 (MQ=255)
                                                                    tACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCaa          <  1:5012879/63‑1 (MQ=255)
                                                                    tACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGAtt  <  1:438125/71‑1 (MQ=255)
                                                                    tACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGAtt  <  1:3305654/71‑1 (MQ=255)
                                                                     |                                                                     
CCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATT  >  minE/822200‑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: