Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,903,164 G→T 16.7% 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,1640GT16.7% 94.3 / 16.7 54intergenic (‑128/‑140)PP_5175/PP_5176HlyD family secretion protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/26);  new base T (5/4);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-02

AGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAATGCCA  >  NC_002947/5903017‑5903309
                                                                                                                                                   |                                                                                                                                                  
aGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGAGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGa                                                                                                                                                   <  1:29748/149‑1 (MQ=255)
    gACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGCGccc                                                                                                                                               >  4:232183/1‑149 (MQ=255)
    gACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGccc                                                                                                                                               <  2:29073/149‑1 (MQ=255)
         aTCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGaa                                                                                                                                          <  1:172616/149‑1 (MQ=255)
            aaCGCCATGCTGGTTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGAaga                                                                                                                                        <  2:334571/148‑1 (MQ=255)
            aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGAaga                                                                                                                                        <  5:261898/148‑1 (MQ=255)
            aaCGCCATGCGGGCTCCTTGCTTGTGCACATATTCTTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGCGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAaga                                                                                                                                        <  2:286996/148‑1 (MQ=255)
                          tCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGa                                                                                                                                           <  6:126615/131‑1 (MQ=255)
                          tCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGa                                                                                                                                           >  5:126615/1‑131 (MQ=255)
                          tCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCATAGCCGGCGAGCCCCCCAAGAGGCCACCACAccc                                                                                                                           >  4:91570/1‑147 (MQ=255)
                            cttggtggTGCAAAAATGATTTGTGGGAGTCGACTGGCAGCGCAGGGCAGGGTTCGATAAACATCACTGAGGGGGGCCAGGCGGGCATTTTTGGTGCCTGCGGGGGCCGCTCCGCGGGCGAGCCCGCGAAGAGGCCGCCACAGCCAt                                                                                                                         <  8:20783/147‑1 (MQ=255)
                            cttgcttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGGTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCAt                                                                                                                         <  5:253266/147‑1 (MQ=255)
                            cttgcttgTGCAAAATTGATTTGTGGGATGCGCCTGGGAGCGTAGTTCAGGGTTCGATAAACATCCCTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCGCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCAt                                                                                                                         <  3:91570/147‑1 (MQ=255)
                               gcttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAagag                                                                                                                                       >  6:245054/1‑130 (MQ=255)
                               gcttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAagag                                                                                                                                       <  5:245054/130‑1 (MQ=255)
                                 ttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCAAATAGCGGGCGAGCCCGCGAAAGGGCCGACACACCACTGTAtaaa                                                                                                                  >  4:186917/1‑146 (MQ=255)
                                        caAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGc                                                                                                           <  3:232183/149‑1 (MQ=255)
                                         aaaTGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGc                                                                                                           >  1:101182/1‑148 (MQ=255)
                                              ttttGTGTGAGGCGACTGTGAGCGTTGTTCTGGGTCCAAAAATCATCCCTGCGGGGGGCCTTGCGGTCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCatat                                                                                                      <  8:243246/147‑1 (MQ=255)
                                                        tGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCa                                                                                                   <  2:41201/141‑1 (MQ=255)
                                                        tGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCa                                                                                                   >  1:41201/1‑141 (MQ=255)
                                                        tGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCACATACCGGGCGAGCCCGCGAAGAGGCCTCCACCGCCAAGTATAAGTATATGCCATAACCAGAAACCa                                                                                            >  4:279532/1‑148 (MQ=255)
                                                           gACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGAag                                                                                                                                         <  4:238557/100‑1 (MQ=255)
                                                           gACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGAag                                                                                                                                         >  3:238557/1‑100 (MQ=255)
                                                           gACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTAGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCACTATCTGCCATAACCAGAAACCATCAg                                                                                        >  7:79814/1‑149 (MQ=255)
                                                                  aGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGGGAGCCCGCGAAGAGGCCTCCCCAGCCATGTTTCAGTATATGCCCTATCCAGAAAACATc                                                                                          >  4:242522/1‑140 (MQ=255)
                                                                  aGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATc                                                                                          <  3:242522/140‑1 (MQ=255)
                                                                     gTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCAAATAGCGGGCGCGCCCGCCCAGACGCCGCCACAGCCCAGTTTAACTATATGCCAAAAACAAAAACCAACAAGCAAAGGcc                                                                               >  8:10717/1‑148 (MQ=255)
                                                                         ttCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTAGCCGGCGCGCCCGCGGAGAGGCCTCCCCCGCCAAGTATCAGTATCTGCCATATACAGAAACCCTCAGGCAACGGCCATccc                                                                          >  8:196103/1‑149 (MQ=255)
                                                                           cAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATATTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCGCCACAGCCATGTATCACTATCTGCCATAACCAGAACCCAACCGGCAACCGCCAAcccc                                                                         >  6:262466/1‑148 (MQ=255)
                                                                              ggTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGCGCCCGCGAAGAGGCCACCCCAACCCTGTATCAGTATCTGCCATATCCAGAAACCAACAGCCATACGCCCACCCCggt                                                                      >  8:154806/1‑146 (MQ=255)
                                                                                        aCATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTAc                                                          <  7:10717/149‑1 (MQ=255)
                                                                                            cACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGTGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCttt                                                      >  8:47582/1‑148 (MQ=255)
                                                                                             aCTGCCGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCttt                                                      <  2:281222/148‑1 (MQ=255)
                                                                                             aCTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCttt                                                      <  3:279532/148‑1 (MQ=255)
                                                                                                     gggCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAg                                             <  7:208100/149‑1 (MQ=255)
                                                                                                         cATGGTGGCATTTTTGGGGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTc                                         <  7:154806/149‑1 (MQ=255)
                                                                                                          aTGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCCAACCGCTACCTTTCACAGGGAGGCTc                                         <  1:178572/149‑1 (MQ=255)
                                                                                                             cTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTcaca                                      <  3:186917/148‑1 (MQ=255)
                                                                                                               ggCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTc                                   <  2:101182/149‑1 (MQ=255)
                                                                                                                gCATTTTTGGTGCCGGCGGGGGCCTCTTCGCGGGGGAGCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTc                                   <  5:27766/148‑1 (MQ=255)
                                                                                                                      ttGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGaa                             <  7:196103/148‑1 (MQ=255)
                                                                                                                                 ggtggCGTCTTCGCGGGCGAGCCCTCGAAGATGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGcc                 <  8:79814/146‑1 (MQ=255)
                                                                                                                                    ggCCTCTTCGCGGGCGAGCCCTCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACGATCGGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTGACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGtt              <  5:262466/149‑1 (MQ=255)
                                                                                                                                     gCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGt               >  4:336649/1‑147 (MQ=255)
                                                                                                                                     gCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTa             <  7:25301/149‑1 (MQ=255)
                                                                                                                                       ctctTCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGc           >  7:140716/1‑149 (MQ=255)
                                                                                                                                        tctTCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGcc          >  7:100108/1‑149 (MQ=255)
                                                                                                                                        tctTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGcc          >  7:92991/1‑149 (MQ=255)
                                                                                                                                        tctTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGcc          >  5:107083/1‑149 (MQ=255)
                                                                                                                                        tctTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGcc          >  1:161738/1‑149 (MQ=255)
                                                                                                                                          ttCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTTACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGa        >  3:276737/1‑149 (MQ=255)
                                                                                                                                          ttCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATTCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGa        >  7:152942/1‑149 (MQ=255)
                                                                                                                                                 gCGATCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACATGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAATGCCa  <  8:299068/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAATGCCA  >  NC_002947/5903017‑5903309

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: