Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,647,850 T→G 7.4% D430A (GAC→GCC)  PP_4112 ← sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,647,8500TG 7.4% 165.3 / 4.7 54D430A (GAC→GCC) PP_4112sulfate transporter
Reads supporting (aligned to +/- strand):  ref base T (22/28);  new base G (2/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

CTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGAG  >  NC_002947/4647708‑4647992
                                                                                                                                              |                                                                                                                                              
cttcttGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCgg                                                                                                                                          <  4:126866/149‑1 (MQ=255)
   cttGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGc                                                                                                                                         <  1:187886/147‑1 (MQ=255)
     tGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCgtg                                                                                                                                     >  1:253256/1‑149 (MQ=255)
               cGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTgg                                                                                                                            <  5:312315/148‑1 (MQ=255)
                gACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGtt                                                                                                                          >  2:276326/1‑149 (MQ=255)
                  cAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGa                                                                                                                        <  8:38563/149‑1 (MQ=255)
                  cAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCCGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGa                                                                                                                        <  5:105515/149‑1 (MQ=255)
                     aTGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTcccc                                                                                                                                            <  2:281258/126‑1 (MQ=255)
                     aTGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTcccc                                                                                                                                            >  1:281258/1‑126 (MQ=255)
                      tGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGaa                                                                                                                    <  4:311150/149‑1 (MQ=255)
                         accacGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAAc                                                                                                                   >  3:86855/1‑147 (MQ=255)
                          ccacGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAg                                                                                                                <  2:225711/149‑1 (MQ=255)
                          ccacGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAg                                                                                                                <  8:215397/149‑1 (MQ=255)
                                gcAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTc                                                                                                          >  6:112941/1‑149 (MQ=255)
                                 cAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCAGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTc                                                                                                          >  8:102072/1‑148 (MQ=255)
                                    aTGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAgg                                                                                                      <  4:86855/149‑1 (MQ=255)
                                          gCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAg                                                                                                 >  8:173313/1‑148 (MQ=255)
                                           cccGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGa                                                                                                        <  3:183460/140‑1 (MQ=255)
                                           cccGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGa                                                                                                        >  4:183460/1‑140 (MQ=255)
                                            ccGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAg                                                                                                       >  5:183648/1‑140 (MQ=255)
                                            ccGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAg                                                                                                       <  6:183648/140‑1 (MQ=255)
                                                cTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCg                                                                                                         <  6:35943/134‑1 (MQ=255)
                                                cTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCg                                                                                                         >  5:35943/1‑134 (MQ=255)
                                                           gcTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAg                                                                                                                <  2:199240/116‑1 (MQ=255)
                                                           gcTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCAGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGt                                                                               <  7:102072/149‑1 (MQ=255)
                                                                 gcAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATaa                                                                         <  1:276326/149‑1 (MQ=255)
                                                                                 cgcAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCcg                                                         >  1:58482/1‑149 (MQ=255)
                                                                                   cAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCcgcg                                                       >  6:183289/1‑149 (MQ=255)
                                                                                   cAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCcgcg                                                       <  2:253256/149‑1 (MQ=255)
                                                                                     aCGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTc                                                     <  8:237898/149‑1 (MQ=255)
                                                                                        gCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCt                                                   >  1:85587/1‑148 (MQ=255)
                                                                                        gCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTc                                                  <  5:112941/149‑1 (MQ=255)
                                                                                         ccGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTgc                                                                  <  8:185475/132‑1 (MQ=255)
                                                                                              taatCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGgt                                            <  7:173313/149‑1 (MQ=255)
                                                                                                atCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGgtgt                                          <  5:183289/149‑1 (MQ=255)
                                                                                                         gCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATa                                   >  4:45992/1‑147 (MQ=255)
                                                                                                           tCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGc                               >  8:296509/1‑149 (MQ=255)
                                                                                                           tCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGc                               >  5:299070/1‑149 (MQ=255)
                                                                                                                  cTGGCAATCCAGAGTGACCTGGGCCGTGTCCCCGGCCGTGTCGTACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGTGTGTCGGCATACAGCTCCcgc                         <  2:58482/148‑1 (MQ=255)
                                                                                                                    ggCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCa                      >  2:296500/1‑149 (MQ=255)
                                                                                                                     gCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATgg                                              >  4:151556/1‑124 (MQ=255)
                                                                                                                     gCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATgg                                              <  3:151556/124‑1 (MQ=255)
                                                                                                                           cAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCAtgctgc                >  5:25142/1‑148 (MQ=255)
                                                                                                                             gCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAg             >  3:78188/1‑149 (MQ=255)
                                                                                                                               gTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGc           <  6:299070/149‑1 (MQ=255)
                                                                                                                                  gacacggCCGGGGCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCa          >  3:103146/6‑147 (MQ=255)
                                                                                                                                  aCCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCaa         <  7:296509/148‑1 (MQ=255)
                                                                                                                                    cTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGt      >  6:257508/1‑149 (MQ=255)
                                                                                                                                    cTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGt      >  7:234772/1‑149 (MQ=255)
                                                                                                                                     acggCCGGGGACCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCAtgct                  >  6:208538/3‑136 (MQ=255)
                                                                                                                                     acggCCGGGGACCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCAtgct                  <  5:208538/134‑1 (MQ=255)
                                                                                                                                     acggCCGGGGACCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGtt     <  4:103146/147‑1 (MQ=255)
                                                                                                                                       ggCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGa   <  6:25142/149‑1 (MQ=255)
                                                                                                                                        gCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGAg  <  3:45992/149‑1 (MQ=255)
                                                                                                                                              tCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGtt     >  4:252410/1‑140 (MQ=255)
                                                                                                                                              tCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGtt     <  3:252410/140‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
CTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGAG  >  NC_002947/4647708‑4647992

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