Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,804,688 T→C 17.4% D142G (GAC→GGC)  PP_4233 ← oxidoreductase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,804,6880TC17.4% 44.1 / 3.7 23D142G (GAC→GGC) PP_4233oxidoreductase small subunit
Reads supporting (aligned to +/- strand):  ref base T (10/9);  new base C (2/2);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

GCTACGTACGAAGCCCATGTTCAACCCTCCTTGACCAGGCCGAGATCGTTGAGCACCTTGCGGGTGGCATCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGAAGCCTGGTGCGCAGTAGCTGCACTGGAAGGCGAACAGGTCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCG  >  NC_002947/4804556‑4804837
                                                                                                                                     |                                                                                                                                                     
gcTACGTACGAAGCCCATGTTCAACCCTCCTTGACCAGGCCGAGATCGTTGAGCACCTTGCGGGTGGCATCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCtcttc                                                                                                                                      <  2:225122/150‑1 (MQ=255)
                  gTTCAACCCTCCTTGACCAGGCCGAGATCGTTGAGCACCTTGCGGGTGGCATCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTgg                                                                                                                    <  1:58116/150‑1 (MQ=255)
                      aaCCCTCCTTGACCAGGCCGAGATCGTTGAGCACCTTGCGGGTGGCATCGTAGTCAAGCCCATAGCCTGTGCATAGCGGTATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCg                                                                                                                >  1:316662/1‑150 (MQ=255)
                          ctcctTGACCAGGCCGAGATCGTTGAGCACCTTGCGGGTGGCATCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTc                                                                                                            <  2:288724/150‑1 (MQ=255)
                                    aGGCCGAGATCGTTGAGCACCTTGCGGGTGGCATCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGctc                                                                                                                                         >  1:136623/1‑111 (MQ=255)
                                    aGGCCGAGATCGTTGAGCACCTTGCGGGTGGCATCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGctc                                                                                                                                         <  2:136623/111‑1 (MQ=255)
                                        cGAGATCGTTGAGCACCTTGCGGGTGGCATCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAg                                                                                              >  2:441621/1‑150 (MQ=255)
                                             tCGTTGAGCACCTTGCGGGTGGCATCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCa                                                                                         >  1:68133/1‑150 (MQ=255)
                                               gTTGAGCACCTTGCGGGTGGCATCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCAtt                                                                                       <  2:403616/150‑1 (MQ=255)
                                                     cACCTTGCGGGTGGCATCGTAGTAACGCACATAGCCTGTGCAGCGGCGGATCGCAGTACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGaa                                                                                 >  1:291141/1‑150 (MQ=255)
                                                                 ggCATCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGAAGCCTGGTGCGCa                                                                     <  2:68133/150‑1 (MQ=255)
                                                                     tCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGCCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCAc                                                                                                   <  1:306337/116‑1 (MQ=255)
                                                                     tCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGCCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCAc                                                                                                   <  1:306290/116‑1 (MQ=255)
                                                                     tCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGCCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCAc                                                                                                   >  2:306290/1‑116 (MQ=255)
                                                                     tCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGCCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCAc                                                                                                   >  2:306337/1‑116 (MQ=255)
                                                                        tagtaACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGAAGCCTGGTGCGCAGTAGCTg                                                              >  2:101196/1‑150 (MQ=255)
                                                                                               gCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGAAGCCTGGTGCGCAGTAGCTGCCCTGGAAGAGGAACAgcgctac                                       >  2:267959/1‑144 (MQ=255)
                                                                                                   gaagctttggtgACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGAAGCCTGGTGCGCAGTAGCTGCACTGGAAGGCGAACAGGTCGACGaa                                    <  2:291141/143‑1 (MQ=255)
                                                                                                            gtgACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGAAGCCTGGTGCGCAGTAGCTGCACTGGAAGGCGAACAGGTCGACGAACTTctgctg                           >  2:81037/1‑150 (MQ=255)
                                                                                                               aCCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGAAGCCTGGTGCGCAGTAGCTGCACt                                                          <  2:162776/116‑1 (MQ=255)
                                                                                                               aCCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGAAGCCTGGTGCGCAGTAGCTGCACt                                                          >  1:162776/1‑116 (MQ=255)
                                                                                                                 cGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGAAGCCTGGTGCGCAGTAGCTGCACTGGAAGGCGAACAGGTCGACGAACTTCTGCTGGATCg                      >  2:464047/1‑150 (MQ=255)
                                                                                                                    ggCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGAAGCCTGGTGCGCAGTAGCTGCACTGGAAGGCGAACAGGTCGACGAACTTCTGCTGGATCGGGt                   <  2:316662/150‑1 (MQ=255)
                                                                                                                                     tcctccACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGAAGCCTGGTGCGCAGTAGCTGCACTGGAAGGCGAACAGGTCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCg  <  1:441621/150‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                     
GCTACGTACGAAGCCCATGTTCAACCCTCCTTGACCAGGCCGAGATCGTTGAGCACCTTGCGGGTGGCATCGTAGTAACGCACATAGCCTGTGCAGCGGCAGATGTGGTGACCGAGGCTGGCCTCGATGCGGTCCTCCACTTCGCTCTTCTTCAACGGCTGGCGCTGGGCCGTTTCCACCAGCACGGTAGCGGCATTGACGAAGCCTGGTGCGCAGTAGCTGCACTGGAAGGCGAACAGGTCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCG  >  NC_002947/4804556‑4804837

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