Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,624,902 A→G 16.7% L210P (CTG→CCG) ‡ PP_3195 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,624,9020AG16.7% 50.5 / 4.3 24L210P (CTG→CCG) ‡PP_3195hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/12);  new base G (2/2);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

TCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGCCTTCATACGGGATCT  >  NC_002947/3624764‑3625043
                                                                                                                                          |                                                                                                                                             
tCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTGCCTGGCTGGCcaggcgaggc                                                                                                                                   <  2:611227/150‑6 (MQ=255)
    cACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCTGGCcaggcgagg                                                                                                                                    >  1:600106/1‑142 (MQ=255)
    cACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCTGGCcaggcgagg                                                                                                                                    <  2:600106/146‑5 (MQ=255)
    cACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCaggcaggcaggc                                                                                                                                <  2:193507/150‑1 (MQ=255)
                      ccACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttaattt                                                                                                              >  2:525216/1‑150 (MQ=255)
                         cTATGCGGCTGGTCCAGTGTTCGCGGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCt                                                                                                           <  2:80266/150‑1 (MQ=255)
                           aTGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTg                                                                                                         <  1:76872/150‑1 (MQ=255)
                            tGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCttag                                                                                                        >  2:85000/1‑148 (MQ=255)
                             gCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGt                                                                                                       <  1:299399/150‑1 (MQ=255)
                                     cccAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGTCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttca                                                                                                                 <  1:457365/132‑1 (MQ=255)
                                     cccAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGTCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttca                                                                                                                 >  2:457365/1‑132 (MQ=255)
                                                  gTCATGGCGTGCCGTCATGTCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTg                                                                                  <  2:404453/150‑1 (MQ=255)
                                                          gTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTGTGTTGAATCCAGCCCGTTGACTTGCACcagca                                                                          <  1:119508/150‑1 (MQ=255)
                                                                   tGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCg                                                                 >  2:346250/1‑150 (MQ=255)
                                                                              cGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACg                                                      >  2:3720/1‑150 (MQ=255)
                                                                                aaGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGcc                                                    <  1:369549/150‑1 (MQ=255)
                                                                                          ccGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGt                                                                     <  2:253303/123‑1 (MQ=255)
                                                                                          ccGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGt                                                                     >  1:253303/1‑123 (MQ=255)
                                                                                          ccGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTgcgcgc                                          <  1:3720/150‑1 (MQ=255)
                                                                                                 cGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACcagca                                   >  2:268219/1‑150 (MQ=255)
                                                                                                           gACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTAGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGagag                         <  1:85000/150‑1 (MQ=255)
                                                                                                              aaTCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCg                                                                 >  2:55790/1‑107 (MQ=255)
                                                                                                              aaTCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCg                                                                 <  1:55790/107‑1 (MQ=255)
                                                                                                                                  gcctGGCTGGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGCCTTCATACGGGAtct  >  1:482922/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                             
TCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGCCTTCATACGGGATCT  >  NC_002947/3624764‑3625043

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