Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,808,1750C.45.5% 15.9 / 16.4 11intergenic (‑284/‑271)PP_4234/dsbD‑IIAldehyde oxidase and xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base . (3/2);  total (6/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGACATGCTTGATCCTCTGGGCAGTTGGGGCCTGCCCTTTCTTTTGATTCTTGGCCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCT  >  AE015451/4807886‑4808468
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                     
ggACATGCTTGATCCTCTGGGCAGTTGGGGCCTGCCCTTTCTTTTGATTCTTGGCCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgcggcgct                                                                                                                                                                                                                                                                                            <  1:828162/301‑8 (MQ=255)
                    gCAGTTGGGGCCTGCCCTTTCTTTTGATTCTTGGCCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgcggcgcttcgcgggtaaccccg                                                                                                                                                                                                                                                                             >  1:397310/1‑274 (MQ=255)
                                   ccTTTCTTTTGATTCTTGGCCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgcggcgctcgatct                                                                                                                                                                                                                                                                                      <  1:224382/272‑14 (MQ=255)
                                   ccTTTCTTTTGATTCTTGGCCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgcggcgctcgatct                                                                                                                                                                                                                                                                                      >  2:224377/1‑259 (MQ=255)
                                                                                                                 aaaGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGccgccccct                                                                                                                                                                                                                                                                                                   >  2:11786/1‑178 (MQ=255)
                                                                                                                 aaaGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGccgccccct                                                                                                                                                                                                                                                                                                   <  1:11788/181‑4 (MQ=255)
                                                                                                                   aGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTcc                                                                                                                                                                                                                     <  1:227982/243‑1 (MQ=255)
                                                                                                                   aGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTcc                                                                                                                                                                                                                     >  2:227977/1‑243 (MQ=255)
                                                                                                                                  gtgtCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCAcagctcgga                                                                                                                                                                                                           >  2:689625/1‑230 (MQ=255)
                                                                                                                                  gtgtCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCAc                                                                                                                                                                                                                   <  1:689650/230‑1 (MQ=255)
                                                                                                                                                                                                                                   aTTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTc                                                                                                                                                                                                                 <  1:920811/135‑1 (MQ=21)
                                                                                                                                                                                                                                   aTTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTc                                                                                                                                                                                                                 >  2:920785/1‑135 (MQ=21)
                                                                                                                                                                                                                                            gcCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTc                                                                                >  1:472587/1‑255 (MQ=255)
                                                                                                                                                                                                                                                                         aaacaaacgtaggagcgggtttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGa                             >  2:579856/24‑291 (MQ=255)
                                                                                                                                                                                                                                                                          aacaaacgtaggagcgggtttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGcgcg                  <  1:623885/279‑1 (MQ=255)
                                                                                                                                                                                                                                                                            caaacgtaggagcgggtttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGt                 >  1:705343/21‑300 (MQ=255)
                                                                                                                                                                                                                                                                                           gtttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCCCCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTTCCAAGGCGAAGACATGCGCGTTCGCCTGCTCTTCCt  >  1:434094/6‑300 (MQ=255)
                                                                                                                                                                                                                                                                                             ttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAGCGCCGCCTCCCCTCGAAAATCCAACCACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGaccc                                                                                                         <  2:743493/192‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                     
GGACATGCTTGATCCTCTGGGCAGTTGGGGCCTGCCCTTTCTTTTGATTCTTGGCCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCT  >  AE015451/4807886‑4808468

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