Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 898,346 T→A 16.0% V76D (GTC→GAC)  PP_0781 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947898,3460TA16.0% 60.6 / 5.1 25V76D (GTC→GAC) PP_0781hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (13/8);  new base A (2/2);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCTTGGCGCGGGGTAGGCGCATTCGCCGGCGGCGCCGG  >  NC_002947/898218‑898477
                                                                                                                                |                                                                                                                                   
ccAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTc                                                                                                                              <  1:433665/136‑1 (MQ=255)
ccAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTc                                                                                                                              <  1:2969/136‑1 (MQ=255)
  aGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGa                                                                                                                            >  1:542003/1‑136 (MQ=255)
  aGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCAAGCGCTTCGTCGACCTCGa                                                                                                                            >  2:491813/1‑136 (MQ=255)
            cAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGc                                                                                                                  >  1:297893/1‑136 (MQ=255)
                 tCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCAtcgt                                                                                                             >  1:573113/1‑136 (MQ=255)
                   aaTGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCAtcgtcg                                                                                                           <  2:542003/136‑1 (MQ=255)
                                ccTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAg                                                                                              >  1:18753/1‑136 (MQ=255)
                                  tCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAggg                                                                                            >  1:619957/1‑136 (MQ=255)
                                   cGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGt                                                                                           <  1:111328/136‑1 (MQ=255)
                                       gggCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTg                                                                                       >  2:450807/1‑136 (MQ=255)
                                                           tCGACGCTGTGCGCGCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGAtcgt                                                                   <  1:453153/136‑1 (MQ=255)
                                                                       gcTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCAc                                                       >  1:593194/1‑136 (MQ=255)
                                                                            gCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGACGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTa                                                  >  1:244983/1‑136 (MQ=255)
                                                                            gCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGACGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTa                                                  >  1:241885/1‑136 (MQ=255)
                                                                              cTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTAcc                                                >  2:64696/1‑136 (MQ=255)
                                                                                 cTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTcgac                                                                                                         <  1:25950/76‑1 (MQ=255)
                                                                                 cTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTcgac                                                                                                         >  2:25950/1‑76 (MQ=255)
                                                                                          tGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCtt                                    >  1:224935/1‑136 (MQ=255)
                                                                                           gAATGGCAGACGGCGCACATGCTACGCAAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCTTg                                   <  1:491813/136‑1 (MQ=255)
                                                                                             aTGGCAGACGGCGCACATGCTACGCCAGCGCTTCGACGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCTTGgc                                 <  2:244983/136‑1 (MQ=255)
                                                                                             aTGGCAGACGGCGCACATGCTACGCCAGCGCTTCGACGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCTTGgc                                 <  2:241885/136‑1 (MQ=255)
                                                                                                                       aGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCTTGGCGCGGGGTAGGCGCATTCGCCggcggc       <  2:18753/136‑1 (MQ=255)
                                                                                                                            ttCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCTTGGCGCGGGGTAGGCGCATTCGCCGGCGGCGCCgg  >  2:327541/1‑136 (MQ=255)
                                                                                                                            ttCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCTTGGCGCGGGGTAGGCGCATTCGCCGGCGGCGCCgg  >  1:369951/1‑136 (MQ=255)
                                                                                                                                |                                                                                                                                   
CCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCTTGGCGCGGGGTAGGCGCATTCGCCGGCGGCGCCGG  >  NC_002947/898218‑898477

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