Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,443,8250TC36.4% 74.4 / 79.7 77K437K (AAA→AAGilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base T (33/16);  new base C (0/28);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 7.72e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAG  >  minE/2443764‑2443891
                                                             |                                                                  
ggCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATg                                                           >  1:995224/1‑71 (MQ=255)
ggCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATg                                                           >  1:468400/1‑71 (MQ=255)
ggCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATg                                                           >  1:2231432/1‑71 (MQ=255)
ggCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATg                                                           >  1:4726078/1‑71 (MQ=255)
ggCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATg                                                           >  1:3595796/1‑71 (MQ=255)
ggCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATg                                                           >  1:3686239/1‑71 (MQ=255)
     tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAg                                                        <  1:77511/69‑1 (MQ=255)
     tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAg                                                        <  1:3438269/69‑1 (MQ=255)
     tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAg                                                        <  1:58512/69‑1 (MQ=255)
     tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAg                                                        <  1:2452915/69‑1 (MQ=255)
     tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAg                                                        <  1:1697054/69‑1 (MQ=255)
     tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAg                                                        <  1:3789686/69‑1 (MQ=255)
     tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGcc                                                      >  1:1041811/1‑71 (MQ=255)
     tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGcc                                                      >  1:3066121/1‑71 (MQ=255)
     tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGcc                                                      >  1:2963267/1‑71 (MQ=255)
     tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGcc                                                      >  1:2940503/1‑71 (MQ=255)
     tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGcc                                                      >  1:4661111/1‑71 (MQ=255)
     tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGcc                                                      >  1:3748545/1‑71 (MQ=255)
                     tgccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa                                      <  1:2620903/70‑1 (MQ=255)
                     ggccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa                                      <  1:4649532/71‑1 (MQ=255)
                     ggccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa                                      <  1:2638649/71‑1 (MQ=255)
                     ggccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa                                      <  1:4454948/71‑1 (MQ=255)
                     ggccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa                                      <  1:4505832/71‑1 (MQ=255)
                      gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAg                                         >  1:588065/1‑67 (MQ=255)
                      gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAg                                         >  1:3572857/1‑67 (MQ=255)
                      gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAg                                         >  1:4290689/1‑67 (MQ=255)
                      gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAg                                         >  1:4498068/1‑67 (MQ=255)
                      gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGtag                                      >  1:4685860/1‑68 (MQ=255)
                      gccggAGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCa                                                         >  1:1696116/1‑51 (MQ=255)
                       ccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc                                    <  1:2960280/71‑1 (MQ=255)
                       ccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc                                    <  1:2706261/71‑1 (MQ=255)
                       ccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc                                    <  1:3093379/71‑1 (MQ=255)
                       ccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc                                    <  1:4613248/71‑1 (MQ=255)
                       ccggTGAATTTGAGGAAGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc                                    <  1:1682618/71‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:3570949/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:883823/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:826293/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:82374/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:3834846/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:78479/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:68449/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:4087833/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:4625977/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:4323995/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:525205/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:2313096/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:1097462/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:1157673/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:1257426/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:1277433/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:1545348/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:1587670/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:2199953/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:3789365/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:2329256/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:2481977/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:2652028/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:2871520/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:3599102/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:3282539/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:3015535/69‑1 (MQ=255)
                          gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa                                 <  1:3178736/69‑1 (MQ=255)
                            gAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTa                                 >  1:229448/1‑69 (MQ=255)
                            gAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTa                                 >  1:1626976/1‑69 (MQ=255)
                            gAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTa                                 >  1:798064/1‑69 (MQ=255)
                            gAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTa                                 >  1:1169459/1‑69 (MQ=255)
                            gAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTa                                 >  1:113447/1‑69 (MQ=255)
                                    ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTcc                                               >  1:3387701/1‑47 (MQ=255)
                                    ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGc                        >  1:2227309/1‑70 (MQ=255)
                                    ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGc                        >  1:4266210/1‑70 (MQ=255)
                                    ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGc                        >  1:4187558/1‑70 (MQ=255)
                                    ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGc                        >  1:3388268/1‑70 (MQ=255)
                                    ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGc                        >  1:4081607/1‑70 (MQ=255)
                                    ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGc                        >  1:3587132/1‑70 (MQ=255)
                                    ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCTAGAGCACCGc                        >  1:788576/1‑69 (MQ=255)
                                    ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAACTTACCg                                   >  1:4585782/1‑59 (MQ=255)
                                                         ccGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAg  >  1:1516587/1‑71 (MQ=255)
                                                             |                                                                  
GGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAG  >  minE/2443764‑2443891

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: