Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 512,643 A→G 64.2% 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,6430AG64.2% 37.7 / 30.2 41intergenic (+711/‑150)recR/rrs1recombination protein RecR/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (10/5);  new base G (11/15);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCA  >  CP000730/512503‑512773
                                                                                                                                            |                                                                                                                                  
aGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta                                                                                                                                    <  1:137180/141‑1 (MQ=255)
   agaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTtgata                                                                                                                                 <  2:138410‑M1/27‑1 (MQ=255)
    ttattttctACTTTTGTCAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGACGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAACATTAATTTTAAAAAGTTATTGACTTAAATGTtaataa                                                                                                                                <  1:214592/137‑1 (MQ=255)
                gttgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatta                                                                                                                    <  2:234767/141‑1 (MQ=255)
                cttagcgtcgatggtaatcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatta                                                                                                                    <  1:169177‑M1/40‑1 (MQ=255)
                       tcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatccaattatggctaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTg                                                                                                             <  2:95026‑M1/47‑1 (MQ=255)
                               aaaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                     <  2:108277/140‑1 (MQ=255)
                                gtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATGAAATGTATAATTAATTATTGTCGGTAAGa                                                                                                    <  1:101551‑M1/56‑1 (MQ=255)
                                   gaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa                                                                                                 >  2:248663‑M1/83‑141 (MQ=255)
                                         acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc                                                                                           <  1:174690‑M1/65‑1 (MQ=255)
                                         acgttccgctagagtagaacgtgaccatacactgtttccatcaaattttggttaaaaaataaaatggacaagataaAAAAAATTAATGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATgaaa                                                                                           <  1:67911‑M1/65‑2 (MQ=255)
                                          cgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggataaagaaaaaaatgggcaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTATTGTAGGTAAGAAGAATGAACa                                                                                          <  2:156523‑M1/66‑1 (MQ=255)
                                             tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                       >  2:235689‑M1/73‑141 (MQ=255)
                                                     tgttgAAGCAAACAAAACAACTAAAATTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCAAGTAAGTAAAAAAAATGAACATTGAAAACTGa                                                                               <  2:251395/141‑1 (MQ=255)
                                                      gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa                                                                              >  2:27283/1‑141 (MQ=255)
                                                      gttgAGGAACAGAAAGCAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGac                                                                              >  2:275423/1‑140 (MQ=255)
                                                       ttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAAt                                                                             >  1:62452/1‑141 (MQ=255)
                                                        tcacgCAAACAAAACAACTCAACTTAAGTAATCACTCGAATTACTCACAAAGATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAACATATATAATTAATTTTTGTCTGTAAGAAAGAGGCACATTAaaaaatcgatt                                                                            >  2:123998/5‑134 (MQ=255)
                                                         gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAACCTGAATGa                                                                           <  1:167270‑M1/81‑1 (MQ=255)
                                                             gttgccaggcagttttttaatcaaattttgggtaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGGCGGTAAGAAAAAGGAACATTGaaaactgaagtgcaat                                                                       <  1:142383‑M1/85‑8 (MQ=255)
                                                                     gtagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                                                               >  2:112949‑M1/49‑141 (MQ=255)
                                                                     acaacTCGACTTAAGAAATAACTTGAATTACTAACGCAAATTAAATTTAAAAAGTTATTGACTTACAGGTTAATAACATGTATAATTAATTCTCGTCGGTAAGAAAAAGGAACATTGCACACTCAATGACAACATGGCAAc                                                               >  2:49645/1‑141 (MQ=255)
                                                                     acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaag                                                               >  2:273691/1‑140 (MQ=255)
                                                                        ttttttaaaacaaattttgattaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                                                            >  2:275078‑M1/46‑141 (MQ=255)
                                                                           caactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                                                         <  1:186232/141‑1 (MQ=255)
                                                                              taatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc                                                      <  1:116316‑M1/102‑1 (MQ=255)
                                                                                  caaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAGTTTGTACATTGAAAACTGAATGACAATATGTCGACGTTAATTCCaaaa                                                  <  1:21814‑M1/106‑1 (MQ=255)
                                                                                  caaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGAGAATATGTGAACGTTCATtccacaa                                                  <  1:4753‑M1/106‑4 (MQ=255)
                                                                                          ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                          <  1:85497‑M1/114‑1 (MQ=255)
                                                                                              aaTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTCAATAGAATGTATAATTAATTTTTGTTGGTAATAAAAATGAACATTGAAATCTGAATGACAATATGTCAACGTTAATACCAAAAAACGTAACtata                                      >  2:280515/1‑141 (MQ=37)
                                                                                                     aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                               >  1:138410‑M1/17‑141 (MQ=255)
                                                                                                          ggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAGAAAAATCCAAAAAACGTAACTATAAGTTACAAACat                          >  1:247678‑M1/12‑141 (MQ=255)
                                                                                                            acaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTACAAAAAACGTAACTATAAGTTACAAACatta                        >  1:197078‑M1/10‑141 (MQ=255)
                                                                                                            acaagataaAAAAAGTAATTGACATAAATGTTGATAAAATGAATAATTAATTCTTGTCAGTAAGAAAAATGAACATTGAAAACTTAATGACCATATGTCAAAAAAAAAAACAAAAAACGTAACTATAAGTTACAAACatta                        >  2:231731‑M1/10‑141 (MQ=255)
                                                                                                              taattTTAAAAAGTTAGTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACattatt                      <  2:84251/141‑1 (MQ=255)
                                                                                                              aagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACattatt                      <  2:85290‑M1/134‑1 (MQ=255)
                                                                                                                attTTAAAATGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa                    >  1:208124/1‑141 (MQ=255)
                                                                                                                 ataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                   >  2:236400‑M1/5‑141 (MQ=255)
                                                                                                                     aaaaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAtt               <  2:61303/141‑1 (MQ=255)
                                                                                                                           ttATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg         <  1:96879/141‑1 (MQ=255)
                                                                                                                               tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATAAACTaa     >  2:191636/1‑141 (MQ=255)
                                                                                                                                gACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt    >  1:268879/1‑141 (MQ=255)
                                                                                                                                gACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt    >  2:220583/1‑141 (MQ=255)
                                                                                                                                  cTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCa  >  2:148383/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                  
AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCA  >  CP000730/512503‑512773

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: