Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,313,1130AT33.3% 14.7 / 15.9 18intergenic (+17/+280)rny/USA300HOU_1219possible hydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (6/6);  new base T (2/4);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAAATCTTATCGATTGGCTAGAGATATTAAAAATCAGATTGAAGATGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGTAAGACGGACCTT  >  CP000730/1312979‑1313219
                                                                                                                                      |                                                                                                          
aaaaTCTTATCGATTGGCTAGAGATATTAAAAATCAGATTGAAGATGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAtttg                                                                                                      <  2:308161/141‑3 (MQ=255)
        aTCGATTGGCTAGAGATATTAAAAATCAGATTGAAGATGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAg                                                                                              >  1:150188/1‑141 (MQ=255)
         tCGATTGGCTAGAGATATTAAAAATCAGATTGAAGATGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGc                                                                                             <  2:318070/141‑1 (MQ=255)
                     agaTATTAAAAATCAGATTGAAGATGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAATTTGTGAGGGAgacaaaaattattt                                                                                 <  1:369293/141‑14 (MQ=38)
                                     aTTGAAGATGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGtt                                                                 >  2:331667/1‑141 (MQ=255)
                                             ggAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGGGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCa                                                         <  2:319006/140‑1 (MQ=255)
                                                         tCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATAGGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGtaa                                             >  1:29233/1‑141 (MQ=255)
                                                                cATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACtat                                      >  1:149813/1‑141 (MQ=255)
                                                                             cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAtt                                                                                                        >  1:240888/1‑62 (MQ=255)
                                                                             cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAtt                                                                                                        >  1:240960/1‑62 (MQ=255)
                                                                             cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAtt                                                                                                        <  2:240888/62‑1 (MQ=255)
                                                                             cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAtt                                                                                                        <  2:240960/62‑1 (MQ=255)
                                                                             cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtagg                         <  1:283627/141‑2 (MQ=255)
                                                                             cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCCCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGCTAATAACTATATTAAACagtagt                         <  1:96337/141‑1 (MQ=255)
                                                                                      gagagaCTAGAGCAGTAGAATATGTGAAATTATTTTTTTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGaaa                <  2:293927/141‑1 (MQ=255)
                                                                                                agCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGTAAGACGGa      <  1:10349/141‑1 (MQ=255)
                                                                                                  cAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGTAAGACGGAcc    >  1:206146/1‑141 (MQ=255)
                                                                                                    gTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGTAAGACGGACCtt  >  1:75716/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                          
AAAATCTTATCGATTGGCTAGAGATATTAAAAATCAGATTGAAGATGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGTAAGACGGACCTT  >  CP000730/1312979‑1313219

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