Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,443,7240T.20.0% 53.7 / 10.0 20intergenic (‑297/+56)USA300HOU_1341/brnQ2IS1272 transposase/LIVCS family branched chain amino acid:cation symporter
Reads supporting (aligned to +/- strand):  ref base T (6/10);  new base . (3/1);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGACCGAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCTAATATCGTTGCAATAA  >  CP000730/1443627‑1443853
                                                                                                 |                                                                                                                                 
aGACCGAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCtt                                                                                        <  1:193437/141‑1 (MQ=34)
 gACCGAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCttt                                                                                       <  2:71528/141‑1 (MQ=35)
            agaCGGCACACTATGAAAGCGAAGCCATTGAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCg                                                                            <  2:71210/141‑1 (MQ=255)
                 gCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTAtt                                                                       <  1:202812/141‑1 (MQ=255)
                             aGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATAtt                                                           >  2:265828/1‑141 (MQ=255)
                                     cATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGa                                                   <  2:364669/141‑1 (MQ=255)
                                       ttCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATc                                                 >  2:263876/1‑141 (MQ=255)
                                              cGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGtttt                                          >  1:122305/1‑141 (MQ=255)
                                               gAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTa                                         <  2:287240/141‑1 (MQ=255)
                                                aaGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTAc                                        <  1:106375/141‑1 (MQ=255)
                                                         ataAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGcc                               <  1:245548/141‑1 (MQ=255)
                                                            aaTAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGAc                            <  1:123196/141‑1 (MQ=255)
                                                                  gaACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATa                      >  1:29848/1‑141 (MQ=255)
                                                                        cagTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCTaa                >  1:156723/1‑141 (MQ=255)
                                                                           tAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCTAatat             <  1:103016/141‑1 (MQ=255)
                                                                                aattttcaattagAAAATATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCt                  <  2:192135/118‑1 (MQ=38)
                                                                                aattttcaattagAAAATATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCt                  >  1:192135/13‑130 (MQ=38)
                                                                                   tttcaattagAAAATATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCTAATATCGTTGCaat    >  2:312537/10‑141 (MQ=255)
                                                                                     tcaattagAAAATATCTTACTGCTGTTTTTTTTTTTTTTTTTTTTACAACCAGATTCGTATTATTCTTTATTCCGGTAGATGGTTAATTGGATCTTGTTTTACAATGATTGAGAGTATATAGCATAATATCGTAGCaatta  >  2:88178/8‑139 (MQ=255)
                                                                                      cTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGTTAGCCTAATATCGTTGCaataa  >  1:97165/1‑141 (MQ=255)
                                                                                                 |                                                                                                                                 
AGACCGAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATATTTTCTAATTGAAAATTATCTTACTGCTGTTTTTTTAGGGATTTATGTCCCAACCTCTTTTATATTATTCGTTATTCCTGTTGATATTTAATTGGATCTTGTTTTACAAATATGCCGACTAGATAGCCTAATATCGTTGCAATAA  >  CP000730/1443627‑1443853

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