Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,063,3040AC54.2% ‑3.1 / 12.6 27intergenic (+54/+90)aroP_1/CJ019_01001Aromatic amino acid transport protein AroP/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/10);  new base C (1/12);  total (2/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAG  >  ACB122/1063169‑1063432
                                                                                                                                        |                                                                                                                                
ctcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAg                                                                                                                                 <  6:369428/137‑1 (MQ=255)
    gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGcgcc                                                                                                                             <  3:343834/137‑1 (MQ=255)
        tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTAGGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCcgcgc                                                                                                                         <  5:128680/137‑1 (MQ=255)
                ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATAGAGCCCCGCccgggggggg                                                                                                                 >  6:371174/1‑135 (MQ=255)
                  ggCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCAGTTTGATAGCGCGCCCCCCGGGCGgggcg                                                                                                               >  2:53140/1‑136 (MQ=255)
                                ggggggggggggCTTGCACGCTGAACACGGCCGGCAGGCCACGGGCTTATCGGGTTTGTCCGTACGGGCTATGCCGTTAGGTTTCCGGCGCCTGGCCCACCCCGCGCCGCGCGGGCGGCGCTCTATCTCACAGGCGa                                                                                                <  4:82576/131‑1 (MQ=255)
                                        gtggCTTTCACGCTGACCACGGCGCGCCCGGCCAAGGGCGGAGCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAGCGCGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                         <  8:278937/137‑1 (MQ=37)
                                                       aagacgGCGCGCCAGGCCAAGGGCGGAGCGGGTTTGTTCGTACGGGGTAGGCCTTTAGGTTTTCGGCGCCTGTGAGAGCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                          <  2:8409/137‑1 (MQ=25)
                                                                  cccggCCAAGGGGGGAGCGGGTTTGTGCGAACGGGGAAGGCCTTTAGGTTTGCGGCGCCTGTCAGAGCGCGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                                               <  4:157102/134‑1 (MQ=25)
                                                                            ggCTGATCGGGTTTGTTCTAACGGGGTAGCCCGTTCGGTTTTCGGCGCCTTGCAGACCGCGCGCCGGGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTccc                                                     <  8:208326/137‑1 (MQ=37)
                                                                               tGATCGGGTTTGTGCGTACGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGCGCCCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                                  <  4:88468/137‑1 (MQ=255)
                                                                                 ttggggTTTGTTCGTAGGGGGTAGCCGTTTAGGTTTGCGGCGCGTGGGAGAGCCCGGCCGGCGGGGGGGGCGCGCCATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                                <  3:164874/134‑1 (MQ=255)
                                                                                 tcggggTTTGTGCGTACGGGGTATGCGGTTTGGTTTGCGGCGCCTGTCCGCTCCGGCGCGGCGCGGGCGGCGGTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                                <  1:231490/134‑1 (MQ=255)
                                                                                 cgcGGGTTTGTTCGTCAGGGGCAGCTCGTTTGGTTTGCGGCGCCTGTGCGAGCGCCCGCCGCGCGGGCGGCGCGCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                                <  4:35913/135‑1 (MQ=255)
                                                                                 agcGGGTTTGTTCGTACGGGGTATGCCGTAAGGTTTGCGGCGCCTGGGAGAGCCCGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                                <  1:53140/135‑1 (MQ=38)
                                                                                                          ccGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCa                       <  6:146745/137‑1 (MQ=255)
                                                                                                              tttggTTTTCGGCCCCGGTGAGAGGGTGCGCCGCGGGGGCGGGGGTCTATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgc                    <  5:85887/134‑1 (MQ=255)
                                                                                                              tttggTTTGCGCCGCCTGTGAGATCGCGCCCCGCGCGGGCGGCGCCCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATACCCCTCAATGCTGCCTCCGCCAAgc                    <  7:15369/134‑1 (MQ=255)
                                                                                                                tgtttttCGGCGCCTGTCTCAGCGCGCGCCGCGCGGGCGGCGAGCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac                  <  7:60610/134‑1 (MQ=255)
                                                                                                                  gTTTTGGGCGCCTGGGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc                <  1:215503/137‑1 (MQ=255)
                                                                                                                   ttttCGGCGCCTGGGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca               <  1:87303/137‑1 (MQ=255)
                                                                                                                   ttttCCGCGCCTTTGAGAGCGGGCGCGGCGCGGGCGGCGGTCTATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCATTGCTGCCTCCTCCAAgcacgca               <  5:80866/137‑1 (MQ=255)
                                                                                                                        ggCGCCTGGGAGCGCGACCGCCGCGCGGGCGGCGCGCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg          <  4:133783/137‑1 (MQ=255)
                                                                                                                        ggCGCCTGGGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg          <  5:371174/137‑1 (MQ=255)
                                                                                                                          cgcCGGTGACATCGAGCGCCGGGCGGGCGGCGCGCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggt        <  1:98987/137‑1 (MQ=255)
                                                                                                                             cggtgACATCCCGCGCCGCGCGGGCGGCGCGCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa     <  3:170133/135‑1 (MQ=255)
                                                                                                                                tgAGATCGCGCGCGGCGCGGGCGGCGCGCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAg  <  4:318942/137‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                
CTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAG  >  ACB122/1063169‑1063432

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