Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,259,361 G→A P97S (CCA→TCA)  tdcG ← L‑serine deaminase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,259,3610GA81.0% 32.1 / 2.4 22P97S (CCA→TCA) tdcGL‑serine deaminase III
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (10/7);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGT  >  NC_000913/3259116‑3259595
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          
cagcagTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGTCAATAAACCCGCCGCCGACAGAGTAATATGTTTTACTTAAGAGCTCTTCCTGTCCCTTCCAGGCAGGGATCCGCATTCCGTTCTCATGGCGAGGCAACATTGCGGGATGGAAGATAATGTGCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGACAgn                                                                                                                                                                                                                                        <  1:41179/250‑2 (MQ=255)
                                           tCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGAAATAGGTTTTACTTAATAGCTCTTCCTGTCACTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGAGAGGCAACATTTCGGTATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAACTGCAGGGATCtcat                                                                                                                                                                                             >  1:56494/1‑250 (MQ=255)
                                                                    gtTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTgg                                                                                                                                                                                                                                            <  1:133923/178‑1 (MQ=255)
                                                                            tcGACAATAAACCCGCCGCCGACAGAGTACTAGGTTTTCATTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACAGTTCGGGATGGAAGATAATGTTCTTTGATACAGGAAAATCAACAATATGCGCAACAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGAGCTCATCAATGACAACATCCTGCGGACTGTTTCCTGcn                                                                                                                                                             <  1:58056/249‑2 (MQ=255)
                                                                              gACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTACTGTCCCTTCCAGGCAGCGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTAAAGGAAAATCAACAATCTGCGCACCAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGa                                                                                                                                                                                       >  1:234339/1‑221 (MQ=255)
                                                                                                                     tAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCg                                                                                                                  >  1:4389/1‑251 (MQ=255)
                                                                                                                                              gCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCGACAATATGCGCACCAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCTATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCg                                                                                                                  <  1:319928/226‑1 (MQ=255)
                                                                                                                                                                    cATNNCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGACAGccgcacgc                                                                                                                                                                                                                                 <  2:19771/93‑1 (MQ=255)
                                                                                                                                                                       ggCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAATTGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGCggn                                                                 <  1:61222/250‑2 (MQ=255)
                                                                                                                                                                             gCAACATTTCGGGATGGAAGATAATGTTCATTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACAGCCTGCGGACTGTTTCCTGCCAGTCCCATGATGAGGGCGACATCCGGGGCATGGCCTTGGCCCGTCAGTTACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGCGGTTAATAAg                                                          <  1:280157/251‑1 (MQ=255)
                                                                                                                                                                                         gATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGa                                                                                                           >  1:256900/1‑190 (MQ=255)
                                                                                                                                                                                           tGGAAGTTAATGTTCTTTGCTACAGGAAAATCAACAATATCCGCACCAGATGCCACTGACAGCCGCACGGTGCGCGTTACTAACTCTATAAATGCAAGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCAGTACAGATACAACACAATATGGCTCGTCAGGGGTAATAAACCGCTACTTTCCa                                             >  1:74949/1‑250 (MQ=255)
                                                                                                                                                                                            ggAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGCGGTTAATAAGCCGCTACTTTCCAGc                                           <  1:274628/251‑1 (MQ=255)
                                                                                                                                                                                            ggAAGATAATGTTATTTGCTACAGGAAAATCAACAATCTGCGCACCAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGCCTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGCGGTTAATAAGCCGCTACTTTCCAGc                                           >  1:268285/1‑251 (MQ=255)
                                                                                                                                                                                             naaGATAATGTTCTTTGCTACAGGCAAATCAACAATATGCGCATCAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGCGGTTAATAAGCCGCTACTTTCCAGcc                                          >  1:85460/2‑251 (MQ=255)
                                                                                                                                                                                                   aaTGTTCTTTGCTACAGGAAAAACAACAATATGCGCGCCAGATGGCACGGTCGGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTg                                                                                                                                  <  1:302470/157‑1 (MQ=255)
                                                                                                                                                                                                         cTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACCCTATAAATGCAGGGATCTCATCAATGACAACAGCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGAc                                                                                                          <  1:46818/175‑1 (MQ=255)
                                                                                                                                                                                                              ntACAGGAAAATCAACAATATGCGCACCAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACttttt                         >  1:80240/2‑251 (MQ=255)
                                                                                                                                                                                                                 cAGGAAAATCAACAATATGCGCACCAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAatat                                                                              >  1:123702/1‑195 (MQ=255)
                                                                                                                                                                                                                           aacaaTATGCGCACCAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGcc                                                                                                                                                             <  1:97972/106‑1 (MQ=255)
                                                                                                                                                                                                                           aacaaTATGCGCACCAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAGACGTTTGCCGGcattcatt             <  1:107117/250‑1 (MQ=255)
                                                                                                                                                                                                                                     gcACCAGATGCCACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATAAATGAGAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGCAGTTAATAAGACGCTACTTTCCCGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCACCACGGt  >  1:229933/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ccACTGACAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGACCaa                                    >  1:73214/1‑207 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          
CAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGT  >  NC_000913/3259116‑3259595

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