Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,903,164 G→T 11.8% intergenic (‑128/‑140) PP_5175 ← / → PP_5176 HlyD family secretion protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,903,1640GT11.8% 65.1 / 2.6 34intergenic (‑128/‑140)PP_5175/PP_5176HlyD family secretion protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (10/20);  new base T (2/2);  total (12/22)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

ATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAATGC  >  NC_002947/5903026‑5903307
                                                                                                                                          |                                                                                                                                               
aTCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGaa                                                                                                                                       <  2:315115/149‑1 (MQ=255)
   aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAagag                                                                                                                                    >  5:291022/1‑149 (MQ=255)
   aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAaga                                                                                                                                     <  8:242760/148‑1 (MQ=255)
   aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAaga                                                                                                                                     <  2:95089/148‑1 (MQ=255)
   aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTCGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGAaga                                                                                                                                     <  5:280047/148‑1 (MQ=255)
             tGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAg                                                                                                                          >  1:39090/1‑149 (MQ=255)
                   cttgcttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCAt                                                                                                                      <  2:39090/147‑1 (MQ=255)
                      gcttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTAt                                                                                                                  <  7:23583/148‑1 (MQ=255)
                       cttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCa                                                                                                                <  7:65612/149‑1 (MQ=255)
                              acaAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTATTAGCGGGCGAGCCCGGGAAGAGGCCTCCACAGCCATGTATCAGtttatg                                                                                                         >  3:101690/1‑144 (MQ=255)
                                aaaTGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGcc                                                                                                       <  1:270658/149‑1 (MQ=255)
                                      tttGTGGTAGGCCACTTTGCGCTTTGTTCAGGGTTCTATAAACATCACTGCCGGGGGCCTGGCGGGAATTTTTGGGGCCGGGGGGGGCCTCTTCGCGGGCGCGCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATc                                                                                                  <  4:159346/148‑1 (MQ=255)
                                               tGCGACTGGGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATGTTTGGTGCCTGCGTGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCAt                                                                                        <  8:251794/149‑1 (MQ=255)
                                                  tagtgtgAGTGTAGTTCTGGGTTCGATAACCATCATTGACGGGGGCCTTGCGGGCATTTTTGGGGCCTGCGGGGGCCTCTTCGCGGGGGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAg                                                                                     <  2:290509/146‑1 (MQ=255)
                                                  gACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCCCAGCCATGTATCAGTATATGCCATATACAGAAACCCTCAg                                                                                     >  1:95252/1‑149 (MQ=255)
                                                          gCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATc                                                                                >  6:244629/1‑146 (MQ=255)
                                                          gCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATc                                                                                <  5:244629/146‑1 (MQ=255)
                                                           cGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGGGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGcc                                                                            <  6:291022/149‑1 (MQ=255)
                                                           cGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCACTTCGCGGGCGCGCCCGCGAAGAGGCCGCCACAGCCATGTATCAGTAAATGGCATATCAAGAAACCAACAAGCATCGGc                                                                             >  6:231732/1‑148 (MQ=255)
                                                           cGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCACTTAGAGGGCGAGCCCGGCAAGCGGCCGCCACAGCCCTGTCTAAGTAACTGCCATTTCCAGAAACCCTCCGGCCTCGGc                                                                             >  6:148975/1‑148 (MQ=255)
                                                            gTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGCGCCCGCGAAGAGGCCGCCACCGCCATGTATCAGTATCTGCCCTATCCAGAACCCAACACGCCTCCGcc                                                                            >  2:250712/1‑148 (MQ=255)
                                                                   agtgTTCTATAAAAATCAATTCCGTGGGCCATGGTGGCATTTTTGTTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGt                                                                    <  8:220306/146‑1 (MQ=255)
                                                                                     cTGACGGGGTCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCGCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTca                                                  <  1:250712/149‑1 (MQ=255)
                                                                                            gggCCATGGTGGCATTTTTGGTGCCGGCGGGGGCGTCTTCGCGGGCGAGCCCGGGATGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAg                                           <  3:156145/149‑1 (MQ=255)
                                                                                                             ttGGTGCCTGCGGGGTCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGaa                           <  6:307885/148‑1 (MQ=255)
                                                                                                                           ggCCTCTTCGCGGGCGAGCCGGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCg              <  8:192812/147‑1 (MQ=255)
                                                                                                                           ggCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCg              >  7:192812/1‑147 (MQ=255)
                                                                                                                               tctTCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCCGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAg                          >  7:205158/1‑131 (MQ=255)
                                                                                                                               tctTCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAg                          <  8:205158/131‑1 (MQ=255)
                                                                                                                               tctTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTca                                      >  3:117268/1‑119 (MQ=255)
                                                                                                                               tctTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTca                                      <  4:117268/119‑1 (MQ=255)
                                                                                                                                 ttCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCg       >  8:245917/1‑148 (MQ=255)
                                                                                                                                   cgcgGGCGAGCCCTCGAAGAGGCCTCCACATCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGaa     <  1:53077/148‑1 (MQ=255)
                                                                                                                                     cgGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAATGc  <  5:231732/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                               
ATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAATGC  >  NC_002947/5903026‑5903307

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