Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,702,363 G→A A36V (GCC→GTC)  recO ← recombination mediator protein RecO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,702,3630GA95.7% 65.8 / ‑3.8 23A36V (GCC→GTC) recOrecombination mediator protein RecO
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (9/13);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

GTGCAAGTAATCGAAAAAAAGTTCAGAGAAGCGCGTCTCGTATTCCAGTACGCGGGAGAGAAGTTCGTTGATGTACAGACCGCTGTAAAGCGTGATACCGCTTAATGGCAGCGCCAGCGAGACGGCTTCAGCACTGCGCAGCGTTTTGACTTCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGGCAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTTACGCGCTT  >  NC_000913/2702125‑2702608
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
gTGCAAGTAATCGAAAAAAAGTTCAGAGAAGCGCGTCTCGTATTCCAGTACGCGGGAGAGAAGTTCGTTGATGTACAGACCGCTGTAAAGCGTGATACCGCTTAATGGCAGCGCCAGCGAGACGGCTTCAGCACTGCGCAGCGTTTTGACTTCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGacgc                                                                                                                                                                                                                                            >  1:35595/1‑250 (MQ=255)
             aaaaaaaGTTCAGAGAAGCGCGTCTCGTATTCCAGTACGCGGGAGAGAAGTTCGTTGATGTACAGACCGCTGTAAAGCGTGATACCGCTTAATGGCAGCGCCAGCGAGACGGCTTCAGCACTGCGCAGCGTTTTGACTTCGCCCCGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATtc                                                                                                                                                                                                                              <  1:218556/251‑1 (MQ=255)
                  aaGTTCAGAGAAGCGCGTCTCGTATTCCAGTACGCGGGAGAGAAGTTCGTTGATGTACAGACCGCTGTAAAGCGTGATACCGCTTAATGGCAGCGCCAGCGAGACGGCTTCAGCACTGCGCAGCGTTTTGACTTCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCg                                                                                                                                                                                                                         >  1:172148/1‑251 (MQ=255)
                                        gagcccAGTCCGCGGGAGATAAGTTCGTGGATGTACAGACCGCTGTAAAGCTTGATACCGCTTAATGGCAGCGCCAGCGAGCCGGCTTCAGCACTGCGCAGCGTTTTGACTGCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGTCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCt                                                                                                                                                                                                   <  2:87582/247‑1 (MQ=255)
                                         atccaaGTACGCGGGAGGGAAGTTCGTTGATGTACCTCACGCTGTAAAGCGTGAAACTGCGTATTGGCAGCGCGAGCTAGACGGCTGCATCCTGTCGCAGCGTTTTGACTTCGCCACGTCCGCCAAAACGTAGCAAGAGAGGTGTGAAAGGCTGTAATGCACCTGGCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACACCCAGACGCACGCGACCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCt                                                                                                                                                                                                   <  2:111081/245‑1 (MQ=255)
                                           tCCAGTACGCGGGAGATAAGTTCGTTGATGTACAGACCGCTGTAAAGCGGGAGACCGCTTAATGGCAGCGCCAGCGAGACGGCTTCAGCACTGCGCATCGTTTTGACTTCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCATGGTAGAGCGTTTAGAGCGTGCGCCTTTGGCAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTgg                                                                                                                                                                                                 <  2:100571/250‑1 (MQ=255)
                                                        agagaAGTTCTTTGATGTACAGACCGCTGAAAAGCGTGATACCGCTTAATGGCAGCGCCAGCGAGACGGCTTCAGCACTGCGCAGCGTTTGGACTTCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACggg                                                                                                                                                                                   <  1:277184/251‑1 (MQ=255)
                                                             aGTTCGTTGATGTACAGACCGCTGTAAAGCGTGATACCGCTTAATGGCAGCGCCAGCGAGACGGCTTCAGCACTGCGCAGCGTTTTGACTTCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACt                                                                                                                                                                              <  1:233218/251‑1 (MQ=255)
                                                                tCGTTTATGTACAGACCGCTGTAAAGCGTAATACCGCTTAATGGAAGCGCCAACGAGTCGGCTTCAGCACTTCGCAGCGTTTTGACTTCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTAt                                                                                                                                                                            <  2:7861/250‑1 (MQ=255)
                                                                                                                      gagaCGGCTTCAGCACTGCGCAGCGTTTTGACTTCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCCTCGGAGTTACTCTTa                                                                                                                               >  2:236622/1‑241 (MQ=255)
                                                                                                                      gagaCGGCTTCAGCACTGCGCAGCGTTTTGACTTCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTa                                                                                                                               <  1:236622/241‑1 (MQ=255)
                                                                                                                      gagaCGGCTTCAGCACTGCGCAGCGTTTTGACTTCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCa                                                                                                                     <  1:181803/251‑1 (MQ=255)
                                                                                                                                    ccTGCGCAGCGTTTTGACTTCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACt                                                                                                       <  1:107159/250‑1 (MQ=255)
                                                                                                                                     cgtgcgcaggttgtagACTTCGCCACTCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCCCCTTTCAGGGTCGAGCGTTTAGAGCGGGCGCCTTTTACAACCAGACGCCCGCGCCCCGATTCCGCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGGAACCGAGACt                                                                                                       <  2:194340/236‑1 (MQ=255)
                                                                                                                                            ggcgtttggttcgcgccccgcccGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCGTGGACAACCAGACGCACGCGCCCCGAGTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTGCGCTCCACGGTCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGGGACTGCGCAGt                                                                                                <  2:35595/238‑1 (MQ=255)
                                                                                                                                                                       aaaCGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGAGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCAAACCGGAt                                                                     >  1:249054/1‑251 (MQ=255)
                                                                                                                                                                         aCGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGAttt                                                                   <  1:93190/251‑1 (MQ=255)
                                                                                                                                                                              gCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACtt                                                              >  1:124803/1‑251 (MQ=255)
                                                                                                                                                                                aaGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGAGTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCATGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACtttt                                                            <  1:186617/251‑1 (MQ=255)
                                                                                                                                                                                                   gTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGa                                          >  2:294472/1‑250 (MQ=255)
                                                                                                                                                                                                                cAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCCCGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCACCGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCGCAGTTACAGGTGAACAGGCACTTCGa                             >  2:205823/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ccTTTGACAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTTCCGCGCt   >  2:294807/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cTTTGACAACCAGAAGCACGCGCCCCGCTTCATCCGTGAAGACGTCCAGCATCAGGCTGGCTTCGCACCACGGGCTACTATGCCGGCCAACTGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGGGCGCGTTCATCGTCGGCCCAGCCGGATCTCACTTTTACCACCCGTTCAAGGTGAACACGAGCTTCGCACATTTCCTGCATGGTTTTACGAGCtt  >  2:181976/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
GTGCAAGTAATCGAAAAAAAGTTCAGAGAAGCGCGTCTCGTATTCCAGTACGCGGGAGAGAAGTTCGTTGATGTACAGACCGCTGTAAAGCGTGATACCGCTTAATGGCAGCGCCAGCGAGACGGCTTCAGCACTGCGCAGCGTTTTGACTTCGCCACGCCCGCCAAAACGTAGCAAGAGAGGGGTGAAAGGCTGTAATGCACCTTTCAGGGTAGAGCGTTTAGAGCGTGCGCCTTTGGCAACCAGACGCACGCGCCCCGATTCCTCCGTGAAGACGTCCAGCATCAGGCTGGTTTCGCTCCACGGGCGACTATGCAGGACAAATGCGCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTTACGCGCTT  >  NC_000913/2702125‑2702608

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