Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 81,984 G→C 15.5% E98D (GAG→GAC PP_0071 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294781,9840GC15.5% 82.1 / 8.1 33E98D (GAG→GACPP_0071hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/10);  new base C (2/3);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTC  >  NC_002947/81844‑82130
                                                                                                                                            |                                                                                                                                                   
cTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCGTCGACTCGATGCCa                                                                                                                                            <  2:323959/150‑1 (MQ=255)
         cgcAGCCCTTACCAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACGCGATGGCGTCGACTCGATGCCATCgagtatc                                                                                                                                   <  2:273792/150‑7 (MQ=255)
                           ggCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCa                                                                                                                 >  2:359016/1‑150 (MQ=255)
                               aCAAAGCTTTAGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGAATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAgc                                                                                                             <  2:73173/150‑1 (MQ=255)
                                      tttGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCCTCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTa                                                                                                      >  1:66932/1‑150 (MQ=255)
                                      tttGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTa                                                                                                      >  2:209689/1‑150 (MQ=255)
                                                 tGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCCATCGCCGTCACGCCCAGCCTCCAAGGCACAGCGCTGCTATCGAGGCCGACTCAGTACCTGTggg                                                                            >  2:79479/1‑150 (MQ=255)
                                                    gTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCa                                                                                       <  1:359016/150‑1 (MQ=255)
                                                         gaagGCCGTGGCATGTTCTGGTAGAAACTGAGTGCGAATGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTAcc                                                                                  >  2:376815/1‑150 (MQ=255)
                                                              ccGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGc                                                                                                                           <  2:258014/105‑1 (MQ=255)
                                                              ccGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGc                                                                                                                           >  1:258014/1‑105 (MQ=255)
                                                               cGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTggg                                                                            >  1:271752/1‑150 (MQ=255)
                                                                              tCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCcgc                                                             >  2:576801/1‑150 (MQ=255)
                                                                              tCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCcgc                                                             <  1:488846/150‑1 (MQ=255)
                                                                                         tGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAac                                                         >  2:13272/1‑143 (MQ=255)
                                                                                         tGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAac                                                         <  1:13272/143‑1 (MQ=255)
                                                                                                 agaTGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGt                                          <  1:56638/150‑1 (MQ=255)
                                                                                                 agaTGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACCCAGTACCTGTGGGAGCGCGCGTTGCCGCGAACACCGGCAAAGCCGGt                                          >  2:558949/1‑150 (MQ=255)
                                                                                                  gaTGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTACCAAGCCAACGAAAGTACGTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTg                                         >  2:172284/1‑150 (MQ=38)
                                                                                                                tgATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCgg                                                    >  2:473738/1‑125 (MQ=255)
                                                                                                                tgATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCgg                                                    <  1:473738/125‑1 (MQ=255)
                                                                                                                         gtcgaCTCGATGCCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCttcttc                  >  2:436050/5‑150 (MQ=255)
                                                                                                                          gATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACAcc                                                      >  1:8730/1‑113 (MQ=255)
                                                                                                                          gATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACAcc                                                      <  2:8730/113‑1 (MQ=255)
                                                                                                                               tCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGc                                                   >  2:94887/1‑111 (MQ=255)
                                                                                                                               tCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGt            >  2:609274/1‑150 (MQ=255)
                                                                                                                               tCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCACAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGc                                                   <  1:94887/111‑1 (MQ=255)
                                                                                                                                 gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCcatccatc                                <  2:238434/128‑1 (MQ=255)
                                                                                                                                 gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCcatccatc                                >  1:238434/1‑128 (MQ=255)
                                                                                                                                 gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCCGGTaa          >  1:375991/1‑150 (MQ=38)
                                                                                                                                 gATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGCaa          >  2:243803/1‑150 (MQ=255)
                                                                                                                                        tCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCt   >  1:135776/1‑150 (MQ=255)
                                                                                                                                         cGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTc  <  1:376815/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTC  >  NC_002947/81844‑82130

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: