Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,087,755 C→T 9.3% intergenic (+404/‑22) PP_3597 → / → PP_3598 amino acid ABC transporter ATP‑binding protein/peptidase C26

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,087,7550CT 9.3% 151.1 / 3.4 53intergenic (+404/‑22)PP_3597/PP_3598amino acid ABC transporter ATP‑binding protein/peptidase C26
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base T (2/3);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGACTACCAGGTCATGACCCACAAATACATTACCCC  >  NC_002947/4087614‑4087901
                                                                                                                                             |                                                                                                                                                   
cgtcgtTTCACATATCGAAAGGCTCTCGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTg                                                                                                                                              <  2:324598/149‑1 (MQ=255)
    gtTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTTCGACCTGTcgt                                                                                                                                          <  8:135962/149‑3 (MQ=255)
    gtTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTTCGACCTGTcgt                                                                                                                                          <  3:145965/149‑3 (MQ=255)
     tttGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTTCGACCTGt                                                                                                                                             >  4:148081/1‑145 (MQ=255)
     tttGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTTCGACCTGt                                                                                                                                             <  3:148081/145‑1 (MQ=255)
         aCATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGt                                                                                                                                     >  4:114418/1‑149 (MQ=255)
           atatCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGATTCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGa                                                                                                                                   <  8:281332/149‑1 (MQ=255)
             atCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACa                                                                                                                                 >  2:270609/1‑149 (MQ=255)
              tCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACat                                                                                                                                <  1:270609/149‑1 (MQ=255)
              tCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACa                                                                                                                                 >  2:318193/1‑148 (MQ=255)
                    ggCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGa                                                                                                                                   <  5:250834/140‑1 (MQ=255)
                    ggCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACATCTCCAAAGCTGACaactcgg                                                                                                                           >  6:250834/1‑142 (MQ=255)
                      ctctGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCaa                                                                                                                        >  5:307775/1‑149 (MQ=255)
                      ctctGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCaa                                                                                                                        <  1:95346/149‑1 (MQ=255)
                                cTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACAtt                                                                                                              <  6:307775/149‑1 (MQ=255)
                                  aTCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTgg                                                                                                            >  3:206404/1‑149 (MQ=255)
                                     tACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTcc                                                                                                                          <  7:276346/132‑1 (MQ=255)
                                     tACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTcc                                                                                                                          >  8:276346/1‑132 (MQ=255)
                                      aCAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCaa                                                                                                                        <  6:199621/133‑1 (MQ=255)
                                      aCAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCaa                                                                                                                        >  5:199621/1‑133 (MQ=255)
                                               gACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCt                                                                                               >  2:3775/1‑149 (MQ=255)
                                                    gAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCATCAGCAACATTGGcaac                                                                                                        >  4:102991/1‑135 (MQ=255)
                                                    gAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGcaac                                                                                                        <  3:102991/135‑1 (MQ=255)
                                                      aTCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTg‑c                                                                                       <  1:318193/149‑1 (MQ=255)
                                                       tcgtcgCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACTAGCAACCTTCCCTg‑ct                                                                                      >  4:49676/1‑148 (MQ=255)
                                                              ctctGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTg                                                                                                             <  4:144642/120‑1 (MQ=255)
                                                              ctctGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTg                                                                                                             >  3:144642/1‑120 (MQ=255)
                                                               tctGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCcg                                                                              >  7:272667/1‑149 (MQ=255)
                                                               tctGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCcg                                                                              >  6:177387/1‑149 (MQ=255)
                                                                   aCGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCcgtc                                                                            <  8:272667/147‑1 (MQ=255)
                                                                         cGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCtgat                                                                    >  3:116511/1‑149 (MQ=255)
                                                                           gCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCtgatga                                                                  >  7:306862/1‑149 (MQ=255)
                                                                                 cTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTGTCGCAAACCCCTCGTCCTGATGACCAt                                                              <  1:3775/148‑1 (MQ=255)
                                                                                  tCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATggg                                                           >  6:242650/1‑149 (MQ=255)
                                                                                            aGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAagag                                                 <  4:206404/149‑1 (MQ=255)
                                                                                               cagcaCCCAAGGCCCCGAGCCGTGCACCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAgcgc                                              <  8:306862/149‑1 (MQ=255)
                                                                                                        aGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCAt                                     <  5:242650/149‑1 (MQ=255)
                                                                                                        aGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCAt                                     <  4:116511/149‑1 (MQ=255)
                                                                                                            cccGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTc                                                                                                                           >  5:262023/1‑60 (MQ=255)
                                                                                                            cccGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTc                                                                                                                           <  6:262023/60‑1 (MQ=255)
                                                                                                            cccGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGAcc                                 >  5:55194/1‑148 (MQ=255)
                                                                                                                   cGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGACTACCAGGt                          >  1:306912/1‑149 (MQ=255)
                                                                                                                    gTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGACTACCAGGTc                         <  3:49676/149‑1 (MQ=255)
                                                                                                                       cgcCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGACTACCAGGTCATg                      <  2:306912/149‑1 (MQ=255)
                                                                                                                        gcCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGACTACCAGGTCATGa                     <  3:114418/149‑1 (MQ=255)
                                                                                                                            gCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTcc                                                                        <  8:309702/94‑1 (MQ=255)
                                                                                                                            gCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTcc                                                                        >  7:309702/1‑94 (MQ=255)
                                                                                                                                 cttgggACAGGTTCGACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGACTACCAGGTCATGACCCACaaa             >  5:286082/6‑148 (MQ=255)
                                                                                                                                  caacGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGa                                   <  2:321594/125‑1 (MQ=255)
                                                                                                                                  caacGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGa                                   >  1:321594/1‑125 (MQ=255)
                                                                                                                                     cGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGACTACCAGGTCATGACCCACAAATACAt        >  8:338410/1‑149 (MQ=255)
                                                                                                                                           gTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGACTAcc                              >  2:306504/1‑121 (MQ=255)
                                                                                                                                           gTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGACTACCAGGTCATGACCCACAAATACATTAcccc  >  1:286624/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGAGCCGGCAGTGCTCGGGGCTACAGGCAGCACCCAAGGCCCCGAGCCGTGCGCCTGCACAACAACGACAGGTCGAACCTGTCCCAAGGTGACATATGTCCAACAGCAACATTGGCAACAAGCAACCTTCCCTG‑CGCAAACCCGTCGTCCTGATGACCATGGGCAGCCAAGAGCGCAAAGGCCATGACTACCAGGTCATGACCCACAAATACATTACCCC  >  NC_002947/4087614‑4087901

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