Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,419,194 T→G 10.0% Y51S (TAC→TCC)  PP_5553 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,419,1940TG10.0% 99.4 / 4.3 40Y51S (TAC→TCC) PP_5553hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/22);  new base G (2/2);  total (16/24)
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.62e-01

TCTGCAGCGTCGCCTTGGCGCGCGCGCCGTTCGATTCGTCGCGAAGGATGTGGGTGAATTGCCCTGTCAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGCTGCTGGATGGACATGTAGTGA  >  NC_002947/3419052‑3419330
                                                                                                                                               |                                                                                                                                        
tCTGCAGCGTCGCCTTGGCGCGCGCGCCGTTCGATTCGTCGCGAAGGATGTGGGTGAATTGCCCTGTCAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAg                                                                                                                                    <  2:135221/149‑1 (MQ=255)
    cAGCGTCGCCTTGGCGCGCGCGCCGTTCGATTCGTCGCGAAGGATGTGGGTGAATTGCCCTGTCAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATACGATGGAgg                                                                                                                                   <  1:185726/146‑1 (MQ=255)
      gCGTCGCCTTGGCGCGCGCGCCGTTCGATTCGTCGCGAAGGATGTGGGTGAATTGCCCTGGCAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCATTGCGTAGCTCATGCCCGCGTGCCCTCCATACGATGGAgggcac                                                                                                                               <  1:15954/149‑5 (MQ=255)
       cgagctcTTGGCGCGGGCGCCGTGCGATTCGTCGCGAAGGATGTGGGGGAATTGCCCTGGCAGGTTGGGCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCGCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGc                                                                                                                              <  5:190897/142‑1 (MQ=255)
        gTCGCCTTGGCGCGCGCGCCGTTCGATTCGTCGCGAAGGATGTGGGTGAATTGCCCTGTCAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAg                                                                                                                            <  7:236164/149‑1 (MQ=255)
         tCGCCTTGGCGCGCGCGCCGTTCGATTCGTCGCGAAGGATGTGGGTGAATTGCCCTGTCAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAgg                                                                                                                           <  7:57646/149‑1 (MQ=255)
              ttGGCGCGCGCGCCGTTCGATTCGTCGCGAAGGATGTGGGTGCATTGCCCTGTCAGGGTGGTCTGCGCCTTAAGGCAGGCAAACGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCCTCGGAGGGAGGCGTGCacgagcg                                                                                                                       >  6:190897/1‑142 (MQ=255)
                   gcgcgcgcCGTTCGATTCGTCGCGAAGGATGTGGGTGAATTGCCCTGTCAGGTTGGTCTGCGCCTTTAGGCATGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCggg                                                                                                                 <  2:216953/149‑1 (MQ=255)
                                     gtcgCGAAGGATGTGGGTGAATTGCCCTGTCAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGttt                                                                                               >  1:164439/1‑149 (MQ=255)
                                     gtcgCGAAGGATGTGGGTGAATTGCCCTGTCAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGttt                                                                                               <  2:34/149‑1 (MQ=255)
                                         cgAAGGATGTGGGTGAATTGCCCTGTCAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTa                                                                                           <  8:99754/149‑1 (MQ=255)
                                                   gggTGAATTGCCCTGTCAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATg                                                                                 >  7:113239/1‑149 (MQ=255)
                                                                tGTCAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGATCATG‑CCGCGTGACCTCCAACGTATGGAGGCGTGCAGGTTGAAGCGGTGGACGACGTAGCGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTg                                                                    >  8:218345/1‑149 (MQ=255)
                                                                   cAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCtt                                                                 <  4:8372/149‑1 (MQ=255)
                                                                    aGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTg                                                                <  4:118959/149‑1 (MQ=255)
                                                                       ttGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGc                                                              >  5:89736/1‑148 (MQ=255)
                                                                       ttGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGc                                                              >  7:32418/1‑148 (MQ=255)
                                                                               cgcCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCg                                                     >  6:173717/1‑149 (MQ=255)
                                                                                 ccTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGa                                                   <  3:222995/149‑1 (MQ=255)
                                                                                      aggctggcAAAAGCGTTGTGTCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTAGATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGGAGCCGTGTTGCTTGGCCAGTCGCCGGACttt                                               <  7:218345/148‑1 (MQ=255)
                                                                                         caggcaAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAg                                           <  8:113239/149‑1 (MQ=255)
                                                                                           ggcaAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGt                                         >  6:231729/1‑149 (MQ=255)
                                                                                            gcaAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGtt                                        <  2:47720/149‑1 (MQ=255)
                                                                                               aaaGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAgg                                     >  4:206421/1‑149 (MQ=255)
                                                                                                              aGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGgct                      <  4:259745/149‑1 (MQ=255)
                                                                                                                   gTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGc                                                              >  7:200782/1‑104 (MQ=255)
                                                                                                                   gTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGc                                                              <  8:200782/104‑1 (MQ=255)
                                                                                                                   gTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGctgctg                  >  4:41534/1‑148 (MQ=255)
                                                                                                                    tAGCTCATG‑CCGCGTGCCCTCCATCGTTTGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGCTGCtgga                <  3:31203/149‑1 (MQ=255)
                                                                                                                      gCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGCTGCtggat               <  2:164439/148‑1 (MQ=255)
                                                                                                                          aTG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGCTGCTGGATGGAc           <  6:89736/148‑1 (MQ=255)
                                                                                                                            g‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGCTGCTGGATGGACATg        >  1:114689/1‑149 (MQ=255)
                                                                                                                              ccGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGCTGCTGGATGGACATGt       <  5:173717/149‑1 (MQ=255)
                                                                                                                                   acgccTCCATACGATGGAGGCGTGCAGGTTGATGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGCTGCTGGATGGACATGTAGTGa  >  5:50078/4‑149 (MQ=255)
                                                                                                                                     cccTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGCTGCTGGATGGACATGTAGt    <  1:130456/145‑1 (MQ=255)
                                                                                                                                     cccTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGCTGCTGGATGGACATGTAGt    >  2:130456/1‑145 (MQ=255)
                                                                                                                                       cTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGCTGCTGGATGGACATGTAGTGa  >  2:18016/1‑145 (MQ=255)
                                                                                                                                       cTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGCTGCTGGATGGACATGTAGTGa  <  1:18016/145‑1 (MQ=255)
                                                                                                                                        tCCATCGTATGGAGGCGTGCAGGTTGAAGC‑GGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAg                                           <  5:171758/102‑1 (MQ=255)
                                                                                                                                        tCCATCGTATGGAGGCGTGCAGGTTGAAG‑AGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAg                                           >  6:171758/1‑102 (MQ=255)
                                                                                                                                               |                                                                                                                                        
TCTGCAGCGTCGCCTTGGCGCGCGCGCCGTTCGATTCGTCGCGAAGGATGTGGGTGAATTGCCCTGTCAGGTTGGTCTGCGCCTTAAGGCAGGCAAAAGCGTTGTGGCTCAGTGCGTAGCTCATG‑CCGCGTGCCCTCCATCGTATGGAGGCGTGCAGGTTGAAGCGGGGGACGACGTAGGGCGTTTGGTAGTGACGAATGCGCAGCCGTGTTGCTTGGCCAGTCGCCGGACTTTGAAGGTTAGGTCGATGGTTGCGGCTGCTGGATGGACATGTAGTGA  >  NC_002947/3419052‑3419330

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