Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,871,819 C→T 11.8% S185L (TCG→TTG)  PP_2524 → HD domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,871,8190CT11.8% 83.2 / 3.9 34S185L (TCG→TTG) PP_2524HD domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (15/15);  new base T (2/2);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

ATTGCCCGGTGTTTTTACGTCGCCGGGCGAATGGGGTCGGCGCTCTATGACTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACGAGATTACCGGGGCAGGTTAGTCATGCCGGGCTTGTACCGTCGTGCTTGAGC  >  NC_002947/2871680‑2871953
                                                                                                                                           |                                                                                                                                      
aTTGCCCGGTGTTTTTACGTCGCCGGGCGAATGGGGTCGGCGCTCTATGACTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTaaa                                                                                                                              >  1:42741/1‑150 (MQ=255)
 ttGCCCGGTGTTTTTACGTCGCCGGGCGAATGGGGTCGGCGCTCTATGACTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTaaaa                                                                                                                             >  2:610556/1‑150 (MQ=255)
      cGGTGTTTTTACGTCGCCGGGCGAATGGGGTCGGCGCTCTATGACTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCaa                                                                                                                        >  1:455906/1‑150 (MQ=255)
       ggTGTTTTTACGTCGCCGGGCGAATGGGGTCGGCGCTCTATGACTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACGAGGACCGGAACTATTGCAGTGGGCGCTTCTAGACCAAACTGCTCAAGCTCGCCTTGGGCTTTAAAACCCaag                                                                                                                       >  1:317034/1‑150 (MQ=255)
                  gTCGCCGGGCGAATGGGGTCGGCGCTCTATGACTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGgcgcgc                                                                                                            >  2:420032/1‑150 (MQ=255)
                             aaTGGGGTCGGCGCTCTATGACTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGa                                                                                                 <  2:455906/150‑1 (MQ=255)
                                      ggCGCTCTATGACTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATAc                                                                                        >  1:187613/1‑150 (MQ=255)
                                         gcTCTATGACTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTTGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCg                                                                                     >  1:219990/1‑150 (MQ=255)
                                              aTGACTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAggg                                                                                                                  >  2:115469/1‑116 (MQ=255)
                                              aTGACTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAggg                                                                                                                  <  1:115469/116‑1 (MQ=255)
                                                 aCTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGgc                                                                                                                >  1:367176/1‑115 (MQ=255)
                                                 aCTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGgc                                                                                                                <  2:367176/115‑1 (MQ=255)
                                                 aCTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGaaa                                                                             <  1:610556/150‑1 (MQ=255)
                                                 aCTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGaaa                                                                             <  1:420032/150‑1 (MQ=255)
                                                    tCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTTGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCCGCTTGAAACCt                                                                          <  2:317034/150‑1 (MQ=255)
                                                     cGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCtt                                                                         <  2:257380/150‑1 (MQ=255)
                                                        gaATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCt                                                                      <  2:780227/150‑1 (MQ=255)
                                                           tCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGCCCCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCttcc                                                                   <  2:760965/150‑3 (MQ=255)
                                                            ccTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCCAACTGCTCAAGCTCGCCTCGGGCTTTaaaa                                                                                                                             <  1:23142/91‑1 (MQ=255)
                                                            ccTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTaaaa                                                                                                                             >  2:23142/1‑91 (MQ=255)
                                                              tttCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCa                                                                >  2:124704/1‑150 (MQ=255)
                                                                            gggCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACgaga                                                  <  1:124704/150‑1 (MQ=255)
                                                                                   ccAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACGAGATTACCgg                                           <  2:42741/150‑1 (MQ=255)
                                                                                                          aGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACGAGATTACCGGGGCAGGTTAGTCATGCCGGGCtt                    <  1:411927/150‑1 (MQ=255)
                                                                                                               ttCCAGACCAAACTGCTCAAGCTCGCCTTGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACGAGATTACCGGGGCAGGTTAGTCATGCCGGGCTTGTAcc               <  2:219990/150‑1 (MQ=255)
                                                                                                                    gACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGAc                                                      <  2:36857/106‑1 (MQ=255)
                                                                                                                    gACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGAc                                                      >  1:36857/1‑106 (MQ=255)
                                                                                                                    gACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACGAGATTACCGGGGCAGGTTAGTCATGCCGGGCTTGTACcgtcgt          >  1:376452/1‑150 (MQ=255)
                                                                                                                       cAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACGAGATTACCGGGGCAGGTTAGTCATGCCGGGCTTGTACCGTCGTGCt       >  1:76169/1‑150 (MQ=255)
                                                                                                                            tGCTCAAGCTCGCCTCGGGCGTTAAAACTCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACGAGATTACCGGGGCAGGTTAGTCATGCCGGGCTTGGACCGTCGTGCTTGAGc  >  1:285950/1‑150 (MQ=255)
                                                                                                                             gCTCAAGCTAGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACGAGATTACCGGGGCAGGTTAGTCATGCCGGGCTTg                   >  2:36316/1‑132 (MQ=255)
                                                                                                                             gCTCAAGCTAGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACGAGATTACCGGGGCAGGTTAGTCATGCCGGGCTTg                   <  1:36316/132‑1 (MQ=255)
                                                                                                                                     tCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACga                                                    <  2:839257/91‑1 (MQ=255)
                                                                                                                                     tCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACga                                                    >  1:839257/1‑91 (MQ=255)
                                                                                                                                           |                                                                                                                                      
ATTGCCCGGTGTTTTTACGTCGCCGGGCGAATGGGGTCGGCGCTCTATGACTTCGAGAATCCTTTCGCGCAGCAGAGGGCGTACCAGGACGGGGCGTATGCCATCGAGCACTTCCAGACCAAACTGCTCAAGCTCGCCTCGGGCTTTAAAACCCAAGAAGGGGCGCGCCTTGCTGCCGAGCGTCATACCCGGCTTGAAACCTTTCTTACCCAGTTCATGGACGAGATTACCGGGGCAGGTTAGTCATGCCGGGCTTGTACCGTCGTGCTTGAGC  >  NC_002947/2871680‑2871953

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: