Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,971,420 G→T 8.7% intergenic (‑71/‑57) PP_5236 ← / → PP_5237 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,971,4200GT 8.7% 159.7 / 2.4 53intergenic (‑71/‑57)PP_5236/PP_5237hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/19);  new base T (3/2);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGTGCGTGCGGGCGAGCGTTTGCACGGAGCAATACCTCAAGTCCATTTGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAAGCGGCATTTTCCAGCAACAACGGAGTGTTCATGTCTC  >  NC_002947/5971274‑5971564
                                                                                                                                                  |                                                                                                                                                
ccTGGTGCGTGCGGGCGCGCGTTTGCACGGAGCAATACCTCAAGTCCATTTGACAGAAGACGCAGTGCCATCCCATGGCATCTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGccc                                                                                                                                               >  2:1722696/1‑149 (MQ=255)
 ctggtgcGTGCGGGCGAGCGTTTGCAAGGAGCAATACCTCAAGTCCATTTGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGa                                                                                                                                              >  1:1422101‑M2/8‑150 (MQ=255)
   ggtgcGTGCGGGCGAGCGGTTGCACGGAGCAATACCTCAAGTCCATTTGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTGAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAAc                                                                                                                                            >  1:1425560‑M2/6‑150 (MQ=255)
                  gCGTTTGCACGGAGCAATACCTCAAGTCCATTTGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAg                                                                                                                             <  1:789397/150‑1 (MQ=255)
                    gTTTGCACGGAGCAATACCTCAAGTCCATTTGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTTCCGAACCGCATGAGACCGCAggg                                                                                                                           >  1:84604/1‑150 (MQ=255)
                    gTTTGCACGGAGCAATACCTCAAGTCCATTTGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTTCCGAACCGCATGAGACCGCAggg                                                                                                                           >  1:84898/1‑150 (MQ=255)
                    gTTTGCACGGAGCAATACCTCAAGTCCATTTGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAggg                                                                                                                           >  1:1200043/1‑150 (MQ=255)
                    gTTTGCACGGAGCAATACCTCAAGTCCATTTGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAggg                                                                                                                           >  1:1368605/1‑150 (MQ=255)
                                            tCCATTTGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCac                                                                                                   <  2:1728838/150‑1 (MQ=255)
                                              cATTTGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCacga                                                                                                 >  1:738974/1‑150 (MQ=255)
                                               ctgtGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCacgac                                                                                                >  1:860649/4‑150 (MQ=255)
                                               aTTTGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCacgac                                                                                                <  2:738974/150‑1 (MQ=255)
                                                  tGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGca                                                                                             >  2:423093/1‑150 (MQ=255)
                                                     aagaagACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCacac                                                                                          >  2:1596706/1‑150 (MQ=255)
                                                     aagaagACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCacac                                                                                          >  1:833183/1‑150 (MQ=255)
                                                      agaagaCGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAAc                                                                                                                      >  2:1459804/1‑121 (MQ=255)
                                                      agaagaCGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAAc                                                                                                                      <  1:1459804/121‑1 (MQ=255)
                                                        aagaCGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCTTTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACttg                                                                                       <  1:765391/150‑1 (MQ=255)
                                                        aagaCGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAgg                                                                                                                            <  2:1080939/113‑1 (MQ=255)
                                                        aagaCGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAgg                                                                                                                            >  1:1080939/1‑113 (MQ=255)
                                                               aGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTGAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGcc                                                                                <  2:1425560/150‑1 (MQ=255)
                                                                     aTCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGtt                                                                          <  2:833183/150‑1 (MQ=255)
                                                                      tCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTg                                                                         >  2:1003293/1‑150 (MQ=255)
                                                                                  tGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCt                                                             >  2:1823190/1‑150 (MQ=255)
                                                                                   gCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTa                                                            >  2:356884/1‑150 (MQ=255)
                                                                                     ggTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTAtt                                                          >  1:1746203/1‑150 (MQ=255)
                                                                                     ggTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTAtt                                                          >  2:373514/1‑150 (MQ=255)
                                                                                       ttACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCt                                                        <  1:373514/150‑1 (MQ=255)
                                                                                           ccAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCa                                                                                                                              <  2:608163/76‑1 (MQ=255)
                                                                                           ccAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCa                                                                                                                              >  1:608163/1‑76 (MQ=255)
                                                                                           ccAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGc                                                    >  1:587617/1‑150 (MQ=255)
                                                                                                   acTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCCCAATTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGc                                            >  1:903224/1‑150 (MQ=255)
                                                                                                   acTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCACCACCCTGAAGCACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGc                                            >  1:228763/1‑150 (MQ=21)
                                                                                                    cTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCa                                           <  1:356884/150‑1 (MQ=255)
                                                                                                     tAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAg                                          >  2:936613/1‑150 (MQ=255)
                                                                                                     tAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGACGATGCTTTGCCCAGCCGCAg                                          >  1:592707/1‑150 (MQ=255)
                                                                                                           aTACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTTCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAAGCg                                    <  2:84604/150‑1 (MQ=255)
                                                                                                           aTACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTTCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAAGCg                                    <  2:84898/150‑1 (MQ=255)
                                                                                                           aTACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCa                                                 >  1:1130499/1‑137 (MQ=255)
                                                                                                           aTACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCa                                                 <  2:1130499/137‑1 (MQ=255)
                                                                                                            tACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAATCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAAGCgg                                   <  1:1454187/150‑1 (MQ=255)
                                                                                                             aCCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAACCGGc                                  >  2:1915359/1‑150 (MQ=255)
                                                                                                                     ccTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAAGCGGCATTTTCCa                          <  2:1368605/150‑1 (MQ=255)
                                                                                                                     ccTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAAGCGGCATTTTCCa                          <  2:592707/150‑1 (MQ=255)
                                                                                                                      cTGCGCAAAGCCTATTGGATCCATTTTTTTTTTTTTTTTTGAGACCGCAGGGAAAACCCATGCACTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCCGTGTCTATTCTTTGCCCAGCCGCAGGAAGCGGCATTTTCCAg                         >  2:1326082/1‑150 (MQ=11)
                                                                                                                                ccTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTTGCTATTCTTTGCCCAGCCGCAGGAAGCGGCATTTTCCAGCAACAACGGa               <  2:860649/150‑1 (MQ=255)
                                                                                                                                 cTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGc                                                    >  2:1926331/1‑112 (MQ=255)
                                                                                                                                 cTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGc                                                    <  1:1926331/112‑1 (MQ=255)
                                                                                                                                  tATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAAGCGGCATTTTCCAGCAACAACGGAgt             >  1:693207/1‑150 (MQ=255)
                                                                                                                                     tGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAAGCGGCATTTTCCAGCAACAACGGAGTGtt          >  1:871622/1‑150 (MQ=255)
                                                                                                                                       gATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAAGCGGCATTTTCCAGCAACAACGGAGTGTTCa        >  1:1646640/1‑150 (MQ=255)
                                                                                                                                        aTCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAAGCGGCATTTTCCAGCAACAACGGAGTGTTCAt       <  1:1596706/150‑1 (MQ=255)
                                                                                                                                             tttttGCCGAACCACACGACGCCCAAGGGGAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAAGCGGCATTTTCCAGCAACAACGGAGTGTTCATGtctc  <  2:903224/150‑1 (MQ=14)
                                                                                                                                                  |                                                                                                                                                
CCTGGTGCGTGCGGGCGAGCGTTTGCACGGAGCAATACCTCAAGTCCATTTGAAAGAAGACGCAGTGCCATCCCATGGCAACTGCGGTTACCCAACAACACTAAGTCATACCTGACGCCTGCGCAAAGCCTATTGGATCCATTTTTGCCGAACCGCATGAGACCGCAGGGAAAACCCCTGCCCTGAAAAACCACGACGCACACTTGTTGCGCCGGCGTTGGCAGCACTGGCTATTCTTTGCCCAGCCGCAGGAAGCGGCATTTTCCAGCAACAACGGAGTGTTCATGTCTC  >  NC_002947/5971274‑5971564

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