Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,409 G→A V29V (GTG→GTA deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,4090GA96.0% 67.8 / ‑3.1 25coding (87/780 nt)deoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/9);  minor base . (0/1);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

AATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGT  >  NC_000913/4617172‑4617630
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                              
aaTTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGtc                                                                                                                                                                                                                  <  2:26497/251‑1 (MQ=255)
                   cATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACCCCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGt                                                                                                                                                                                               >  1:136003/1‑251 (MQ=255)
                                     ttGTACGTAGATATCCTTCATTGTGCCGTATCTCGCAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGCATTTTATTGCCGACAAGCCAGGAGAATGAAACGCCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATCGATGGACCTGACCACCCTGAATGACGCCGACACCGACGAGAAAGTAATTGCTCTATGTCATCACGCCCAAACCCCGGTCGGCACTACCGCCGCTAt                                                                                                                                                                             >  2:109881/1‑251 (MQ=255)
                                      tGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATc                                                                                                                                                                            >  2:23897/1‑251 (MQ=255)
                                        tAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTg                                                                                                                                                                          >  2:103105/1‑251 (MQ=255)
                                                            gaaggtaaaTCGAAGTGGGTTGCGGTGAAGAATTGAGAAAAATAAAAAACTCGAAAGGTGATTTTGATGGGCGAAAAGCCAGGAAAATGAAATTACTGATCTAAAAGCAAACAGCCTGCGTGCACTGAAATTGATGGACCTGACCCCCCTAAAGAACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGCCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATcc                                                                                                                                                      <  1:185664/243‑1 (MQ=255)
                                                              aTGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCg                                                                                                                                                    >  1:176271/1‑251 (MQ=255)
                                                                         ggtgtgtTGCGGAGTAGATGTTAAAAAACTAACAAACTCGCAAGTTAAATTTTATTGGCGACAAGCCAGGAGAATAAAATGACTAATCTGAAAGCAAGCAGCCGGCGTGCATTGAAATTGAGGGACCGGCCCACCCTAAATGACGCCGACACCGACGAGAAAGTATTTCCTCTGTTTCATCAGGCCAAACCCCCGTTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCaa                                                                                                                                         <  1:114402/250‑1 (MQ=255)
                                                                                         gATGTTAGAATACTAACAAACTCGCAAGGTGATTTTTATTGGCGACAAGCCAGGAGCATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGc                                                                                                                           >  2:185852/1‑249 (MQ=255)
                                                                                                         aaaaCTCGAAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATTACTGATCTGAAAGCAATAAGCCTGCGTGCACTGAAATTGATGGACCTGACCCCCCTGAATGACGACGACCCCGACGAGAAAGTAATTGTTTTGTGTCATCGGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATcc                                                                                                         <  1:21014/250‑1 (MQ=255)
                                                                                                                   aGGTGAATTTTCTTGGCGACACGCCAGGAGAATGAAATGACGGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGCCACCGACGAGAATGTAATTGATCTGTGTCATCAGGCACAAACCCCGGTCGACAATACCGCCGCTATCTGTATCAATACTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGCCAATCCGTAttgatac                                                                                               >  2:197886/1‑245 (MQ=255)
                                                                                                                   aGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTCAAAGAAAGCAGCCTGCGTGCACTGCAATTGATGGACCTGACCAACCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATAAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAACGAGCAGGGCACCCCGGAAATCCGTATTGCTAc                                                                                               >  1:172401/1‑251 (MQ=255)
                                                                                                                                       aaGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACg                                                                           >  2:134937/1‑251 (MQ=255)
                                                                                                                                           cAGGAGAATGAAATGACTGATCTAAAAGCAAGCAGCCGGCGTGCACTAAAATTGATGGACCTGACCCCCCTAAATGACGACGACACCGACGAAAAAGTAATTCCTCTGTGTCATCAGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTaa                                                                       <  1:23897/251‑1 (MQ=255)
                                                                                                                                                  atgaaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGACCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCCCCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACgatgac                                                                >  1:130476/1‑251 (MQ=255)
                                                                                                                                                             aatcggaaagcaaaaggcatGCGGGCACTGAAATTGATGGACCTGCCCCCCCTAAAAGACGCCACCCCCGACAAGAAAGTAATGGCTCTGTGTCAACAGGCAAAACCCCCGGCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACATCCCACACGGTAACGATgacatcgaaatc                                                       <  2:136003/230‑1 (MQ=255)
                                                                                                                                                                   aaaGCAACAAGCCTGCGTGCCGGGAAATTGATGGACCTGCCCCCCCTGAATGACAACGACACCGCCGAGAAAGTAATTGCTCTGTTTCACCAGCCCAAACCCCCGGCCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGc                                               <  2:15109/251‑1 (MQ=255)
                                                                                                                                                                          gcagcCTGCGTGCACTGAAATTGCTGGACCTGACCCCCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTCTCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCACAACACTGAAAGCGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGCCATCGCCATCGCGCTTGCAGAAAcc                                        >  1:12427/1‑251 (MQ=255)
                                                                                                                                                                                tGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGCACCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCg                                  >  1:226463/1‑251 (MQ=255)
                                                                                                                                                                                      cACTGAAAAGGATGGACCGGCCAACCCGGAATGCCGACGACACCGACTAGAAAGTAATTGCTCGGTGCCATCAGGCCAACTCCCGGGGCGCCAAAACCGCCGCTATCTGTATCTATCCGCGCGTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACAGCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                            <  2:12427/251‑1 (MQ=255)
                                                                                                                                                                                                    gACCTGACCCCCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                                          <  1:223821/223‑1 (MQ=255)
                                                                                                                                                                                                    gACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                                          >  2:223821/1‑223 (MQ=255)
                                                                                                                                                                                                    gACCTGACAACCCTGAATGACGACGACACCGAAGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGCCGGCAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGACAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAAACAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGccgg              >  2:163162/1‑248 (MQ=255)
                                                                                                                                                                                                         gACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATACGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCcgatga           >  1:41477/1‑249 (MQ=255)
                                                                                                                                                                                                                cgaaaTGACGACAACACCGACGAAAAAGTAATTGCTCTGGGTCATCAGGCAAAACCCCCGGTCGGCAAAACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGt  <  1:130736/248‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                              
AATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGT  >  NC_000913/4617172‑4617630

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