Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE111,5160CA25.0% 61.7 / 15.6 33Q196K (CAA→AAA) USA300HOU_RS00520putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (13/11);  new base A (2/6);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAA  >  USA300TCH1516_ALE/111379‑111651
                                                                                                                                         |                                                                                                                                       
aGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                                                                                                      <  1:227814/141‑1 (MQ=255)
        aaaGAAGAAAACATAAAAACAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTa                                                                                                                              <  1:132943/141‑1 (MQ=255)
          agaagaAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAc                                                                                                                            <  1:212570/141‑1 (MQ=255)
          agaagaAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAATAAc                                                                                                                            <  2:26201/141‑1 (MQ=255)
          aaaagaaAACATAAAAAAATAAATAGAAAACTATAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAc                                                                                                                            <  1:161614/139‑1 (MQ=255)
                  aaaTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgatg                                                                                                                    <  1:58042/139‑1 (MQ=255)
                  aCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                                                                    >  1:125560/1‑141 (MQ=255)
                   cATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTTTTTTTTTTAGAAATTTTAAAGATTTAATGAATTATAAAAATGGAGATATATCATATAATCTAGAGTTGCCGAGTTATTCTGATAAATATCAATTAACTATtgatga                                                                                                                   >  2:124157/1‑141 (MQ=14)
                     tAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGAtt                                                                                                                 <  1:41326/141‑1 (MQ=255)
                     tAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGAtt                                                                                                                 <  1:103404/141‑1 (MQ=255)
                         aaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATaa                                                                                                             <  2:113165/141‑1 (MQ=255)
                            aGAAATCGAAAACTTTAAGTTTTTTGCGGAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGt                                                                                                          <  2:48061/141‑1 (MQ=255)
                              gaattcaaaaCTAGAAGTTATTTTCACAATATAGCTATTTTAAAGATTTAATTAATTATAAAGAGTGAGAAATATCATATAATCCAGAGGTGCCGAGTTATTCGCCTAAATATCAATTAATTAATGATGATTATAATGTaa                                                                                                        <  2:140841/135‑1 (MQ=14)
                                       aCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                                                                               >  1:133605/1‑141 (MQ=255)
                                            aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCTATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                                                          <  1:193031/141‑1 (MQ=17)
                                             aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaaaa                                                                                         <  1:215332/141‑1 (MQ=255)
                                                          atatAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAcc                                                                            <  1:184520/141‑1 (MQ=14)
                                                            ataGCAATTTTATAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAccga                                                                          <  1:366/141‑3 (MQ=14)
                                                                         aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTAATATGCCTAAAAATTTTTTCCTCTCGTCAACAACAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAACCGcc                                                             >  2:231960/1‑141 (MQ=12)
                                                                             ttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAg                                                         <  1:40498/141‑1 (MQ=255)
                                                                               aaTGAATTATAAAGATGGCGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATACTGTAAAACAATTAAGTAACAGATATGCCATACCCACACAGACACCGCCGCAGcc                                                       >  2:201889/1‑140 (MQ=17)
                                                                                aTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtg                                                      <  1:234532/141‑1 (MQ=255)
                                                                                         aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACTAACAAATAAAGCGCAGAAGATGTTTTTGaaa                                             >  1:43527/1‑141 (MQ=255)
                                                                                                   atatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGa                                   >  2:33151/1‑141 (MQ=255)
                                                                                                            ataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAg                          >  1:125632/1‑141 (MQ=255)
                                                                                                                 ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGAACAGGGAATTTAAAAGGTtca                     >  2:135196/1‑141 (MQ=255)
                                                                                                                 ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCCAAGCAGTTGGTGATAGGGACAGGGAATTTAAAAGGTtca                     >  2:178393/1‑141 (MQ=17)
                                                                                                                 ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATCACAAAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAACGGTtca                     >  2:85717/1‑141 (MQ=255)
                                                                                                                   agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGGATTTAAAAGGTtcatc                   >  1:153826/1‑141 (MQ=255)
                                                                                                                         gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTgg             >  2:149905/1‑141 (MQ=255)
                                                                                                                            gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTAAAGGGAATTTAAAAGGTGCATCAGTTGGata          >  2:121170/1‑141 (MQ=255)
                                                                                                                                   tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaaa   >  2:187655/1‑141 (MQ=255)
                                                                                                                                    cAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaaaa  >  1:179117/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                       
AGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAA  >  USA300TCH1516_ALE/111379‑111651

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