Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 512,643 A→G 18.8% intergenic (+711/‑150) recR → / → rrs1 recombination protein RecR/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730512,6430AG18.8% 67.3 / 7.2 33intergenic (+711/‑150)recR/rrs1recombination protein RecR/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (11/15);  new base G (2/4);  total (13/20)
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.00e+00

GGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATC  >  CP000730/512504‑512779
                                                                                                                                           |                                                                                                                                        
ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta                                                                                                                                          >  1:113766/1‑140 (MQ=255)
      aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCACACAGAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATATTACAAAGTTATTGACTTAAATGTTAATaaaa                                                                                                                                    >  1:40944/1‑140 (MQ=255)
         tttACTGTTGTTAAGCAGTTTCAAAGACTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGt                                                                                                                                 <  2:134725/140‑1 (MQ=255)
           tACTGTTGTTAAGCAGTTTGAATGCCTGTATAGTATTTATTTGTTTAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtat                                                                                                                               <  2:4293/140‑1 (MQ=255)
                ttgttAAGCAGTTTGAAGGCCGTTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAACAACTCATCCCAGATGGCGGACTTAAATGTTAATAAAATGTAtaatta                                                                                                                          >  2:224400/1‑140 (MQ=255)
                   ttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATaattaatt                                                                                                                       <  2:283334/140‑1 (MQ=255)
                     aaGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCt                                                                                                                     <  1:197410/140‑1 (MQ=255)
                                taGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACTAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                                                                          <  2:253015/139‑1 (MQ=255)
                                 aGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCggcacgga                                                                                                         >  2:79213/1‑134 (MQ=255)
                                 aGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                                                                         >  1:85551/1‑140 (MQ=255)
                                          aGTAATTATTTGTTTAGGCAAACTAAACAATTCTACTTAAGATATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                                                                <  1:66431/140‑1 (MQ=255)
                                                  tttGTTGAGGCAAACACAACAACTCAACTTACGAAATAACTTGAATTACTAACGAAAATTAATCTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAc                                                                                        >  2:91868/1‑140 (MQ=255)
                                                    tgttgAGGTAAACAAAACAACTAAACTTAAGAAATAATTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTAACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg                                                                                      <  2:108133/140‑1 (MQ=255)
                                                    tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCNGTANGNAAAATGAACATTGAAAAATg                                                                                      >  1:226194/1‑140 (MQ=255)
                                                    tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg                                                                                      >  1:97959/1‑140 (MQ=255)
                                                      ttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTTGGTGAGAAAAAGGAACATTGAAAACTGaa                                                                                    >  1:72369/1‑140 (MQ=255)
                                                      ttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa                                                                                    >  2:205341/1‑140 (MQ=255)
                                                       tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAAt                                                                                   >  1:188427/1‑140 (MQ=255)
                                                              aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                                                                            <  2:20809/140‑1 (MQ=255)
                                                                caaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATg                                                                          <  2:168549/140‑1 (MQ=255)
                                                                              ttAAGACATAATTTGAATTACTAACGAAAATTAATTTTAAAAAGTAATTGACTAAAATGTTAATAAAATGTATAATGAATCCTTGTCGGTAAGAAAAAGGAACATTGAAAACTGAATGACAATATGTCAAAGTTAATTcc                                                            <  2:40944/140‑1 (MQ=38)
                                                                                aataaaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                          <  1:96131/137‑1 (MQ=255)
                                                                                 aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                         <  1:105302/140‑1 (MQ=255)
                                                                                                 aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATCGTTATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGt                                         <  1:298175/140‑1 (MQ=255)
                                                                                                     aCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAACCTGTATGATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                     <  1:231446/140‑1 (MQ=255)
                                                                                                     aCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTCTTTCTAAAAAACGTAACTATAAGTTACa                                     <  1:77041/140‑1 (MQ=255)
                                                                                                        aaaaaaaaataaaaaaaaTTATTGACATAAATGTTGATAAAATGTAAAAGAAATTCTTGTAGGTAAGAAAAATGAAAATTGAAAACTGAATGACAATATGTAAAAAAAAAAAAAAAAAAACGTAACTATAAGTTACAAAc                                  <  2:234744/128‑1 (MQ=16)
                                                                                                         ggacaagataaaaaaaGTTATTGACTTAAATGTTGATAAGATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACa                                 <  1:226786/129‑1 (MQ=25)
                                                                                                           acaagataaaaaaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatt                               <  2:13792/131‑1 (MQ=255)
                                                                                                                           tATTGACTTAAATGATGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAAGATGTCAACGTTACTTCCAAAAAACGTAACTATAAGTTAAAACAATTATTTAGTATTTAtacg               >  2:218439/1‑137 (MQ=21)
                                                                                                                            aTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTAATGAGc              >  1:38279/1‑140 (MQ=255)
                                                                                                                               gACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa           >  2:264570/1‑140 (MQ=255)
                                                                                                                                 cTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATc         <  2:557/140‑1 (MQ=255)
                                                                                                                                        gTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatc  <  2:97959/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
GGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATC  >  CP000730/512504‑512779

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