Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,930,569 C→A 21.0% D42Y (GAC→TAC)  PP_1725 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,930,5690CA21.0% 43.2 / 2.6 19D42Y (GAC→TAC) PP_1725hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (9/6);  new base A (2/2);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

GTTGGTTGGACAGGCGGTTGATGTTGTTGTGCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTAGGCGTGCAGGTGGCCCATCGCATGGTGGGCGACCTGATAGTTGAGGCCGTCGAGCAGGACCAGGATGACGTTGTGTTGCATGTTTTCTCCAGTCATGAACCGCGTTGGCACCTTCGCGGG  >  NC_002947/1930434‑1930701
                                                                                                                                       |                                                                                                                                    
gttggttgGCCAGGCGGTTGAGGTTGTTGTGCACGATCCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTc                                                                                                                                      <  1:443980/136‑1 (MQ=255)
gttggttgGACAGGCGGTTGATGTTGTTGTGCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTc                                                                                                                                      >  2:277543/1‑136 (MQ=255)
       ggACAGGCGGTTGATGTTGTTGTGCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTAGGCCTCg                                                                                                                               >  1:185740/1‑136 (MQ=255)
       ggACAGGCGGTTGATGTTGTTGTGCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTAGGCCTCg                                                                                                                               >  1:185930/1‑136 (MQ=255)
        gACAGGCGGTTGATGTTGTTGTGCACGATGCCGCTGTCGATAGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCAGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGa                                                                                                                              <  2:108677/136‑1 (MQ=255)
            ggCGGTTGATGTTGTTGTGCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTa                                                                                                                          >  1:284256/1‑136 (MQ=255)
                gTTGATGTTGTTGTGCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTAGGCg                                                                                                                      >  2:114934/1‑136 (MQ=255)
                  tGATGTTGTTGTGCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTAGGCCTCGACGTAGGCGTg                                                                                                                    <  2:185740/136‑1 (MQ=255)
                  tGATGTTGTTGTGCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTAGGCCTCGACGTAGGCGTg                                                                                                                    <  2:185930/136‑1 (MQ=255)
                       ttgttgtGCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTAGGCGTGCAGGt                                                                                                               <  2:474567/136‑1 (MQ=255)
                            gtgCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTAGGCGTGCAGGTGGccc                                                                                                          >  1:164918/1‑136 (MQ=255)
                              gCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTAGGCGTGCAGGTGGCCCAt                                                                                                        >  1:301195/1‑136 (MQ=255)
                              gCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTAGGCGTGCAGGTGGCCCAt                                                                                                        >  1:80145/1‑136 (MQ=255)
                                       cGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTAGGCGTGCAGGTGGCCCATCGCAtggtg                                                                                               >  1:11955/1‑136 (MQ=255)
                                                               tGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTAGGCGTGCAGGTGGCCCATCGCATGGTGGGCGACCTGATAGTTGAGGCCGTc                                                                       <  1:314986/136‑1 (MQ=255)
                                                                     tGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACATAGGCGTGCAGGTGGCCCATCGCATGGTGGGCGACCTGATAGTTGAGGCCGTCGAGCAg                                                                 <  1:223692/136‑1 (MQ=255)
                                                                          tCGTACAGCGGCCGCGACAGCGATGGCAGGTCGCATTCCACGCGGTACTGCGCGGCGCGGTCGGCCTCGACGTAGGCGTGCAGGTGGCCCATCGCATGGTGGGCGACCTGATAGTTGAGGCCGTCGAGCAGGACCa                                                            <  2:365603/136‑1 (MQ=255)
                                                                          tCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTAGGCGTGCAGGTGGCCCATCGCATGGTGGGCGACCTGATAGTTGAGGCCGTCGAGCAGGACCa                                                            <  2:62186/136‑1 (MQ=255)
                                                                                     ccGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTAGGCGTGCAGGTGGCCCATCGCATGGTGGGCGACCTGATAGTTGAGGCCGTCGAGCAGGACCAGGATGACGTtg                                                 >  2:192825/1‑136 (MQ=255)
                                                                                                          gCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTAGGCGTGCAGGTGGCCCATCGCATGGTGGGCGACCTGATAGTTGAGGCCGTCGAGCAGGACCAGGATGACGTTGTGTTGCATGTTTTCTCCAGTc                            >  1:334309/1‑136 (MQ=255)
                                                                                                                                    ggTCGGCCTCGACGTAGGCGTGCAGGTGGCCCATCGCATGGTGGGCGACCTGATAGTTGAGGCCGTCGAGCAGGACCAGGATGACGTTGTGTTGCATGTTTTCTCCAGTCATGAACCGCGTTGGCACCTTCGCggg  >  1:503641/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                    
GTTGGTTGGACAGGCGGTTGATGTTGTTGTGCACGATGCCGCTGTCGATCGGTGGCACGCCGGTGAGGATGCATTCGTACAGCGGCCGCGACAGCGATGGCAGCTCGCATTCCACGCGGTACAGCGCGGCGCGGTCGGCCTCGACGTAGGCGTGCAGGTGGCCCATCGCATGGTGGGCGACCTGATAGTTGAGGCCGTCGAGCAGGACCAGGATGACGTTGTGTTGCATGTTTTCTCCAGTCATGAACCGCGTTGGCACCTTCGCGGG  >  NC_002947/1930434‑1930701

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