Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730512,6430AG69.6% 22.7 / 16.5 23intergenic (+711/‑150)recR/rrs1recombination protein RecR/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base G (12/4);  total (16/7)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGATGACAAAGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA  >  CP000730/512494‑512783
                                                                                                                                                     |                                                                                                                                            
aGATGACAAAGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATATTTTTTATTTGTTAAGTTTAACAAATTAATTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa                                                                                                                                             <  1:42269/151‑1 (MQ=255)
         cacagaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGta                                                                                                                                    <  2:252941‑M1/34‑1 (MQ=255)
           cagaagttaagctccttagcgtcgatggtagtcgaacttacgttgcgcttgagtagaacgttgccaggcagttttttagtcaaattgtgggtaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTAGATAACATGtata                                                                                                                                  <  2:452861‑M1/36‑1 (MQ=255)
                   aagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTc                                                                                                                          >  1:116530‑M1/108‑151 (MQ=255)
                          ttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCgg                                                                                                                   >  2:353046‑M1/101‑151 (MQ=255)
                            gttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGAAAACAAAACAACTCAACTTAAGAAATAACTTGACTTACTAAAGAACATTATTTTTAAAAAATTATTGAATTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTa                                                                                                                 >  1:294221/1‑151 (MQ=255)
                                     tttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATg                                                                                                        <  2:111289/151‑1 (MQ=255)
                                                 tacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAc                                                                                            >  2:158682‑M1/78‑151 (MQ=255)
                                                       ccgctagagtagaacgttgccaggcagtattttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTCTTGACGTAAATGTTGATAAAATGTATAATTAATTCTAGTCGGTAAGAAAAATGAACATTTAAAACTGAATg                                                                                      <  2:279499‑M1/80‑1 (MQ=255)
                                                                      caaacaaaACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                                                                       >  1:308347/1‑151 (MQ=255)
                                                                           caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAATTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                                                                  >  1:385071‑M1/52‑151 (MQ=255)
                                                                                      ttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt                                                       <  1:186666‑M1/111‑1 (MQ=255)
                                                                                            gAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTAAAAAACAAAAAACGTAACtat                                                 <  2:308347/151‑1 (MQ=255)
                                                                                             aaaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                                <  1:162470/151‑1 (MQ=255)
                                                                                                      taaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                       >  1:314355‑M1/25‑151 (MQ=255)
                                                                                                                      caagataaAAAAAGTTATTGGCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAATATTGAAAACTGAATGATAATATGTCAACGTTAATTACAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTa                       >  2:195245‑M1/9‑151 (MQ=255)
                                                                                                                         attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa                    >  2:102924/1‑151 (MQ=255)
                                                                                                                           taaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGc                  >  1:280564‑M1/4‑151 (MQ=255)
                                                                                                                                     tATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAca        >  2:106286/1‑151 (MQ=255)
                                                                                                                                     naTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAca        >  2:434143/2‑151 (MQ=255)
                                                                                                                                         gACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTCTGAGCTAATCAAAcatcat    >  2:224224/1‑151 (MQ=255)
                                                                                                                                         gACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat    >  2:332257/1‑151 (MQ=255)
                                                                                                                                          aCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATa   >  2:301664/1‑151 (MQ=255)
                                                                                                                                           cTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATaa  >  2:306256/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
AGATGACAAAGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA  >  CP000730/512494‑512783

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