Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 585,643 C→A 21.0% M377I (ATG→ATT PP_0496 ← sodium/alanine symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947585,6430CA21.0% 44.7 / 2.8 19M377I (ATG→ATTPP_0496sodium/alanine symporter
Reads supporting (aligned to +/- strand):  ref base C (9/6);  new base A (2/2);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

GCAGGGCCAGGGCCACCAGGTTGACGAAGGCCAGGCAGGTCATGGTGATGTCGGCGAAGGCGAACACGGTCGACAGGTCTTGCATCGAACCCCACACGACCAGTGCCAGCACCAGGCCACGGAACACCATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAGGTAGTAGTTGTAGAGGATGCAGGTGAAGACGAACAGCGACAGCGCCACGCTGACGAACACGCGCCCCCAGTCACC  >  NC_002947/585516‑585768
                                                                                                                               |                                                                                                                             
cctgggccagggccaCCAGGTTGACGACGGCCAGGCAGGTCATGGTGATGTCGGCGAAGGCGAACACGGTCGACAGGTCTTGCATCGAACCCCACACGACCAGTGCCAGCACCAGGCCACGGAACACCATCAacac                                                                                                                       <  1:177263/133‑1 (MQ=255)
     gccagggccaCCAGGTGGACGAAGGCCAGGCTGGTCTTGGTGATGTCGGCGAAGGCTAACACGGTCGACAGGTCTTGCATCGAACCCCACACGACCAGTGCCAGCACCAGGCCACGGAACACCATCAACACTAcgc                                                                                                                  <  2:207725/136‑1 (MQ=255)
                  ggTTGACGAAGGCCAGGCAGGTCATGGTGATGTCGGCGAAGGCGAACACGGTCGACAGGTCTTGCATCGAACCACACACGACCAGTGCCAGCACCAGGCCACGGAACACCATCAACACTACGCGGTTACGGCTGAg                                                                                                     >  1:284900/1‑136 (MQ=255)
                       aCGAAGGCCAGGCAGGTCATGGTGATGTCGGCGAAGGCGAACACGGTCGACAGGTCTTGCATCGAACCCCACACGACCAGTGCCAGCACCAGGCCACGGAACACCATCAACACTACGCGGTTACGGCTGAGGAAtt                                                                                                >  1:105585/1‑136 (MQ=255)
                                         aTGGTGATGTCGGCGAAGGCGAACACGGTCGACAGGTCTTGCAACGAACCCCACACGACCAGTGCCAGCACCAGGCCACGGAACACCATCAACACAACGCGGTTACGGCTCAGGAATAGCAGGCTGTTCTCGCCGa                                                                              >  1:354576/1‑136 (MQ=255)
                                           ggTTATGTCGGCGAAGGCGAACACGGTCGACAGGTCTTGCATCGAACCCCACACGACCAGTGCCAGCACCAGGCCACGGAACACCATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAgg                                                                            <  2:301144/136‑1 (MQ=255)
                                           ggTGTTGTCGGCGAAGGCGAACACGGTCGACAGGTCTTGCATCGAACCCCACACGACCAGTGCCGGCACCAGGCCACGGAACACCATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAgg                                                                            <  1:302096/136‑1 (MQ=255)
                                                               acacGGTCGACAGGTCTTGCATCGAACCCCACACGACCAGTGCCAGCACCAGGCCACGGAACACCATCAACACTACGCGGTTACGGCTGAGGAATTGTAGGCTGTTCTCGCCGAGGTAGTAGTTGTAGAGGATGCa                                                        >  1:179046/1‑136 (MQ=255)
                                                                    gTCGACAGGTCTTGAATCGAACCCCACACGACCAGTGCCAGCACCAGGCCACGGAACACCATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAGGTAGTAGTTGTAGAGGATGCAGGTGa                                                   <  2:382520/136‑1 (MQ=255)
                                                                                               acGACCAGTGCCAGCACCAGGCCACGGAACACAATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAgg                                                                            >  2:56327/1‑84 (MQ=255)
                                                                                               acGACCAGTGCCAGCACCAGGCCACGGAACACAATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAgg                                                                            >  2:53743/1‑84 (MQ=255)
                                                                                               acGACCAGTGCCAGCACCAGGCCACGGAACACAATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAgg                                                                            <  1:56327/84‑1 (MQ=255)
                                                                                               acGACCAGTGCCAGCACCAGGCCACGGAACACAATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAgg                                                                            <  1:53743/84‑1 (MQ=255)
                                                                                                          cAGCACCAGGCCACGGAACACCATCACCACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCGCGCCGAGGTAGTAGTTGTAGAGGATGCAGGTGAAGACGAACAGCGACAGCGCCACGCTGACGAACACgcgc             <  1:456185/136‑1 (MQ=255)
                                                                                                          cAGCACCAGGCCACGGAACACCATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAGGTAGTAGTTGTAGAGGATGCAGGTGAAGACGAACAGCGACAGCGCCACGCTGACGAACACgcgc             >  2:432274/1‑136 (MQ=255)
                                                                                                                cAGGCCACGGAACACCATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAGGTAGTAGTTGTAGAGGATGCAGGTGAAGACGAACAGCGACAGCGCCACGCTGACGAACACGCGCCCCCAg       >  1:270788/1‑136 (MQ=255)
                                                                                                                    ccACGGAACACCATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAGGTAGTAGTTGTAGAGGATGCAGGTGAAGACGAACAGCGACAGCGCCACGCTGACGAACACGCGCCCCCAGTCAc   >  2:378740/1‑136 (MQ=255)
                                                                                                                    ccACGGAACACCATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAGGTAGTAGTTGTAGAGGATGCAGGTGAAGACGAACAGCGACAGCGCCACGCTGACGAACACGCGCCCCCAGTCAc   >  1:354568/1‑136 (MQ=255)
                                                                                                                     cACGGAACACCATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAGGTAGTAGTTGTAGAGGATGCAGGTGAAGACGAACAGCGACAGCGCCACGCTGACGAACACGCGCCCCCAGTCAcc  >  2:204067/1‑136 (MQ=255)
                                                                                                                               |                                                                                                                             
GCAGGGCCAGGGCCACCAGGTTGACGAAGGCCAGGCAGGTCATGGTGATGTCGGCGAAGGCGAACACGGTCGACAGGTCTTGCATCGAACCCCACACGACCAGTGCCAGCACCAGGCCACGGAACACCATCAACACTACGCGGTTACGGCTGAGGAATTGCAGGCTGTTCTCGCCGAGGTAGTAGTTGTAGAGGATGCAGGTGAAGACGAACAGCGACAGCGCCACGCTGACGAACACGCGCCCCCAGTCACC  >  NC_002947/585516‑585768

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