Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,282,868:1 +T 7.3% intergenic (‑217/+12) PP_2010 ← / ← PP_2011 cytochrome b561/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,282,8681.T 7.3% 204.8 / 10.7 55intergenic (‑217/+12)PP_2010/PP_2011cytochrome b561/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (35/16);  new base T (2/2);  total (37/18)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

TGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTCCAGCACATGCAGCCAG  >  NC_002947/2282744‑2283015
                                                                                                                             |                                                                                                                                                   
tgtgAAAACGGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCGGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCa                                                                                                                            <  1:316817/150‑1 (MQ=255)
                                          gCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGcct                                                                                  >  1:1794661/1‑150 (MQ=255)
                                                  ccACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAg                                                                          >  2:422583/1‑150 (MQ=255)
                                                   cacaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCATTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAgg                                                                         >  1:1133534/1‑150 (MQ=255)
                                                   cacaGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAgg                                                                         >  2:1583828/1‑150 (MQ=25)
                                                       gTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCt                                                                                      >  2:1405793/1‑133 (MQ=25)
                                                       gTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCt                                                                                      <  1:1405793/133‑1 (MQ=25)
                                                          cAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCTACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc                                                                  >  2:1250605/1‑150 (MQ=255)
                                                          cAACTGCACAGTCTTGGAAGCCTCAGCCGTTAGACCATCGACGAGGGCGCGCCCGAGAACGGGCCGAAGCGCGACGGATCACGCAGTGGGCGGATAGGGCCAGGCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc                                                                  >  2:689783/1‑150 (MQ=255)
                                                           aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAg                                                                                                                       <  1:57842/96‑1 (MQ=255)
                                                           aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAg                                                                                                                       >  2:57842/1‑96 (MQ=255)
                                                           aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                                                                 >  2:1045273/1‑150 (MQ=255)
                                                           aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                                                                 >  1:1584726/1‑150 (MQ=255)
                                                           aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                                                                 >  1:1690819/1‑150 (MQ=255)
                                                           aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                                                                 >  2:1450118/1‑150 (MQ=37)
                                                            aCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGc                                                                        <  2:242045/142‑1 (MQ=255)
                                                            aCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGc                                                                        >  1:242045/1‑142 (MQ=255)
                                                             cTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGGCTGGCAACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCgg                                                                >  1:1257038/1‑150 (MQ=17)
                                                              tGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                                                              >  2:103600/1‑150 (MQ=255)
                                                              tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACGAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                                                              >  1:293520/1‑150 (MQ=37)
                                                              tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGc                                                                                    <  1:839850/128‑1 (MQ=25)
                                                              tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGc                                                                                    >  2:839850/1‑128 (MQ=25)
                                                              tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACTATCGGCa                                                              >  2:752798/1‑150 (MQ=37)
                                                              tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc                                                                  <  2:169680/146‑1 (MQ=255)
                                                              tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc                                                                  >  1:169680/1‑146 (MQ=37)
                                                              tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                                                              >  2:1675173/1‑150 (MQ=37)
                                                              tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGGCTGGCAACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                                                                 >  1:1438830/1‑148 (MQ=18)
                                                              tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGGCTGGCAACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                                                                 <  2:1438830/148‑1 (MQ=18)
                                                                  caGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAg                                                                                                                                      <  1:484639/74‑1 (MQ=255)
                                                                  caGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAg                                                                                                                                      >  2:484639/1‑74 (MQ=255)
                                                                     tCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCa                                                       >  1:59456/1‑150 (MQ=255)
                                                                     tCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCa                                                       >  1:253246/1‑150 (MQ=255)
                                                                         ggAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCggg                                                                                                                              <  1:1297499/75‑1 (MQ=255)
                                                                         ggAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCggg                                                                                                                              >  2:1297499/1‑75 (MQ=255)
                                                                                  tgCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCCGGCACCATCGGCATCAGGCACCAGTCGTACCAg                                          >  2:309527/1‑150 (MQ=255)
                                                                                  tgCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAg                                          >  2:621145/1‑150 (MQ=255)
                                                                                     aGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCaa                                       <  1:422583/150‑1 (MQ=255)
                                                                                      gTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGGCTGGCAACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCaa                                       <  2:1257038/150‑1 (MQ=37)
                                                                                          gACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACGAGAAACACCCGTTGTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGGACCAGCAAGTgcg                                  >  2:1032989/1‑150 (MQ=255)
                                                                                           aCTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCATTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTgcgc                                 <  2:1133534/150‑1 (MQ=255)
                                                                                               tgGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATc                                                            <  2:536520/119‑1 (MQ=255)
                                                                                               tgGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATc                                                            >  1:536520/1‑119 (MQ=255)
                                                                                               tgGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTACGCCgaa                             >  1:28559/1‑148 (MQ=255)
                                                                                                    gCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCctg                        <  1:103600/150‑1 (MQ=255)
                                                                                                      gggCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCctgct                      >  2:91200/1‑150 (MQ=255)
                                                                                                             cccGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCATTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAgg                                                         >  1:1482081/1‑108 (MQ=255)
                                                                                                             cccGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAgg                                                         <  2:1482081/108‑1 (MQ=255)
                                                                                                             cccGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTc                  <  2:1263584/147‑1 (MQ=255)
                                                                                                             cccGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTc                  >  1:1263584/1‑147 (MQ=255)
                                                                                                              ccGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTCCAGc              >  1:266880/1‑150 (MQ=255)
                                                                                                               cgcgAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTCCAGca             <  1:1583828/150‑1 (MQ=255)
                                                                                                                        gCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTAc                                             <  2:1744585/109‑1 (MQ=255)
                                                                                                                        gCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTAc                                             >  1:1744585/1‑109 (MQ=255)
                                                                                                                        gCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTCCAGCACATGCAGcc    >  1:1037134/1‑150 (MQ=255)
                                                                                                                          cGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTCCAGCACATGCAGCCAg  >  2:26638/1‑150 (MQ=255)
                                                                                                                             |                                                                                                                                                   
TGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTACCAGCAAGTGCGCCGGTGCCTGCTGCTCCAGCACATGCAGCCAG  >  NC_002947/2282744‑2283015

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