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.

CTGGCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAATGGGGTTGGACGGGACCGTGTAAACCGCCGTGGGCTGGGTGATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGCGAACAGCCC  >  ACB143/3608364‑3608621
                                                                                                                                          |                                                                                                                           
cTGGCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAATGGGGTTGGACGGGACCGTGTAAACCGCCGTGGGCTGGGTGAtt                                                                                                                             <  8:91743/135‑1 (MQ=255)
   gCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCGGGAGGCGGATTTCTGGCAAATGGGGGTGGACGGGACCGTGTAAACCGCCGTGGGCTGGGTGATtact                                                                                                                          <  2:277590/135‑1 (MQ=255)
       cgcgcAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAATGGGGTTGGACGGGACCGTGTAAACCGCCGTGAGCTGGGTGATTACTAccca                                                                                                                     <  5:230207/136‑1 (MQ=255)
                             gTTCGCCTAAGTGGCTGGGTTTTTGGGCCCGGCCTCCCGCGGGGAGGGGGATTGCTGGAAAAGGGGGTGGGCCGGCACCGTGTAACCCGCCGGGGGCGGGGTGATAACTACCCCGCCCAGGACCGCTTTGGTGCAt                                                                                               <  6:94782/136‑1 (MQ=255)
                                 gCCGAACTGGCAGGTATTTTTGGCCAGGCCGCCCGCCGGGAGGCGGATTTCTGGCAAATGGGGTTGGCCGGGACCGTGTAACCCGCCGGGGGCTGGGTGATAACTACCCAGCCCAGGACCGCTTTGGTGCATggg                                                                                            <  2:344562/135‑1 (MQ=255)
                                                                     cgggAGGCGTATTTCTGGCAATTGGGGTGGCACGGCACCGTGAAAACCGCCGTGGGCGGGGTGATTACTCCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCg                                                        <  7:87222/133‑1 (MQ=255)
                                                                          tgctgATTGCTGGCAATTGGGGTTGGAAAGGACGGTGTAAACCGCCGGGGGCTGGGTGATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGa                                                   <  7:290790/131‑1 (MQ=255)
                                                                          ggCGGATTTGTGCAAATGGGGGGGGGCCGGCCCCGTTAAACCCGCGGGGGGGTGTTTGATTACTACCCCCCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGa                                                   <  5:241073/135‑1 (MQ=255)
                                                                               ttttCTGTCAATTGGGGTTGGAGGGGACCGTGAAACCCGCGGGGGGCGGGGTTATTACAACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCa                                              <  8:234485/134‑1 (MQ=255)
                                                                               aTTTCTGCAAATTGGGGGTGGCCGGCCCGGTGAAACCCGCCGGGGGGGGGGTGATTACACCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCa                                              <  7:43098/135‑1 (MQ=255)
                                                                                  tCTGCCAATTGGGGTTGGCCGGGACCGTGTAACCCGCCGTGGGCGGGGTGATAACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCt                                           <  6:191788/135‑1 (MQ=255)
                                                                                        aattGGGGTTGGACGGGACGGTGTAAACCGCGGGGGGCGGGTTTATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                                    <  5:147098/133‑1 (MQ=255)
                                                                                         attGGGGTTGGCGGGGACCTTTAAACCCGCCGTGGGTTGGGTTATAAAAACCCCGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                                    <  7:18483/133‑1 (MQ=255)
                                                                                         attGGGGTTGGACGGGACGGTGAAAACCGCCGTGGGCTGGGTTATAACAACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                                    <  2:93414/133‑1 (MQ=255)
                                                                                         attGGGGTGGCCCGGGACCGTGTAAACCGCCGTGGGCTGGTTTATTACAACCCCGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                                    <  3:291247/133‑1 (MQ=255)
                                                                                         aatgggggTGGGCCGGGACGGTGAAACCCGGGGTGGGCTGGGTTATAACAACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                                    <  5:109511/132‑1 (MQ=255)
                                                                                         aaTGGGGGTGGCCCGGGACGGTTAAAACCGCCGTGGGCTGGTTTATAACAACCCCCCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                                    <  5:166323/136‑1 (MQ=255)
                                                                                           tGGGGTGGCACGGCACGGTGTAAACCGCGGTGGGCTGGGTTATAACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGc                                  <  1:297608/135‑1 (MQ=255)
                                                                                            ttgggTTGGGCGGGTCCGGTTAAACCGGCGGTGGGGTGGGTTATAACAACCCCGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgtg                                   <  3:406705/133‑1 (MQ=255)
                                                                                             ggggTTGGGAGGGCCCGGTTTAACCCGCCGTGGGCGGGGTTATTACACCCCCGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGc                                  <  2:355262/135‑1 (MQ=255)
                                                                                                          gACGGTGTAAACCGCCGTGGGGTGGTTGATTACTCCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATgg                     <  8:193856/135‑1 (MQ=255)
                                                                                                                     cccgccgGGGGCTGGGTGATAACACCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGc           <  2:75247/134‑1 (MQ=255)
                                                                                                                     cccgccgGGGGCGGGGTGATTACACCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGc           <  2:384705/134‑1 (MQ=255)
                                                                                                                         gccgTGGGCTGGGTGATTACTACCCACCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGTACAGTGCATATGGATTTCTGGGCGaa        >  5:269878/1‑135 (MQ=255)
                                                                                                                               ggCGGGGTGAAAACACCCCCGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGCGAACAGccc  <  8:208465/135‑1 (MQ=255)
                                                                                                                                          |                                                                                                                           
CTGGCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAATGGGGTTGGACGGGACCGTGTAAACCGCCGTGGGCTGGGTGATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGCGAACAGCCC  >  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: