Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE460,6870GT39.3% 113.1 / 79.3 84T79T (ACC→ACAgltAcitrate synthase
Reads supporting (aligned to +/- strand):  ref base G (14/37);  new base T (0/33);  total (14/70)
Fisher's exact test for biased strand distribution p-value = 5.59e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATCAATAAAAG  >  minE/460618‑460751
                                                                     |                                                                
gTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCg                                                                  <  1:3893050/70‑1 (MQ=255)
gTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCg                                                                  <  1:4536433/70‑1 (MQ=255)
gTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCg                                                                  <  1:4850024/70‑1 (MQ=255)
 tCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGt                                                                <  1:1006650/71‑1 (MQ=255)
 tCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGt                                                                <  1:341734/71‑1 (MQ=255)
 tCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGt                                                                <  1:2802170/71‑1 (MQ=255)
 tCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGt                                                                <  1:4518723/71‑1 (MQ=255)
 tCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGt                                                                <  1:2325135/71‑1 (MQ=255)
 tCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGt                                                                <  1:1024438/71‑1 (MQ=255)
         ttCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAg                                                          <  1:3628715/69‑1 (MQ=255)
         ttCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAg                                                          <  1:2179323/69‑1 (MQ=255)
         ttCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAg                                                          <  1:3200425/69‑1 (MQ=255)
         ttCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAg                                                          <  1:4165846/69‑1 (MQ=255)
         ttCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAg                                                          <  1:3017261/69‑1 (MQ=255)
         ttCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAg                                                          <  1:2941349/69‑1 (MQ=255)
         ttCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAg                                                          <  1:2075907/69‑1 (MQ=255)
         ttCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAg                                                          <  1:2009439/69‑1 (MQ=255)
                                     agGATGTAACAAACTTGCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCaa                            <  1:2895807/71‑1 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcga                                                  >  1:2775213/1‑49 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcga                                                  >  1:281416/1‑49 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcga                                                  >  1:2809797/1‑49 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcga                                                  >  1:4027819/1‑49 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcga                                                  >  1:836456/1‑49 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcga                                                  >  1:270480/1‑49 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcga                                                  >  1:4405844/1‑49 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcga                                                  >  1:913579/1‑49 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcga                                                  >  1:1833985/1‑49 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcga                                                  >  1:1786772/1‑49 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcga                                                  >  1:668493/1‑49 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCaa                            <  1:905516/71‑1 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCaa                            <  1:2296151/71‑1 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCaa                            <  1:3504066/71‑1 (MQ=255)
                                     agGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCaa                            <  1:565414/71‑1 (MQ=255)
                                      ggATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCaa                            <  1:4262059/70‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:834369/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:3655453/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:3741473/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:1008504/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:819743/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:3960997/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:4087875/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:4140112/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:870237/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:4795856/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:4645289/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:4627196/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:942248/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:4268334/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:896165/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:449470/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:241045/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:1046426/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:1399144/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:1458364/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:1513779/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:1639355/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:1658330/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:1800171/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:2164462/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:2975992/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:3117565/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:3277669/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:31223/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCTGTCGACAGCTGATCGATCGGGAAACcgc                                      <  1:3115148/53‑1 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCaaaa                          >  1:1406714/1‑65 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCaaaa                          >  1:2605822/1‑65 (MQ=255)
                                             aCAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCaaaa                          >  1:533828/1‑65 (MQ=255)
                                              cAAACTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:4259178/52‑1 (MQ=255)
                                               aaaCTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:3587375/51‑1 (MQ=255)
                                               aaaCTTCCAGGTAGTTAGAATCTGTCGCCAGCTGATCGATCGGGAAACcgc                                      <  1:823338/51‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGTGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:2578684/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:3893899/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:3600517/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:1603910/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:3311383/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:3774604/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:145409/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:113859/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:330898/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:462972/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:2170333/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:2328706/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:3604096/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:437609/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:2794740/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:4263925/68‑1 (MQ=255)
                                                          tAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATc          <  1:4172283/68‑1 (MQ=255)
                                                                     ggTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATCAATAAAAg  >  1:1657266/1‑65 (MQ=255)
                                                                     |                                                                
GTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATCGGGAAACCGCGGTGCAGCAAAATACCTTCATCACCATCAATAAAAG  >  minE/460618‑460751

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: