Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,903,164 G→T 15.2% 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,1640GT15.2% 69.7 / 7.4 35intergenic (‑128/‑140)PP_5175/PP_5176HlyD family secretion protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/16);  new base T (3/2);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-01

CATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAA  >  NC_002947/5903025‑5903304
                                                                                                                                           |                                                                                                                                             
cATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGGGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGa                                                                                                                                      <  1:305061/149‑1 (MQ=255)
 aTCAACGCCTTGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCGCTTAGCCGGGGGGccccccg                                                                                                                                      >  7:44467/1‑144 (MQ=255)
 aTCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGa                                                                                                                                      >  4:217092/1‑148 (MQ=255)
    aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGGGAag                                                                                                                                   >  5:151674/1‑147 (MQ=255)
    aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAaga                                                                                                                                  <  4:113872/148‑1 (MQ=255)
    aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGCGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGAagag                                                                                                                                 <  2:207045/149‑1 (MQ=255)
      cGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGc                                                                                                                               <  4:259186/149‑1 (MQ=255)
              tGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTAGCCGGCGCGCCCCCGCAGCGGccggcccag                                                                                                                       >  1:180559/1‑142 (MQ=255)
                gCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGc                                                                                                                      >  4:101888/1‑148 (MQ=255)
                    gttgcttgTGCACAAATGATTTGTGGGTTGCGACTGTGAGCTTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCAt                                                                                                                   <  4:282465/146‑1 (MQ=255)
                    cttgcttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGCCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCAt                                                                                                                   <  3:79400/147‑1 (MQ=255)
                    attgcttgTCCACAAATTATTTGTGGGATGCGACTGTTAGCGTAGTTCAGGGTTCGATAAACATCACTTACGGTGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCTAAGAGGCCTCCACAGCCAt                                                                                                                   <  6:303093/146‑1 (MQ=255)
                       gcttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTAt                                                                                                               <  3:217092/148‑1 (MQ=255)
                           gtgCACAAATGATTGGGGGGAGGCGACTGGTCGCGTAGTTCTGGGTTCTAAAAACATCAGTGCGGGGGGCCATGGGGGCATTTTTGGTGCGGGGGGGGGCCGCTTCGCGGGCGCGCCGGCGAAGAGGCCTCCACAGCCATGTATCa                                                                                                             <  8:64132/146‑1 (MQ=255)
                              cacaAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCAATGCTGGCATTTTTGGTGCCTGAGGGGGCCTCTTAGCCGGCGAGCACGCGAAGAGGACTCCACAGCCATGTATAAGTAtat                                                                                                       >  5:137224/1‑147 (MQ=255)
                                 aaaTGATTTGTGGTATGCGACTTCTAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCTGGCGAGCCCGCGAaga                                                                                                                                  <  8:245239/119‑1 (MQ=255)
                                 aaaTGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGACATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAaga                                                                                                                                  >  7:245239/1‑119 (MQ=255)
                                      aTTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCCAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATc                                                                                               >  8:230057/1‑149 (MQ=255)
                                         tgtgGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGAaga                                                                                                                                  >  4:95524/1‑111 (MQ=255)
                                         tgtgGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCa                                                                                             <  6:151674/148‑1 (MQ=255)
                                           tgGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGaa                                                                                          <  7:230057/149‑1 (MQ=255)
                                                 gCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTATTAGCGGGCGAGCCCGCGAAGAGGCCTACACAGCCATGTAt                                                                                                               >  1:107183/1‑122 (MQ=255)
                                                          aGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCg                                                                            <  3:101888/148‑1 (MQ=255)
                                                                           gATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATcc                                                                     <  1:195481/138‑1 (MQ=255)
                                                                           gATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATACCCGTTCCAAcc                                                          >  7:255786/1‑149 (MQ=255)
                                                                              ctaCATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGc                                                        <  7:316241/146‑1 (MQ=255)
                                                                                 cATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTAGCCGGCCAGCCCCCGAACCGGCCTCCACACCCAAGAATCAGTATCAGCCCTATCCAGACACCAACAGGCCACCGCCCTACCCCATCCCACCGCAAc                                                     >  8:316241/1‑148 (MQ=255)
                                                                                   tCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAATAAGGCCTCCCCAGCCATGTATCAGTATATGCCATAACCAGAAACCAACAGGCATCGGCCTACCCCGTTCCAACCCCTACCt                                                   >  1:278045/1‑149 (MQ=255)
                                                                                      cTGACGGGGGCCTTGCTGGAATTTTTGGTGCCTGCGGGGGCCGCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTca                                               <  3:294575/149‑1 (MQ=255)
                                                                                                                                tctTCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTa        >  7:147267/1‑146 (MQ=255)
                                                                                                                                tctTCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTa        <  8:147267/146‑1 (MQ=255)
                                                                                                                                  ttCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGACATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCGAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCg    >  5:260527/1‑148 (MQ=255)
                                                                                                                                    cgcgGGCTAGCCCGCGAAGAGGCATCCACAGCCATGTATCAGAATCTGCCATATCCAGAAACCAACAGCCAACGGCCATCCCCGTTCAAACCGCTACCTTTCACAGGGAGGCTAACAATCAAGCAAGATCAAAATGCCGTTAGCCGaa  >  4:125081/1‑148 (MQ=255)
                                                                                                                                       gggCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGTCGaa  <  3:74329/145‑1 (MQ=255)
                                                                                                                                       gggCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGGCGaa  >  4:74329/1‑145 (MQ=255)
                                                                                                                                           |                                                                                                                                             
CATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAA  >  NC_002947/5903025‑5903304

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