Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,608,4980TA47.8% ‑6.3 / 10.7 23intergenic (+21/‑208)tenA/codBAminopyrimidine aminohydrolase/Cytosine permease
Reads supporting (aligned to +/- strand):  ref base T (1/11);  new base A (0/11);  total (1/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAAT‑GGGGTTGG‑ACGGGACC‑GTGTAAACC‑GCCGTGGGCTGGGTGATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGCGAACAGCCC  >  ACB143/3608364‑3608621
                                                                                                                                          |                                                                                                                           
cTGGCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAAT‑GGGGTTGG‑ACGGGACC‑GTGTAAACC‑GCCGTGGGCTGGGTGAtt                                                                                                                             <  8:91743/135‑1 (MQ=255)
   gCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCGGGAGGCGGATTTCTGGCAAAT‑GGGGGTGG‑ACGGGACC‑GTGTAAACC‑GCCGTGGGCTGGGTGATtact                                                                                                                          <  2:277590/135‑1 (MQ=255)
       cgcgcAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAAT‑GGGGTTGG‑ACGGGACC‑GTGTAAACC‑GCCGTGAGCTGGGTGATTACTAccca                                                                                                                     <  5:230207/136‑1 (MQ=255)
                             gTTCGCCTAAGTGGCTGGGTTTTTGGGCCCGGCCTCCCGCGGGGAGGGGGATTGCTGGAAAAG‑GGGGTGGG‑CCGGCACC‑GTGTAACCC‑GCCGGGGGCGGGGTGATAACTACCCCGCCCAGGACCGCTTTGGTGCAt                                                                                               <  6:94782/136‑1 (MQ=255)
                                 gCCGAACTGGCAGGTATTTTTGGCCAGGCCGCCCGCCGGGAGGCGGATTTCTGGCAAAT‑GGGGTTGG‑CCGGGACC‑GTGTAACCC‑GCCGGGGGCTGGGTGATAACTACCCAGCCCAGGACCGCTTTGGTGCATggg                                                                                            <  2:344562/135‑1 (MQ=255)
                                                                     cgggAGGCGTATTTCTGGCAATT‑GGGGTGGC‑ACGGCACC‑GTGAAAACC‑GCCGTGGGCGGGGTGATTACTCCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCg                                                        <  7:87222/133‑1 (MQ=255)
                                                                          tgctgATTGCTGGCAATT‑GGGGTTGG‑AAAGGACG‑GTGTAAACC‑GCCGGGGGCTGGGTGATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGa                                                   <  7:290790/131‑1 (MQ=255)
                                                                          ggCGGATTTGT‑GCAAATGGGGGGGGG‑CCGGCCCC‑GTTAAACCC‑GCGGGGGGGTGTTTGATTACTACCCCCCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGa                                                   <  5:241073/135‑1 (MQ=255)
                                                                               ttttCTGTCAATT‑GGGGTTGG‑AGGGGACC‑GTGAAACCC‑GCGGGGGGCGGGGTTATTACAACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCa                                              <  8:234485/134‑1 (MQ=255)
                                                                               aTTTCTGCAAATT‑GGGGGTGG‑CCGGCCCG‑GTGAAACCC‑GCCGGGGGGGGGGTGATTACACCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCa                                              <  7:43098/135‑1 (MQ=255)
                                                                                  tCTGCCAATT‑GGGGTTGG‑CCGGGACC‑GTGTAACCC‑GCCGTGGGCGGGGTGATAACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCt                                           <  6:191788/135‑1 (MQ=255)
                                                                                        aatt‑GGGGTTGG‑ACGGGACG‑GTGTAAACC‑GCGGGGGGCGGGTTTATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                                    <  5:147098/133‑1 (MQ=255)
                                                                                         att‑GGGGTTGG‑CGGGGACC‑TTTAAACCC‑GCCGTGGGTTGGGTTATAAAAACCCCGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                                    <  7:18483/133‑1 (MQ=255)
                                                                                         att‑GGGGTTGG‑ACGGGACG‑GTGAAAACC‑GCCGTGGGCTGGGTTATAACAACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                                    <  2:93414/133‑1 (MQ=255)
                                                                                         att‑GGGGTGGC‑CCGGGACC‑GTGTAAACC‑GCCGTGGGCTGGTTTATTACAACCCCGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                                    <  3:291247/133‑1 (MQ=255)
                                                                                         aatgggggTGGG‑CCGGGACG‑GTGAAACCC‑GGGGTGGGCTGGGTTATAACAACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                                    <  5:109511/132‑1 (MQ=255)
                                                                                         aaT‑GGGGGTGGCCCGGGACG‑GTTAAAACC‑GCCGTGGGCTGGTTTATAACAACCCCCCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                                    <  5:166323/136‑1 (MQ=255)
                                                                                           t‑GGGGTGGC‑ACGGCACG‑GTGTAAACC‑GCGGTGGGCTGGGTTATAACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGc                                  <  1:297608/135‑1 (MQ=255)
                                                                                            ttgggTTGG‑GCGGGTCC‑GGTTAAACCGGCGGTGGGGTGGGTTATAACAACCCCGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgtg                                   <  3:406705/133‑1 (MQ=255)
                                                                                             ggggTTGG‑GAGGGCCCGGTTTAACCC‑GCCGTGGGCGGGGTTATTACACCCCCGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGc                                  <  2:355262/135‑1 (MQ=255)
                                                                                                          gACG‑GTGTAAACC‑GCCGTGGGGTGGTTGATTACTCCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATgg                     <  8:193856/135‑1 (MQ=255)
                                                                                                                     ccc‑gccgGGGGCTGGGTGATAACACCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGc           <  2:75247/134‑1 (MQ=255)
                                                                                                                     ccc‑gccgGGGGCGGGGTGATTACACCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGc           <  2:384705/134‑1 (MQ=255)
                                                                                                                         gccgTGGGCTGGGTGATTACTACCCACCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGTACAGTGCATATGGATTTCTGGGCGaa        >  5:269878/1‑135 (MQ=255)
                                                                                                                               ggCGGGGTGAAAACACCCCCGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGCGAACAGccc  <  8:208465/135‑1 (MQ=255)
                                                                                                                                          |                                                                                                                           
CTGGCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAAT‑GGGGTTGG‑ACGGGACC‑GTGTAAACC‑GCCGTGGGCTGGGTGATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGCGAACAGCCC  >  ACB143/3608364‑3608621

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