Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,070,2461.G76.9% 69.4 / 16.7 27intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger related‑protein
Reads supporting (aligned to +/- strand):  ref base . (0/6);  new base G (13/7);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 1.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACTGGGCGAGCATCAGCATTACCCACCTCAAATTCCGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGT  >  AE015451/1069952‑1070490
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      
aactggtcaaccagcAGCATTACCCACAGCAGATTCCGCAAGGGCAAGCTGGCCGCCGGGATCACCCCCTTCTACAGAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGGCGCCGCCGGTTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑Gagatc                                                                                                                                                                                                                                                 <  1:499031/287‑5 (MQ=255)
 aCTGGGCGAGCATCAGCATTACCCACCTCAAATTCCGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑Gagatcg                                                                                                                                                                                                                                                <  2:681889/301‑6 (MQ=255)
   tGGGCGAGCATCAGCATTACCCACCTCAAATTCCGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑Gagatcgag                                                                                                                                                                                                                                              <  1:66938/301‑8 (MQ=255)
                                ataccgtcaGGCCAAGCTGGCCGCCGTGATCACCCCCTTCTACACAAGCCTGGGGCACTCGCGTACCAACTACCTGTGGCTGACGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                                                                                                                                                                   <  1:571284/293‑1 (MQ=255)
                                                            aTCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCAGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCCCGGTCGGCGCCCGATCCCCCAGCCGATGAAAGCCTAAAGTCAtccgcctg                                                                                                                                                                                       >  1:348478/1‑294 (MQ=255)
                                                                           aaaaGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTACAGGCATGCGTCTGGCGCGCCTGACGCGt                                                                                                                                                                        >  1:134319/1‑301 (MQ=255)
                                                                                  ggcggCACCCACTGACCAACTACCGGTGCCTGGCGTTCATCGTGCTGGTGCTGGTAGGGATGTACATGACGACGCCGATTCGCATCGCGGTAATGCGGATCCCGGAATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                                                                                                                                                 <  2:409506/298‑1 (MQ=255)
                                                                                                 ccAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTTAGATCGAGCGCCGCGCGGGAGGCGCCCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACTCGTTCCCCTTCCCTGccccctcatt                                                                                                                                                  >  2:146644/1‑295 (MQ=255)
                                                                                                   aaCTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCACATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                                                                                                                                >  1:552090/1‑301 (MQ=255)
                                                                                                              ccTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCACATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                                                                                                     <  2:551904/301‑1 (MQ=255)
                                                                                                                tGGCGTTCAGCGTGCGTATCCTGGTAGGTATGTACCTGACGCCGCCGATTCGCATCTCGGTAAGTCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgc                                                                                                                                   <  2:137080/301‑1 (MQ=255)
                                                                                                                                      ggTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAac                                                                                                            >  1:819814/1‑301 (MQ=255)
                                                                                                                                                      gtgcgtccgATTCGAATGTCGGAAATGTTGTGCCCGGCATGGATGGATGTGCTGTGGGTGACTTTCAAGCTGAAGAAGGCTAGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCg                                                                                            <  1:752255/295‑1 (MQ=255)
                                                                                                                                                             cgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCCCGCACAGTGGTTTACAGCAACACCGTACGGAGCCCCGGAAGTAc                                                                                     >  2:729097/1‑301 (MQ=255)
                                                                                                                                                             cgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCCCCGCCACGCACGCACAGTGGTGTCACGCAACACCGTAGCCAGCAcccgccgcac                                                                                     >  2:630403/1‑293 (MQ=255)
                                                                                                                                                                 tCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCCCAATTCTGCCTCCGCCAAGCACGCACAGTGGGGTACAGCAACACCGTAGCGAGCACCGGCAGTACCTAg                                                                                 >  1:790203/1‑301 (MQ=255)
                                                                                                                                                                          ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATc                                                                        >  1:617124/1‑301 (MQ=255)
                                                                                                                                                                             aaTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGTC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGGTCCCCTTCCATGCCCCTCAATGCTGCTCCCGCCAAGCACGCACCGTGGGGTACAGCAACACCGTAGCGAGCACCGGCAGTATACAGAACAGCATCaaa                                                                     >  1:827483/1‑300 (MQ=255)
                                                                                                                                                                                  ggATCCCGGCATGGATGGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTc                                                                <  2:819524/300‑1 (MQ=255)
                                                                                                                                                                                       ctggatgggatggcgtgtgtgtgGGTTGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGGTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAggc                                                           <  2:616911/282‑1 (MQ=255)
                                                                                                                                                                                                     tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACgcc                                             >  1:767831/1‑299 (MQ=255)
                                                                                                                                                                                                           gtgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCGGT‑‑GA‑‑GATCGAGCGCCGCGCGGGAGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATAccc                                       <  1:146695/301‑1 (MQ=255)
                                                                                                                                                                                                               gtggggttCAAGCCGAAGAAGGCTCGCAAGGCAAAGGGCTGATCGGGTTTGTTCGGAAGGGGGATGCCGTTAGGTTTCCGGCGCC‑‑TGGGA‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATgg                                   <  1:567668/295‑1 (MQ=255)
                                                                                                                                                                                                               gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATgg                                   >  1:521295/1‑301 (MQ=255)
                                                                                                                                                                                                                gggttttCAAGCTGAAGGAGGCTCGCCAGGCCAAGGGGTGATCGGGTTTGTTCGTAGGGGGTATGTTCTAAGGGTTTCTGAGCC‑‑TGCGA‑‑GATCGCGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGACTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGt                                  <  2:437381/297‑1 (MQ=255)
                                                                                                                                                                                                                                            ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGAACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGACACACAAAGGCCGTAGGcc      >  1:171698/1‑301 (MQ=255)
                                                                                                                                                                                                                                                aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCGCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCCTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCATGCGGCCAGCCAACGGCATACCCATGGTAACAGCCGCCACATGAAGGCCGTAgtacagga  >  1:671768/1‑294 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      
AACTGGGCGAGCATCAGCATTACCCACCTCAAATTCCGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGT  >  AE015451/1069952‑1070490

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: