Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 638,745 A→T 15.9% L16Q (CTG→CAG)  PP_0550 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947638,7450AT15.9% 69.1 / 6.5 31L16Q (CTG→CAG) PP_0550hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (14/12);  new base T (3/2);  total (17/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.63e-01

TGCGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGCTTCTTCGGCTT  >  NC_002947/638610‑638892
                                                                                                                                       |                                                                                                                                                   
tGCGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCg                                                                                                                                       <  1:499171/150‑1 (MQ=255)
  cGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGc                                                                                                                                     <  1:241958/150‑1 (MQ=255)
    gaaggaaTCCCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCaa                                                                                                                                   <  1:370082/150‑1 (MQ=255)
           tACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGAc                                                                                                                            <  1:228603/150‑1 (MQ=255)
                        cccACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCTGACCGGTgaa                                                                                                                                           <  2:476415/122‑3 (MQ=255)
                        cccACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCTGACCGGTgaa                                                                                                                                           >  1:476415/1‑120 (MQ=255)
                         ccACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCa                                                                                                              >  1:679411/1‑150 (MQ=255)
                                                   tCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTac                                                                                    >  2:104202/1‑150 (MQ=255)
                                                   tCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGAGGGTGGGATCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTac                                                                                    >  2:174660/1‑150 (MQ=255)
                                                                      tAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGc                                                                 >  1:245134/1‑150 (MQ=255)
                                                                                     gggTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCt                                                  >  2:406159/1‑150 (MQ=255)
                                                                                     gggTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCt                                                  >  1:781236/1‑150 (MQ=255)
                                                                                             gtggGCTCGTTGAGCAGGCGGTCGAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAACTACATGATGTCGCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTc                                          <  2:96649/150‑1 (MQ=255)
                                                                                                  cTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGAc                                                                                                                            <  1:223117/63‑1 (MQ=255)
                                                                                                  cTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGAc                                                                                                                            >  2:223117/1‑63 (MQ=255)
                                                                                                   tCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTg                                    <  2:51976/150‑1 (MQ=255)
                                                                                                                  tCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTa                     <  1:174660/150‑1 (MQ=255)
                                                                                                                   ctggCAGGTGTTCAGCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTAc                    <  2:558847/147‑1 (MQ=255)
                                                                                                                   cTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTAc                    <  1:336759/150‑1 (MQ=255)
                                                                                                                      gCAGGTGTTCACCGTTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCa                 >  2:60509/1‑150 (MQ=255)
                                                                                                                      gCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCa                 >  1:92751/1‑150 (MQ=255)
                                                                                                                      gCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCa                 >  1:834509/1‑150 (MQ=255)
                                                                                                                      gCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCCCCCTAGTTACCCTCAAGTCTAGTGAGCACCGGGCATTTTGGCCTGCTTGCGCACGGTGGCCTGTGTGGCGTTACGCa                 >  1:275180/1‑150 (MQ=255)
                                                                                                                       cAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCaa                <  2:781236/150‑1 (MQ=255)
                                                                                                                             gccACCGGTCTGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCtgtgt                            <  1:529305/130‑1 (MQ=255)
                                                                                                                             gccACCGGTCTGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCtgtgt                            >  2:529305/3‑132 (MQ=255)
                                                                                                                               cACCGGTCTGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTATTGCCTGCTTGCTCACGTTGGCCTGTGTGGGGTTAGGCAACGcttctt        >  2:573419/1‑150 (MQ=255)
                                                                                                                               cACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGcttctt        >  1:500373/1‑150 (MQ=255)
                                                                                                                                    gTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTagt                                                                       >  1:546358/1‑82 (MQ=255)
                                                                                                                                    gTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTagt                                                                       <  2:546358/82‑1 (MQ=255)
                                                                                                                                     tCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGCTTCTTCGGCtt  >  1:569752/1‑150 (MQ=255)
                                                                                                                                       |                                                                                                                                                   
TGCGGAAGGAATACCGCCGGGTCGCCCACCGGGTTTTCCAGGTGCGGGTGCTCGGCCAATGGCCGGCGTGTAGGCGCCAGGACACGGGTGATGGTGGGCTCGTTGAGCAGGCGGTCTAGCAGGTGTTCACCGGTCAGACCGGTGGCTCCGGCAAGCAGGACATGCTGAGGCGTCAAGTACATGATGTCTCTCCCTTGTTACACTCAAGTCTAGTGAGTGCCGGGCATTTTTGCCTGCTTGCTCACGTTGGCCTGTGTGGCGTTACGCAACGCTTCTTCGGCTT  >  NC_002947/638610‑638892

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