| 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 | ||
| * | CP000730 | 512,643 | 0 | A | G | 51.0% | ‑4.9 / 58.9 | 45 | intergenic (+711/‑150) | recR/rrs1 | recombination 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: ‑ |