Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,436,349 G→T 13.8% D8Y (GAC→TAC)  PP_3926 → short‑chain dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,436,3490GT13.8% 94.4 / 2.6 29D8Y (GAC→TAC) PP_3926short‑chain dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (16/9);  new base T (2/2);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

CGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCG  >  NC_002947/4436214‑4436484
                                                                                                                                       |                                                                                                                                       
cGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCg                                                                                                                                         <  2:400823/136‑1 (MQ=255)
  tttCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                                                                                                       >  1:62662/1‑136 (MQ=255)
                         aCCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCTACCTCGACGGCAAGATTGCCTTTGt                                                                                                                >  1:31609/1‑136 (MQ=255)
                                gCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCTACCTCGACGGCa                                                                                                                             >  2:226613/1‑116 (MQ=255)
                                gCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCTACCTCGACGGCa                                                                                                                             <  1:226613/116‑1 (MQ=255)
                                                ttGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCg                                                                                         <  1:488817/136‑1 (MQ=255)
                                                                       ccTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTc                                                                  <  1:453414/136‑1 (MQ=255)
                                                                             aGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGTATCGGCGAGGCCATCGCCCACTTGCTCGCgcag                                                            <  2:55580/136‑1 (MQ=255)
                                                                              gCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCGTCGGCGAGGCCATCGCCCACTTGCTCGCgcagc                                                           >  1:413503/1‑136 (MQ=255)
                                                                                acaGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACATGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCaa                                                         <  1:393263/136‑1 (MQ=255)
                                                                                     ccACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGc                                                    >  1:182620/1‑136 (MQ=255)
                                                                                       acgacgacgGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCaa                                                         >  1:24771/1‑129 (MQ=255)
                                                                                       acgacgacgGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCaa                                                         <  2:24771/129‑1 (MQ=255)
                                                                                       acgacgacgGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGccc                                                  <  1:38629/136‑1 (MQ=255)
                                                                                       acgacgacgGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGccc                                                  >  1:3882/1‑136 (MQ=255)
                                                                                          acgacgGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCAtg                                               <  1:38075/136‑1 (MQ=255)
                                                                                                 aTAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCgg                                        >  1:234078/1‑135 (MQ=255)
                                                                                                 aTAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGt                                        >  1:57177/1‑136 (MQ=255)
                                                                                                 aTAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGt                                        <  2:116364/136‑1 (MQ=255)
                                                                                                 aTAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAACGGGCCCATGTGATCGt                                        >  1:486595/1‑136 (MQ=255)
                                                                                                 aTAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCAGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGt                                        >  1:293700/1‑136 (MQ=255)
                                                                                                       aGGAAAACAGCATGTCCAAGACCCACCTGTTCTACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAg                                  <  2:31609/136‑1 (MQ=255)
                                                                                                           aaaCAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGc                              >  2:397544/1‑136 (MQ=255)
                                                                                                                            cccACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGccg             >  2:48419/1‑135 (MQ=255)
                                                                                                                             ccACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCcag            <  2:9632/136‑1 (MQ=255)
                                                                                                                                  tGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGt       >  2:382036/1‑136 (MQ=255)
                                                                                                                                  tGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGt       >  1:432156/1‑136 (MQ=255)
                                                                                                                                  tGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGt       >  1:245686/1‑136 (MQ=255)
                                                                                                                                  tGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGt       >  1:16846/1‑136 (MQ=255)
                                                                                                                                  tGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGt       >  1:155858/1‑136 (MQ=255)
                                                                                                                                       gACCTCGACGGCAAGATTGCCTTTGTTTCATTAGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACATGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCg  >  2:227036/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                       
CGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCG  >  NC_002947/4436214‑4436484

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