Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,562,063 G→T 7.0% P157T (CCG→ACG)  PP_3145 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,562,0630GT 7.0% 180.2 / 3.8 57P157T (CCG→ACG) PP_3145hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/31);  new base T (2/2);  total (24/33)
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.46e-01

TCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATA  >  NC_002947/3561937‑3562208
                                                                                                                               |                                                                                                                                                 
tCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCgcgg                                                                                                                             >  7:5055/1‑149 (MQ=255)
  acacCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaagg                                                                                                                                  <  6:2842/142‑4 (MQ=255)
  acacCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaagg                                                                                                                                  >  5:2842/1‑139 (MQ=255)
         ccGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGcc                                                                                                                    >  3:93573/1‑149 (MQ=255)
            cgATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCg                                                                                                                 >  2:159385/1‑149 (MQ=255)
               tAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGtt                                                                                                               >  7:131758/1‑148 (MQ=255)
                 aGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                            <  7:20793/149‑1 (MQ=255)
                 aGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                            <  6:162985/149‑1 (MQ=255)
                 aGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGa                                                                                                             <  8:132144/148‑1 (MQ=255)
                  gcgcGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAgg                                                                                                                       >  7:40139/1‑137 (MQ=255)
                  gcgcGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAgg                                                                                                                       <  8:40139/137‑1 (MQ=255)
                      gATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGC                                                                                                       <  1:159385‑M1/149‑1 (MQ=255)
                      gATAGGGATCACCTCGTAAGGCCGGAGGCCCTCGTTACTGGGGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGC                                                                                                       <  7:143815‑M1/149‑1 (MQ=255)
                           ggATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaaggg                                                                                                                                 <  8:170960/118‑5 (MQ=255)
                           ggATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaaggg                                                                                                                                 >  7:170960/1‑114 (MQ=255)
                            gATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttat                                                                                                 >  5:1028‑M1/1‑145 (MQ=255)
                             aTCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGAt                                                                                                <  4:93573/149‑1 (MQ=255)
                                aCCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCCCCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTc                                                                                             >  2:318708/1‑149 (MQ=255)
                                 ccTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGt                                                                                                   >  3:271667/1‑142 (MQ=255)
                                 ccTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGt                                                                                                   <  4:271667/142‑1 (MQ=255)
                                  cTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCt                                                                                           >  1:216897/1‑149 (MQ=255)
                                  cTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCt                                                                                           <  7:11254/149‑1 (MQ=255)
                                      tAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCGACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTg                                                                                          >  1:241985/1‑147 (MQ=255)
                                       aaGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGc                                                                                      >  6:285762/1‑149 (MQ=255)
                                        aGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCGCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGc                                                                                      <  5:285762/148‑1 (MQ=255)
                                          gCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                            <  7:224920/124‑1 (MQ=255)
                                          gCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                            >  8:224920/1‑124 (MQ=255)
                                                 gCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCtgctgc                                                                             >  7:28615/1‑148 (MQ=255)
                                                    cTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCG                                                                                                      >  4:252782‑M1/1‑120 (MQ=255)
                                                    cTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCG                                                                                                      <  3:252782‑M1/120‑1 (MQ=255)
                                                       gTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATa                                                                       <  8:5055/148‑1 (MQ=255)
                                                          aCTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAAccg                                                                   <  5:247132‑M2/149‑4 (MQ=255)
                                                            tggtggACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGAGAAccgcg                                                                 <  8:80284‑M2/149‑6 (MQ=255)
                                                              gtggACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGt                                                               <  3:287060/149‑1 (MQ=255)
                                                               tggACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGtt                                                              <  2:241985/149‑1 (MQ=255)
                                                                 gACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTgg                                                            <  3:153293/149‑1 (MQ=255)
                                                                  acacTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttatctcgcagcagacgctgcaggaaatcacccgctacttca                                                           <  6:1028‑M1/149‑43 (MQ=255)
                                                                      tGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttatctcgcagcagacgctgcaggaaatcac                                                                      <  2:92087‑M1/134‑32 (MQ=255)
                                                                      tGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttatctcgcagcagacgctgcaggaaatcac                                                                      >  1:92087‑M1/1‑103 (MQ=255)
                                                                        ccGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCgg                                                      <  1:309204/148‑1 (MQ=255)
                                                                           tGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTg                                                  >  2:228375/1‑149 (MQ=255)
                                                                                gCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTa                                             <  8:28615/149‑1 (MQ=255)
                                                                                 ccGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTgggg                                               >  7:183866/1‑146 (MQ=255)
                                                                                 ccGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTgggg                                               >  7:43944/1‑146 (MQ=255)
                                                                                   gggCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGaa                                          >  6:253517/1‑149 (MQ=255)
                                                                                    ggCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAAc                                         <  2:216897/149‑1 (MQ=255)
                                                                                       gcACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTg                                                  >  8:101645/1‑137 (MQ=255)
                                                                                       gcACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTg                                                  <  7:101645/137‑1 (MQ=255)
                                                                                                       ccGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGaa                        >  5:6853/1‑148 (MQ=255)
                                                                                                       ccGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAac                       <  8:131758/149‑1 (MQ=255)
                                                                                                         gggCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAacac                     <  8:183866/149‑1 (MQ=255)
                                                                                                             cAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTg                 <  1:228375/149‑1 (MQ=255)
                                                                                                                  cccTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTa            >  8:193783/1‑149 (MQ=255)
                                                                                                                   ccTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTAc           <  6:247132/149‑1 (MQ=255)
                                                                                                                         ccAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTg     <  8:43944/149‑1 (MQ=255)
                                                                                                                            aCGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGata  <  1:318708/149‑1 (MQ=255)
                                                                                                                            aCGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGata  <  6:6853/149‑1 (MQ=255)
                                                                                                                               |                                                                                                                                                 
TCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATA  >  NC_002947/3561937‑3562208

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