Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,290,7490AG22.2% 56.6 / 17.5 36intergenic (‑163/‑63)arg/USA300HOU_2155arginase/transposase
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base G (4/4);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCTATGTAATTAAAGAATAATGTATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTGCGCCTATGTGTAATGATACCTTAGAA  >  CP000730/2290603‑2290832
                                                                                                                                                  |                                                                                   
tctATGTAATTAAAGAATAATGTATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAgc                                                                                 <  2:138850/151‑1 (MQ=255)
  tATGTAATTAAAGAATAATGTATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAACAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTAAAAAAAAAAAAAAgcgc                                                                               <  2:847286/151‑1 (MQ=38)
      tAATTAAAGAATAATGTATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGCAGGTAAAACTTCCTGCTTTTTTTAATACACAAGAAGCGCAAtt                                                                           >  1:154826/1‑151 (MQ=255)
      tAATTAAAGAATAATGTATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAGAAGCGCAAtt                                                                           >  1:154650/1‑151 (MQ=255)
      tAATTAAAGAATAATGTATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAGAAGCGCAAtt                                                                           >  1:380136/1‑151 (MQ=255)
      tAATTAAAGAATAATGTATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAGAAGCGCAAtt                                                                           >  1:380177/1‑151 (MQ=255)
          tAAAGAATAATGTATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTAATTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCct                                                                       <  2:591536/151‑1 (MQ=255)
            aaGAATAATGTATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCctct                                                                     >  2:881911/1‑151 (MQ=255)
            aaGAATAATGTATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCctct                                                                     <  2:713083/151‑1 (MQ=255)
                      tATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTTTTTCATTAAAAGCGCAAATACCTCAATAATTTaaa                                                           >  2:763551/1‑151 (MQ=25)
                      tATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTaaa                                                           >  2:763275/1‑151 (MQ=255)
                       aTGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAg                                                          <  1:937525/151‑1 (MQ=255)
                                ccATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGAccaaac                                                 <  1:729161/151‑1 (MQ=255)
                                 cATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTAAAAAAAAAAAAAAAGCGCAATTGCCTCTATAATTTAAAGTGAccaaacc                                                <  2:522294/151‑1 (MQ=21)
                                   ttATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGAccaaaccca                                              <  2:455047/151‑1 (MQ=255)
                                     aTCAAGCAATAGCTACACATATCATCTGTTTATCTTATTACTTCATAATAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCTTGCTTTTTTTAATACACAAAGAGCGCAATTGCCTCTTTAATTTAAAGTGAccaaacccaaa                                            <  2:253746/151‑1 (MQ=255)
                                     aTCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGAccaaacccaaa                                            >  2:203450/1‑151 (MQ=255)
                                     aTCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCAACTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTACCTCTATAATTTAAAGTGAccaaacccaaa                                            >  2:506682/1‑151 (MQ=255)
                                      tCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGAccaaacccaaac                                           <  1:109977/151‑1 (MQ=255)
                                      tCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGAccaaacccaaac                                           <  2:599534/151‑1 (MQ=255)
                                        aaGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTa                                         >  2:229986/1‑151 (MQ=255)
                                        aaGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTa                                         >  1:399263/1‑151 (MQ=255)
                                           cAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAg                                      <  2:274361/151‑1 (MQ=255)
                                            aaTAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAgg                                     >  1:1037363/1‑151 (MQ=255)
                                                      caTATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAGAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTgc                           <  2:380177/151‑1 (MQ=255)
                                                      caTATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAGAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTgc                           <  2:380136/151‑1 (MQ=255)
                                                      caTATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAGAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTgc                           <  2:154826/151‑1 (MQ=255)
                                                      caTATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAGAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTgc                           <  2:154650/151‑1 (MQ=255)
                                                             tCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTGCGCCTAtg                    <  2:19636/151‑1 (MQ=255)
                                                             tCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTGCGCCTAtg                    <  2:737048/151‑1 (MQ=255)
                                                               tGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTGCGCCTAtgtg                  <  1:996079/151‑1 (MQ=255)
                                                               tGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTGCGCCTAtgtg                  >  1:739909/1‑151 (MQ=255)
                                                                    aTCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTGCGCCTATGTGTAATg             >  1:1015920/1‑151 (MQ=255)
                                                                             cTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTGCGCCTATGTGTAATGATACCTTAg    >  1:918562/1‑151 (MQ=38)
                                                                              ttCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTGCGCCTATGTGTAATGATACCTTAGa   >  2:244255/1‑151 (MQ=38)
                                                                               tCATAAAAAAAGGTTCTTCATCATTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAATGACTAGACCCAAAATAAAGGAGACAATTGGAACTATGTTTAATAATActtcagga  >  2:226270/1‑144 (MQ=11)
                                                                                                                                                  |                                                                                   
TCTATGTAATTAAAGAATAATGTATGCGCTTACCATTATCAAGCAATAGCTACACATATAATCTGTTTATCTTATTACTTCATAAAAAAAGGTTCTTCATCTTTTATGGTGGGAAGGTAAAACTTCCTGCTTTTTTTAATACACAAAAAGCGCAATTGCCTCTATAATTTAAAGTGACCAAACCCAAACTAAAGGAGACAAGTGCGCCTATGTGTAATGATACCTTAGAA  >  CP000730/2290603‑2290832

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