Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 257,830 C→A 12.9% R506S (CGC→AGC)  PP_0205 → oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947257,8300CA12.9% 78.4 / 5.1 31R506S (CGC→AGC) PP_0205oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (14/13);  new base A (2/2);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

GGAAATCGGCCTGGCGCGCTTTGCCGAGATCGAGCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGAAAAAAGGCGAGAACGGTGAAATGACCAGCTTCAAGAAGGCCGTCGAACCCTACCTGATCG  >  NC_002947/257688‑257977
                                                                                                                                              |                                                                                                                                                   
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)
                                                                                                                                             ccGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGAAAAAAGGCGAGAACGGTGAAATGACCAGCTTCAAGAAGGCCGTCGAACCCTACCTGATc   <  7:28267/148‑1 (MQ=255)
                                                                                                                                             ccGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGAAAAAAGGCGAGAACGGTGAAATGACCAGCTTCAAGAAGGCCGTCGAACCCTACCTGATCg  >  5:330521/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GGAAATCGGCCTGGCGCGCTTTGCCGAGATCGAGCGTGATCTGGAGCAGATGAAGGCCAACAACCCTCATGAGCTGATGCGTGCGATGGAAGTGTCGGTGATTCGCGATTGCGCCGAGATGGCGGCGCGGTCCTCGCTGTTCCGCGAGGAAAGCCGCTGGGGCCTGTACCACCACCGGGTGGACTTCCCTGAGCGCAACGATGGCGAGTGGTTCTGCCATTGCCACCTGAAAAAAGGCGAGAACGGTGAAATGACCAGCTTCAAGAAGGCCGTCGAACCCTACCTGATCG  >  NC_002947/257688‑257977

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