Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,310,045 G→A V75V (GTC→GTT kch ← K(+) channel Kch

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,310,0450GA95.5% 63.8 / ‑5.7 22V75V (GTC→GTTkchK(+) channel Kch
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (9/12);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

AGGTTGAGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCAAAATCCAATACTAAATTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGGACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTTAAAAATAAGCAATCCATTTAAAAAGACGGCCAGCGCGAGAATGTCGTGCCGAAGGGTCACCCATAAATTTGTTGCTGTTTGTTTGAATGTAGCCCAGTGACTCACTGA  >  NC_000913/1309804‑1310272
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                   
aGGTTGAGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCAAAATCCAATACTAAATTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACaagaaa                                                                                                                                                                                                                             <  1:263505/250‑1 (MQ=255)
   ttGAGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCAAAATCCAATACTAAATTTCAGCCAGGGATAAAAATGCAGGGTGAAAATTAGCGCTCTCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGACCAAGAAAGACg                                                                                                                                                                                                                         >  1:330751/1‑251 (MQ=255)
                                          aaaTTGTCCCGGCTGCGGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCAAAATCCAATACTAAATTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGgataga                                                                                                                                                                                  <  2:133524/251‑1 (MQ=255)
                                                    ggCTGCGGCGCTACTGTGGGAGAAGTCTTTGCTCAGTATCAGCAAAAACACCAGCGTAAAAAGGCAAAATCCAATACTAAATTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGTGGCGTTTGCCAGCAAGGAGAGATCCAACCCTg                                                                                                                                                                        <  2:69123/251‑1 (MQ=255)
                                                                                               aaaaaCACCAGCGTAAAAATGCAAAATCCAATACTAAATTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGAt                                                                                                                                                                                 <  2:81025/199‑1 (MQ=255)
                                                                                               aaaaaCACCAGCGTAAAAATGCAAAATCCAATACTAAATTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGAt                                                                                                                                                                                 >  1:81025/1‑199 (MQ=255)
                                                                                                 aaaCACCAGCGTAAAAATGCAAAATCCAATACTAAATTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATAc                                                                                                                           >  1:218052/1‑251 (MQ=255)
                                                                                                                 aTGCAAAATCCAATACTAAATTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCAGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTNaa                                                                                                            >  2:419967/1‑250 (MQ=255)
                                                                                                                                 tAAATTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTTACAAATAAGCAATCCAttt                                                                                           >  1:268264/1‑251 (MQ=255)
                                                                                                                                 tAAATTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTTAAAAAt                                                                                                        <  2:34317/238‑1 (MQ=255)
                                                                                                                                 tAAATTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACAGACATACCATAGATTGTTTTAAAAAt                                                                                                        >  1:34317/1‑238 (MQ=255)
                                                                                                                                  aaaTTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTTAAAAATAAGCAATCCAttt                                                                                           >  2:259118/1‑250 (MQ=255)
                                                                                                                                       aaaGCCAGGGATAAAAATGCAGCGTGGAAAATAGCGCTATCAACACCAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCCATGGAGTTCAGAACAAGAAAGACGCCGAGCTTAAAGCGTGGGTCGTTTTACAGCAAGGACAGATCCGACACTGAAAATGCTTTGCAATGGAAAAAATCAAGCAAACTGACAGACATTCCATAGATTTTTTTTAACATAAGCCATCCCTCTaaaaa                                                                                      <  2:100668/248‑1 (MQ=255)
                                                                                                                                                 aTAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTTAAAAATAAGCAATCCATTTAAAAAGACGGCCAgcg                                                                           <  2:2362/251‑1 (MQ=255)
                                                                                                                                                           aGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGATCTGTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTTAAAAATAAGCAATCCATTTAAAAAGACGGCCAGCGCGAGAATGTc                                                                 <  2:330751/251‑1 (MQ=255)
                                                                                                                                                                gtAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTTAAAAATAAGCAATCCATTTAAAAAGACGGCCAGCGCGAGAATGTCGTGcc                                                            <  1:28517/251‑1 (MQ=255)
                                                                                                                                                                  aaaTTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTTAAAAATAAGCAATCCATTTAAAAAGACGGCCAGCGCGAGAATGTCGTGCCg                                                           <  1:259118/250‑1 (MQ=255)
                                                                                                                                                                           cTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGGACAAGAAAGACGCCGAGCATAAATAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAat                                                                                                                                                  <  2:218934/154‑1 (MQ=255)
                                                                                                                                                                                                 gATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTTAAAAATAAGCAATCCATTTAAAAAGACGGCCAGCGCGAGAATGTCGTGCCGAAGGGTCACCCATAAATTTGTTGCtgtttgtt                           >  1:138797/1‑251 (MQ=255)
                                                                                                                                                                                                    tGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAAGATCAAGCAAATTGACCGACATACCATAGATGGTGTTAAAAATAAGCAAGCCATTTAAAAAGACGGCCAGCGCGAGAATGTCGTGCCGAAGGGTCACCGATAAAGTTGTTGCtgttggttag                         <  2:55801/250‑3 (MQ=255)
                                                                                                                                                                                                               gCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTTAAAAATAAGCAATCCATTTAAAAAGACGGCCAGCGCGAGAATGTCGTGCCGAAGGGTCACCCATAAATTTGTTGCTGTTTGTTTGAATGTAGccc               >  1:340302/1‑249 (MQ=255)
                                                                                                                                                                                                                          aaCAGTAAGCCAATGGAGTTCAGAACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTTAAAAATAAGCAATCCATTTAAAAAGACGGCCAGCGCGAGAATGTCGTGCCGAAGGGTCACCCATAAATTTGTTGCTGTTTGTTTGAATGTAGCCCAGTGACTCACTGa  <  2:92744/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                   
AGGTTGAGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCAAAATCCAATACTAAATTTCAGCCAGGGATAAAAATGCAGGGTGTAAATTAGCGCTATCAACAACAAAATGATACTGATTGCCCATGCGAGCTTTGCCCGGAACAGTAAGCCAATGGAGTTCAGGACAAGAAAGACGCCGAGCATAAAGAGTGGGGCGTTTGCCAGCAAGGAGAGATCCAACTCTGAAAATGCTTTGATATGGAAAATATCAAGCAAATTGACCGACATACCATAGATTGTTTTAAAAATAAGCAATCCATTTAAAAAGACGGCCAGCGCGAGAATGTCGTGCCGAAGGGTCACCCATAAATTTGTTGCTGTTTGTTTGAATGTAGCCCAGTGACTCACTGA  >  NC_000913/1309804‑1310272

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