Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 900,535 C→G 6.8% intergenic (+27/+18) PP_0784 → / ← cysZ hypothetical protein/sulfate transport protein CysZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947900,5350CG 6.8% 172.0 / 2.3 55intergenic (+27/+18)PP_0784/cysZhypothetical protein/sulfate transport protein CysZ
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base G (2/2);  total (28/27)
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.61e-01

CAGCCTTGGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGCCAGCGC  >  NC_002947/900387‑900679
                                                                                                                                                    |                                                                                                                                                
cAGCCTTGGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGccc                                                                                                                                                  <  4:227105/149‑1 (MQ=255)
cAGCCTTGGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGccc                                                                                                                                                  >  6:402322/1‑149 (MQ=255)
    cTTGGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCa                                                                                                                                              <  5:402322/149‑1 (MQ=255)
    cTTCGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCTGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAATCAATCCGCAAAAAGGGGCTGCCCTGCa                                                                                                                                              <  8:103496/149‑1 (MQ=255)
               gcggGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTc                                                                                                                                   >  5:67249/1‑149 (MQ=255)
                   gTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTc                                                                                                                               <  4:136792/149‑1 (MQ=255)
                    tCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCa                                                                                                                              <  1:163812/149‑1 (MQ=255)
                               ttcTAACCTCTAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCATTCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAgctccgctcc                                                                                                                    <  2:132431/146‑1 (MQ=255)
                                             aGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGaacaac                                                                                                     <  7:387569/149‑1 (MQ=255)
                                              ggTTAAGGCAGTGGAACATCTAATACTGATTTCTAAGATTGCCACCCTCGAACCCTACAGTTTATCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACg                                                                                                    >  3:60774/1‑149 (MQ=255)
                                               gTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCAcc                                                                                                              <  3:164110/138‑1 (MQ=255)
                                               gTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCAcc                                                                                                              >  4:164110/1‑138 (MQ=255)
                                               gTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGt                                                                                                   <  6:408455/149‑1 (MQ=255)
                                                   cggCTGTGTAACATAAAATACGGTTTTCTCATATTCCCCCCCTCCAACCCAACTGTTTCCCCTAATTAAAAAAACCAAGCCCAAAACCGGGGCGCCCGTCCCCCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGc                                                                                                <  6:67249/147‑1 (MQ=255)
                                                    ggCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTccc                                                                                                                   >  2:191024/1‑128 (MQ=255)
                                                    ggCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTccc                                                                                                                   <  1:191024/128‑1 (MQ=255)
                                                    ggCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGccc                                                                                              >  8:298477/1‑149 (MQ=255)
                                                           gAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCAGGGCAGcccc                                                                                                                                         >  3:263762/1‑99 (MQ=255)
                                                           gAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCAGGGCAGcccc                                                                                                                                         <  4:263762/99‑1 (MQ=255)
                                                            aaCATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCATACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGACCAACGTAGCACAGGCCACt                                                                                      >  6:290957/1‑149 (MQ=255)
                                                                       cTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGa                                                                                                          >  3:388486/1‑118 (MQ=255)
                                                                       cTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGa                                                                                                          <  4:388486/118‑1 (MQ=255)
                                                                           tttCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGcatcat                                                                        >  8:250917/1‑148 (MQ=255)
                                                                            ttCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGcatcat                                                                        >  1:283348/1‑147 (MQ=255)
                                                                               tcAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCa                                                                   >  4:7938/1‑149 (MQ=255)
                                                                                 aGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGc                                                                    <  2:346810/146‑1 (MQ=255)
                                                                                 aGATTCCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACg                                                                 >  8:38897/1‑149 (MQ=255)
                                                                                  gATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGcccc                                                                                                                                         <  6:24133/76‑1 (MQ=255)
                                                                                  gATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGcccc                                                                                                                                         >  5:24133/1‑76 (MQ=255)
                                                                                       ccACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGcatc                                                                          <  2:87190/134‑1 (MQ=255)
                                                                                       ccACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGcatc                                                                          >  1:87190/1‑134 (MQ=255)
                                                                                       ccACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCAc                                                            <  3:7938/148‑1 (MQ=255)
                                                                                            cTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTgg                                                      <  7:250917/149‑1 (MQ=255)
                                                                                               gAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCTCACCCAGAACAACGTAGCCGCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGaa                                                    <  5:290957/148‑1 (MQ=255)
                                                                                                    cTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAAtcatca                                              <  2:283348/149‑1 (MQ=255)
                                                                                                             ctcCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCagg                                      <  2:243519/148‑1 (MQ=255)
                                                                                                              tcCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCagg                                      >  1:253501/1‑147 (MQ=255)
                                                                                                               ccGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTGCATCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCaggtag                                   >  1:108190/1‑149 (MQ=255)
                                                                                                                cGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCaggtagg                                  <  7:298477/149‑1 (MQ=255)
                                                                                                                  aCTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCaggtaggt                                 >  5:195585/1‑148 (MQ=255)
                                                                                                                      cTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATgc                            >  8:384037/1‑149 (MQ=255)
                                                                                                                       tAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATgcc                           >  3:27241/1‑149 (MQ=255)
                                                                                                                         aaaaGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACt                                                                                      >  5:330102/1‑88 (MQ=255)
                                                                                                                         aaaaGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACt                                                                                      <  6:330102/88‑1 (MQ=255)
                                                                                                                            aGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCaa                      <  4:60774/149‑1 (MQ=255)
                                                                                                                               aaTCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATgc                            >  4:394903/1‑140 (MQ=255)
                                                                                                                               aaTCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATgc                            <  3:394903/140‑1 (MQ=255)
                                                                                                                                 tCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGccc                 >  8:58753/1‑149 (MQ=255)
                                                                                                                                    gCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAgg              >  8:180353/1‑149 (MQ=255)
                                                                                                                                    gCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAgg              >  5:283376/1‑149 (MQ=255)
                                                                                                                                         gaGGGGCTGCAGGGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTg         >  8:325143/2‑149 (MQ=255)
                                                                                                                                          aGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGc        <  7:384037/149‑1 (MQ=255)
                                                                                                                                          aGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGc        <  4:27241/149‑1 (MQ=255)
                                                                                                                                              gCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGCCAgc    >  6:103834/1‑149 (MQ=255)
                                                                                                                                                tGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGa                                                     <  7:216430/98‑1 (MQ=255)
                                                                                                                                                tGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGa                                                     >  8:216430/1‑98 (MQ=255)
                                                                                                                                                tGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGCCAgcgc  >  4:349666/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CAGCCTTGGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGCCAGCGC  >  NC_002947/900387‑900679

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