Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 822,278 C→T 11.9% 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,2780CT11.9% 232.7 / 15.1 84intergenic (‑254/+150)ldrC/chaAtoxic polypeptide, small/calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base C (37/37);  new base T (8/2);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 9.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTAACGACC  >  minE/822211‑822345
                                                                   |                                                                   
tgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCt                                                                    <  1:4948895/69‑1 (MQ=40)
tgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCt                                                                    <  1:3742314/69‑1 (MQ=40)
tgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCt                                                                    <  1:1886170/69‑1 (MQ=40)
  gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGTCTTTTc                                                                <  1:3874844/71‑1 (MQ=255)
  gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCtttt                                                                 <  1:1141021/70‑1 (MQ=40)
  gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc                                                                <  1:3260012/71‑1 (MQ=39)
  gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc                                                                <  1:1493152/71‑1 (MQ=39)
  gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc                                                                <  1:2363604/71‑1 (MQ=39)
    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa                                                              <  1:2968349/71‑1 (MQ=39)
     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:3071989/1‑70 (MQ=35)
              cAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTg                                                     <  1:1033075/70‑1 (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)
                   cATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAAc                                                             >  1:2537834/1‑57 (MQ=35)
                    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)
                         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:4011666/1‑71 (MQ=255)
                          acttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTg                                        >  1:4098301/1‑71 (MQ=255)
                          acttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTg                                        >  1:707559/1‑71 (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:2970997/1‑71 (MQ=255)
                              caTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                                    >  1:1634326/1‑71 (MQ=255)
                              caTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                                    >  1:64072/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:3545259/1‑71 (MQ=255)
                                      aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTa                            >  1:2602757/1‑71 (MQ=255)
                                      aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTa                            >  1:2646107/1‑71 (MQ=255)
                                       aaTGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAg                                                           >  1:303373/1‑39 (MQ=14)
                                              gCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGc                       <  1:4691385/68‑1 (MQ=255)
                                              gCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGc                       <  1:560414/68‑1 (MQ=255)
                                              gCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGc                       <  1:1179779/68‑1 (MQ=255)
                                               cAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCg                    >  1:3446764/1‑70 (MQ=255)
                                               cAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCg                    <  1:3172830/70‑1 (MQ=255)
                                               cAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGc                   <  1:2235729/71‑1 (MQ=255)
                                               cAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGc                   <  1:514652/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:3305654/71‑1 (MQ=255)
                                                         tACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGAtt         <  1:438125/71‑1 (MQ=255)
                                                           cATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCg                                          >  1:4995691/1‑36 (MQ=255)
                                                           cATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCg                                          >  1:5628679/1‑36 (MQ=255)
                                                           cATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCg                                          >  1:637144/1‑36 (MQ=255)
                                                           cATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAAtt               >  1:1642168/1‑63 (MQ=255)
                                                           cATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGa           >  1:116961/1‑67 (MQ=255)
                                                           cATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTaa       >  1:1539516/1‑71 (MQ=255)
                                                           cATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTaa       >  1:4928501/1‑71 (MQ=255)
                                                           cATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTaa       >  1:3947769/1‑71 (MQ=255)
                                                           cATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTaa       >  1:4032554/1‑71 (MQ=255)
                                                             ttCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGt                                         <  1:3195176/35‑1 (MQ=255)
                                                              tCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTg                                        <  1:4687405/35‑1 (MQ=255)
                                                               ccGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGa           >  1:4094762/1‑63 (MQ=255)
                                                               ccGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTAACGAc   >  1:3920912/1‑71 (MQ=255)
                                                                cGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTAACGAcc  <  1:5567393/71‑1 (MQ=255)
                                                                   |                                                                   
TGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCAATTCTGATTAACGACC  >  minE/822211‑822345

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: