Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE277,4270GC43.5% 6.6 / 30.4 23A84A (GCG→GCCUSA300HOU_RS01215hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (9/4);  new base C (3/7);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 9.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCCGACATACATAATGAATCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACGCCAGCAATTTGTGACATATC  >  USA300TCH1516_ALE/277296‑277561
                                                                                                                                   |                                                                                                                                      
tCCGACATACATAATGAATCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAatgat                                                                                                                               >  2:242143/1‑141 (MQ=255)
tCCGACATACATAATGAATCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAatgat                                                                                                                               >  1:395218/1‑141 (MQ=255)
tCCGACATACATAATGAATCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAatgat                                                                                                                               >  1:287036/1‑141 (MQ=255)
tCCGACATACATAATGAATCTATCAATAAACAAATGGAAGCCAAAGCGCATGAAACAGCGCAAGATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAatgat                                                                                                                               >  1:103630/1‑141 (MQ=255)
tCCGACATACATAATGAATCTATCAATAAACAAATGGAAGCCAAAGCGCATGAAACAGCGCAAAATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAatgat                                                                                                                               >  1:318830/1‑141 (MQ=255)
         cataATGAATCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTa                                                                                                                      >  2:333137/1‑141 (MQ=255)
               gAATCTATCAATAAACAAATGGAAGCCAAAGCGCATGAAACAGCGCAAAATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAATGATGAAAGCTTAAATCtt                                                                                                                <  2:41229/141‑1 (MQ=255)
                 aTCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTaa                                                                                                              >  1:33555/1‑141 (MQ=255)
                 aTCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTaa                                                                                                              >  2:184759/1‑141 (MQ=255)
                             aCAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGatat                                                                                                  <  2:137276/141‑1 (MQ=255)
                                         caaaGCGCATGAAACAGCGCAAAATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGTTTTCTGAAGt                                                                                      <  2:144451/140‑1 (MQ=255)
                                         caaaGCGCATGAAACAGCGCAAAATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGTTTTCTGAAGt                                                                                      <  1:174003/140‑1 (MQ=255)
                                         caaaGCGCATGAAACAGCGCAAAATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGTTTTCTGAAGt                                                                                      <  1:173985/140‑1 (MQ=255)
                                             gcgcATGAAACAGCGCAAAATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGTTTTCTGAAGTAttt                                                                                  >  1:115885/1‑141 (MQ=255)
                                                                     gATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCa                                                          <  1:242143/141‑1 (MQ=255)
                                                                     gATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCa                                                          <  2:395218/141‑1 (MQ=255)
                                                                                                         gCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACg                      >  1:37373/1‑141 (MQ=255)
                                                                                                         gCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACg                      >  2:39228/1‑141 (MQ=255)
                                                                                                                  tcaatcGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACGCCAGCAAtt             >  1:346582/1‑141 (MQ=255)
                                                                                                                      tcGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACGCCAGCAATTtgtg         <  1:290786/141‑1 (MQ=255)
                                                                                                                             aCTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGGTATGTTTTCTGAAGTATTTAAGCCACATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACGCCAGCAATTTGTGATATATc  <  2:263930/141‑1 (MQ=255)
                                                                                                                             aCTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGTTTTCTGAAGTATTTAAGCCACATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACGCCAGCAATTTGTGATATATc  <  2:267730/141‑1 (MQ=255)
                                                                                                                             aCTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGTTTTCTGAAGTATTTAAGCCACATACTAAAAACGAAGCAGATGAAATATTTATAGAAAAAAAAAAAAAAACTACGCCAGCAATTTGTGATATATc  <  2:1399/141‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                      
TCCGACATACATAATGAATCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACGCCAGCAATTTGTGACATATC  >  USA300TCH1516_ALE/277296‑277561

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