Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 512,643 A→G 51.0% 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,6430AG51.0% ‑4.9 / 58.9 45intergenic (+711/‑150)recR/rrs1recombination protein RecR/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (12/10);  new base G (13/10);  total (25/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 = 5.94e-01

TTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAT  >  CP000730/512511‑512771
                                                                                                                                    |                                                                                                                                
ttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATg                                                                                                                          >  1:217264‑M1/110‑141 (MQ=255)
ttaagctccttagcgtcgatgctagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATg                                                                                                                          >  2:169156‑M1/110‑141 (MQ=255)
ttaagcgccttagcgtcgatggtagtcgaacttacgttcccctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATg                                                                                                                          >  2:160644‑M1/110‑141 (MQ=255)
    tACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGCATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGCATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtata                                                                                                                      >  2:189055/1‑141 (MQ=255)
       tgttgtTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatt                                                                                                                   <  1:128917/141‑1 (MQ=255)
         ttagcgtcgatggtagtcgaacttacgttccgatagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaa                                                                                                                 <  1:117323‑M1/41‑1 (MQ=255)
              gtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCtt                                                                                                            <  2:269385‑M1/46‑1 (MQ=255)
                        aaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGTTAAGa                                                                                                  >  1:160534/1‑141 (MQ=255)
                         tcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                                                                 <  2:66917‑M1/57‑1 (MQ=255)
                          cgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaa                                                                                                >  1:66917‑M1/84‑141 (MQ=255)
                           ccactgccgttccccaagcgtagaacgtgttcaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa                                                                                               >  1:64733‑M1/83‑141 (MQ=255)
                            ccTGTATAGTATTTTTTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaaaa                                                                                              <  1:177735/141‑1 (MQ=255)
                                tacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGaa                                                                                          >  1:188579‑M1/78‑141 (MQ=255)
                                       cgctagagtagaacgttgccgggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCATGTCGGTAAGCAAAATGAACATTGaa                                                                                   <  1:173966‑M1/71‑1 (MQ=255)
                                         taaaagaaaaacattgtcaggcagttttttaaacaaattttggtaaaaaaataaaaagaacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                 >  2:265203‑M1/69‑141 (MQ=255)
                                          atttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAGAATTAATTGTTGTCGGTAAGAAAAATGAACATTGAAAAc                                                                                >  1:118069/1‑141 (MQ=255)
                                                      aaacaaaacaACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAACATGAACATTGAAAACTGCATGACActa                                                                    >  1:12050/1‑141 (MQ=255)
                                                                  tttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTCACTTAAATGTTGATAAAATGTATAGTTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTaa                                                        <  1:185004‑M1/98‑1 (MQ=255)
                                                                   ttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                                                       <  1:13215‑M1/99‑1 (MQ=255)
                                                                   caattTAAAAAATAACTTGAAATACTAAAGAAAATTAATTATAAATAGTTAAAAAAAAAAAAAATAATAAAATGTATAAGTAACTCTTGTCGGTAAGAAAAATGAACATCGAAAACTGAATGTTTATATGTCaaaaaaaaa                                                       <  2:106966/141‑8 (MQ=255)
                                                                      atataaaatttgtgttaaaaaaaaaagtggacaagataaAAAAAGTTAGTGAATTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc                                                    >  2:252833‑M1/40‑141 (MQ=255)
                                                                        atcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                  <  1:222652‑M1/104‑1 (MQ=255)
                                                                        aaaGAAATAATTTGAATTACTAATGAAAATTAATTTTAAAAAGGTATTGACTTAAATGTTAATAAAATGTATAATGAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTTAATGACAATATGTCAACGTTAATTCCaa                                                  <  1:3113/140‑1 (MQ=38)
                                                                          aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATGTGTCAACGTTAATTCCaaaa                                                >  2:62572/1‑141 (MQ=255)
                                                                          aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                >  2:178754/1‑141 (MQ=255)
                                                                          aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                >  2:283105/1‑141 (MQ=255)
                                                                              ttttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAaca                                            <  2:107686‑M1/110‑2 (MQ=255)
                                                                              tgttgtttaaaaaataaaatggacaagataaAAAAAGTTGTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg                                            >  1:231022‑M1/32‑141 (MQ=255)
                                                                               tttggttaaaaaataaaaaaaaaataaaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt                                           >  2:179129‑M1/31‑141 (MQ=255)
                                                                                 tggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCACAGTTAATTCCAAAAAACGtac                                         <  2:13362‑M1/113‑2 (MQ=255)
                                                                                     gAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtaa                                     >  2:130843/1‑140 (MQ=255)
                                                                                        aaaataaaatggaaaagaaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAAAAAAAAAAAAAAAACGTAACTATAAg                                  >  2:145796‑M1/22‑141 (MQ=255)
                                                                                         aaagaaaatagaaaagaaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGt                                 >  2:234127‑M1/21‑141 (MQ=255)
                                                                                          acaaaCGAAAATTAATTTTAAAAAGATATTGACTTAAACGTTAATAAAATGTATAATTAATTCTTATCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATACCAAAAAACGTAACTATAAGtt                                >  2:25977/4‑141 (MQ=37)
                                                                                          aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAGAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                <  2:87501/141‑1 (MQ=255)
                                                                                              aaatggacaagataaAAAAAGTTATTGACTTAAATGTGGATAAAATGTATAATTAATTCGTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa                            >  1:22766‑M1/16‑141 (MQ=255)
                                                                                               cGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                           <  1:75662/141‑1 (MQ=255)
                                                                                                    attaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                      >  1:111839/1‑141 (MQ=255)
                                                                                                        attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa                  >  2:200359/1‑141 (MQ=255)
                                                                                                         ttttAAAAAGTTATTGACTTAAATGTTAATAAAATGTAGAATTAATTCTTGTCTGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                 <  1:122455/141‑1 (MQ=255)
                                                                                                               aaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTa           <  1:261219/141‑1 (MQ=255)
                                                                                                                  gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa        >  1:10573/1‑141 (MQ=255)
                                                                                                                   ttttttaCTTCAAAGTAAATAAACTGTATAATTAATTCTTGTCGATAAGAAAAATGATCATTGAAAACTGAATGACAATATGTCAACGTTAATAACAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg       <  2:146572/135‑1 (MQ=17)
                                                                                                                   ttATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg       <  1:247051/141‑1 (MQ=255)
                                                                                                                        gACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt  <  1:107679/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                
TTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAT  >  CP000730/512511‑512771

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: