Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 257,825 T→G 12.5% L504R (CTG→CGG)  PP_0205 → oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947257,8250TG12.5% 81.0 / 5.3 32L504R (CTG→CGG) PP_0205oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (14/14);  new base G (2/2);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

AGAAGATGGAAATCGGCCTGGCGCGCTTTGCCGAGATCGAGCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGAAAAAAGGCGAGAACGGTGAAATGACCAGCTTCAAGAAGGCCGTCGAACCCTACC  >  NC_002947/257681‑257971
                                                                                                                                                |                                                                                                                                                  
agaagaTGGAAATCGGCCTGGCGCGCTTTGCCGAGATCGAGCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTc                                                                                                                                                <  3:156120/149‑1 (MQ=255)
       ggAAATCGGCCTGGCGCGCTTTGCCGAGATCGAGCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAg                                                                                                                                          >  5:1544/1‑148 (MQ=255)
         aaaTCGGCCTGGCGCGCTTTGCCGAGATCGAGCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGa                                                                                                                                        >  8:274092/1‑148 (MQ=255)
         aaaTCGGCCTGGCGCGCTTTGCCGAGATCGAGCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCGGAACAGCGAGGac                                                                                                                                       <  1:148604/149‑2 (MQ=255)
                 cTGGCGCGCTTTGCCGAGATCGAGCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTg                                                                                                                               <  8:149342/149‑1 (MQ=255)
                                gagaTCGAGCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTaccacca                                                                                                                >  6:289881/1‑149 (MQ=255)
                                      gAGCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGt                                                                                                          >  6:100038/1‑149 (MQ=255)
                                       aGCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCACCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTg                                                                                                         <  1:163311/149‑1 (MQ=255)
                                        gCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTgg                                                                                                        <  7:274092/149‑1 (MQ=255)
                                                  gAGCAGATGAATGCCAACAACCCTCATGAGCTTCTGCGTGCGATGGAAGTGTCGGTGAGTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTccc                                                                                                 <  8:287/146‑1 (MQ=255)
                                                  gAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCt                                                                                                <  2:2147/147‑1 (MQ=255)
                                                    gCAGATCAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAgc                                                                                            >  5:321387/1‑149 (MQ=255)
                                                        aTGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCa                                                                                         >  8:28267/1‑148 (MQ=255)
                                                                  acaacCCTCATGAGCTGATGCGTGCGATGGTAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAg                                                                              >  1:31705/1‑149 (MQ=255)
                                                                    aacCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTg                                                                            <  5:289881/149‑1 (MQ=255)
                                                                        cTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTc                                                                        >  7:13991/1‑149 (MQ=255)
                                                                        cTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCGGAACAGCGAGGac                                                                                                                                       >  5:150092/1‑85 (MQ=255)
                                                                        cTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCGGAACAGCGAGGac                                                                                                                                       <  6:150092/86‑2 (MQ=255)
                                                                         tCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCt                                                                       >  1:25684/1‑149 (MQ=255)
                                                                                cTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCAtt                                                                 <  6:1544/148‑1 (MQ=255)
                                                                                     gcgtgcgATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCAc                                                            <  6:321387/148‑1 (MQ=255)
                                                                                         gcgATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGa                                                        >  6:60906/1‑148 (MQ=255)
                                                                                             tGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGaaaa                                                     >  2:291145/1‑147 (MQ=255)
                                                                                                     tCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAACGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCCCCTGCAAAAAGGCGAGaa                                           >  6:337232/1‑149 (MQ=255)
                                                                                                       ggTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGAAAAAAGGCGAGAACg                                         >  6:131921/1‑149 (MQ=255)
                                                                                                              cgcgATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGAAAAAAGGCGAGAACGGTGaaa                                   <  5:100038/148‑1 (MQ=255)
                                                                                                               gcgATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGtt                                                                         <  7:276480/109‑1 (MQ=255)
                                                                                                               gcgATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGtt                                                                         >  8:276480/1‑109 (MQ=255)
                                                                                                                 gATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGAAAAAAGGCGAGAACGGTGAAATGAc                               <  2:31705/149‑1 (MQ=255)
                                                                                                                        ccGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGAAAAAAGGCGAGAACGGTGAAATGACCAGCTTc                        <  5:131921/149‑1 (MQ=255)
                                                                                                                                         tCCTCGCGGAACAGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGAAAAAAGGCGAGAACGGTg                                      >  4:321871/1‑118 (MQ=255)
                                                                                                                                              gCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGAAAAAAGGCGAGAACGGTGAAATGACCAGCTTCAAGAAGGCCGTCGAACCCTAcc  <  5:60906/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AGAAGATGGAAATCGGCCTGGCGCGCTTTGCCGAGATCGAGCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGAAAAAAGGCGAGAACGGTGAAATGACCAGCTTCAAGAAGGCCGTCGAACCCTACC  >  NC_002947/257681‑257971

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: ‑