Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,188,015:1 +A 11.8% intergenic (+76/‑134) PP_3683 → / → PP_3684 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,188,0151.A11.8% 112.1 / 12.2 34intergenic (+76/‑134)PP_3683/PP_3684hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (16/14);  new base A (2/2);  total (18/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 = 8.22e-01

TCTAGCGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCTG  >  NC_002947/4187870‑4188156
                                                                                                                                                  |                                                                                                                                             
tCTAGCGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTa                                                                                                                                           >  1:481153/1‑150 (MQ=255)
        tggttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCtt                                                                                                                                   <  2:451975/150‑1 (MQ=255)
         ggttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGCAAGCGGCTTg                                                                                                                                  <  1:670383/150‑1 (MQ=255)
                                      aCGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGc                                                                                                                                     >  1:596436/1‑118 (MQ=255)
                                      aCGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGc                                                                                                                                     <  2:596436/118‑1 (MQ=255)
                                        gATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCTTGGCCCGCCCCTGGCGGGGTTGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCg                                                                                                   <  1:428804/150‑1 (MQ=255)
                                        gATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGCAAGCGGCt                                                                                                                                    >  1:229385/1‑117 (MQ=255)
                                        gATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGCAAGCGGCt                                                                                                                                    <  2:229385/117‑1 (MQ=255)
                                        gATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCg                                                                                                   >  1:122040/1‑150 (MQ=255)
                                             cAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCaa                                                                                              >  2:394607/1‑150 (MQ=255)
                                                ggATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGc                                                                                           <  1:336843/150‑1 (MQ=255)
                                                   tGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGCAAGCGGCTTg                                                                                                                                  >  2:485282/1‑108 (MQ=255)
                                                   tGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGCAAGCGGCTTg                                                                                                                                  <  1:485282/108‑1 (MQ=255)
                                                                 ccTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCTAAACGGTGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATACTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGc                                                                          <  1:609689/150‑1 (MQ=255)
                                                                 ccTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGc                                                                          <  2:122040/150‑1 (MQ=255)
                                                                    aGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTc                                                                                 >  1:388496/1‑140 (MQ=255)
                                                                    aGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTc                                                                                 <  2:388496/140‑1 (MQ=255)
                                                                                ctcAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTCCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCa                                                           <  1:38397/150‑1 (MQ=255)
                                                                                        tcGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGa                                                   >  1:277358/1‑150 (MQ=255)
                                                                                        tcGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCGTGCCAAATCAAAAGCTCCTCTCTAGa                                                   >  1:120370/1‑150 (MQ=255)
                                                                                                   tCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTgc                                                                                       <  2:330545/103‑1 (MQ=255)
                                                                                                   tCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTgc                                                                                       >  1:330545/1‑103 (MQ=255)
                                                                                                   tCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGt                                        >  1:134913/1‑150 (MQ=255)
                                                                                                             ccccTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTa                                                     <  2:205496/127‑1 (MQ=255)
                                                                                                             ccccTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTa                                                     >  1:205496/1‑127 (MQ=255)
                                                                                                              cccTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCg                             >  2:212589/1‑150 (MQ=255)
                                                                                                               ccTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGc                            <  1:394607/150‑1 (MQ=255)
                                                                                                                    cGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGc                       >  2:445549/1‑150 (MQ=255)
                                                                                                                                gCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGc           >  1:530175/1‑150 (MQ=255)
                                                                                                                                gCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGc           >  1:393969/1‑150 (MQ=255)
                                                                                                                                  gTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCACCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCt         <  2:278676/150‑1 (MQ=255)
                                                                                                                                         aaGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATACTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCTg  >  1:718318/1‑150 (MQ=255)
                                                                                                                                         aaGCCGCTTAGCAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAg            >  1:176928/1‑140 (MQ=255)
                                                                                                                                         aaGCCGCTTAGCAAGCGGCTTGGCCCTAACGGGCTTTCATACTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAg            <  2:176928/140‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TCTAGCGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCTG  >  NC_002947/4187870‑4188156

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