Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730110,6600TA25.0% 111.5 / 30.2 44I173I (ATT→ATA) USA300HOU_0109staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (12/21);  new base A (1/10);  total (13/31)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAATTTAAAAGG  >  CP000730/110526‑110806
                                                                                                                                      |                                                                                                                                                  
aaaaTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTTAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                                       <  1:289441/148‑1 (MQ=255)
aaaaTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                                       <  2:114773/148‑1 (MQ=255)
                  aTCATTCAAACTAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                                                                                     <  2:266277/148‑1 (MQ=255)
                   tCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAAACGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGCGGTGCCGAGTTATTCGGc                                                                                                                    >  2:280439/1‑148 (MQ=255)
                      ttCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGGGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                                                                                 <  2:376838/148‑1 (MQ=255)
                      ttCCAACAAAAAAAATTAAAGATAAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAAATTAAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                                                                                 <  2:208346/148‑1 (MQ=255)
                            caaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAGTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCa                                                                                                           <  1:263372/148‑1 (MQ=255)
                               aaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                        >  2:266504/1‑148 (MQ=255)
                                    aTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                                   >  1:102112/1‑148 (MQ=255)
                                    aTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                                   <  2:3693/148‑1 (MQ=255)
                                    aTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                                   <  2:278070/148‑1 (MQ=255)
                                    aTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                                   <  2:310813/148‑1 (MQ=255)
                                      tAAAAATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGAATTAAAAAAAAAAAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                                 <  2:310555/148‑1 (MQ=17)
                                        aaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                                               <  1:149204/148‑1 (MQ=255)
                                         aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                              <  2:103559/148‑1 (MQ=255)
                                         aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                              <  1:56254/148‑1 (MQ=255)
                                              aaaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                                                                         <  2:257812/148‑1 (MQ=255)
                                                       aaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                >  1:303806/1‑148 (MQ=255)
                                                       aaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAACACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                >  1:320074/1‑148 (MQ=255)
                                                                                       cAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacga                                                <  2:102112/148‑1 (MQ=255)
                                                                                           atGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                                            <  1:136090/148‑1 (MQ=255)
                                                                                               caattttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCt                                        <  2:215294/144‑1 (MQ=11)
                                                                                                   tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaa                                    >  1:198059/1‑148 (MQ=255)
                                                                                                    ttaaaggtttaatgaaTTAGAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaaa                                   <  2:61471/134‑1 (MQ=14)
                                                                                                    ttAAAAATTTGAAAAATTATAAAGACGGATATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAGTTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaaa                                   <  1:379046/148‑1 (MQ=255)
                                                                                                    ttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATCTAAGCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaaa                                   <  1:196688/148‑1 (MQ=255)
                                                                                                    ttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaaa                                   <  2:33101/148‑1 (MQ=255)
                                                                                                     taaagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAg                                  <  1:260361/135‑1 (MQ=14)
                                                                                                     taaagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAg                                  <  2:149266/135‑1 (MQ=11)
                                                                                                     taaagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAg                                  <  2:254066/135‑1 (MQ=14)
                                                                                                     taaagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAg                                  <  2:255877/135‑1 (MQ=14)
                                                                                                     tAAATATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAg                                  <  1:246927/148‑1 (MQ=255)
                                                                                                     tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAg                                  <  1:115774/148‑1 (MQ=255)
                                                                                                        aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa                               >  2:355713/1‑148 (MQ=255)
                                                                                                        aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa                               >  1:8408/1‑148 (MQ=255)
                                                                                                           tttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTg                            >  1:25701/1‑148 (MQ=255)
                                                                                                            ttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt                           >  2:362534/7‑148 (MQ=255)
                                                                                                                tgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaa                       <  1:251757/146‑1 (MQ=255)
                                                                                                                 aaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaaa                      <  2:13433/148‑1 (MQ=255)
                                                                                                                            aCGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGaa           >  2:291400/1‑148 (MQ=255)
                                                                                                                            aCGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGaa           >  2:33887/1‑148 (MQ=255)
                                                                                                                            aCGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGaa           >  1:233630/1‑148 (MQ=255)
                                                                                                                                agaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAATTTa       <  2:219503/148‑1 (MQ=255)
                                                                                                                                agaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAAAAAAAAAAAAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAATTTa       <  2:135762/148‑1 (MQ=18)
                                                                                                                                     tttCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAATTTAAAAgg  <  1:68853/148‑1 (MQ=255)
                                                                                                                                     tttCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAATTTAAAAgg  <  1:365063/148‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                  
AAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAATTTAAAAGG  >  CP000730/110526‑110806

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