Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,163,553 G→A 14.7% intergenic (‑48/‑134) pssA ← / → PP_3665 CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,163,5530GA14.7% 75.7 / 4.8 33intergenic (‑48/‑134)pssA/PP_3665CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (15/13);  new base A (2/3);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAATGGTCCACCCCGCTT  >  NC_002947/4163407‑4163702
                                                                                                                                                  |                                                                                                                                                     
caggcaggCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCgggg                                                                                                                                                    <  2:326889/150‑1 (MQ=255)
   gcaggcGCGGTAAGCCCTGGCATCGGGCAGCGACTGGAGCGCCGCGGGGGCCAGTGCGTAGCCTGGAAAAGTCTGCAGCATGGCCTTGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCa                                                                                                                                                 >  1:350299/1‑150 (MQ=255)
                 cccTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCCGCCGGCAGGCAGGGc                                                                                                                                                   >  1:171613/1‑134 (MQ=255)
                 cccTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCCGCCGGCAGGCAGGGc                                                                                                                                                   <  2:171613/134‑1 (MQ=255)
                   cTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCAt                                                                                                                                 >  1:416811/1‑150 (MQ=255)
                        aTCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCggg                                                                                                                            >  1:249880/1‑150 (MQ=255)
                                 caCCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGAt                                                                                                                   <  2:416811/150‑1 (MQ=255)
                                                 ggggCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGAcc                                                                                                   <  1:231589/150‑1 (MQ=255)
                                                 ggggCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCAGGAGGACCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGAcc                                                                                                   <  1:502991/150‑1 (MQ=255)
                                                      cATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATaa                                                                                              >  1:149726/1‑150 (MQ=255)
                                                          gCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATaacaac                                                                                          <  1:479490/150‑1 (MQ=255)
                                                                ccTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAGCAAGAAc                                                                                    >  1:182166/1‑150 (MQ=255)
                                                                 cTGGCAAGCACGCGAACATGGTTTCGGGGTTCACTGACGCTCCCCTTCCCTGGGCAATCACTGTCTGCCCTGCATGCCGGCGGGGCCATTCCTGGCCGCTCCATCCGGGGTGTATGATAGACCCGGCAACGACCCATAAAAACAAGAAcc                                                                                   >  2:57043/1‑150 (MQ=255)
                                                                  tGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCgggg                                                                                                                                                    >  1:586826/1‑84 (MQ=255)
                                                                  tGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCgggg                                                                                                                                                    <  2:586826/84‑1 (MQ=255)
                                                                             ccaccaTGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGctct                                                                       <  1:57043/150‑1 (MQ=255)
                                                                                  aTGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCCACGACCCATAACAACAAGAACCACGCGAGGCTCTGCttt                                                                  <  1:423972/150‑1 (MQ=255)
                                                                                      tttCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCa                                                              >  1:494873/1‑150 (MQ=255)
                                                                                                cACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGa                                                    <  2:601336/150‑1 (MQ=255)
                                                                                                 aCTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGac                                                   >  1:25126/1‑150 (MQ=255)
                                                                                                       gCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGa                                                    >  2:61341/1‑143 (MQ=255)
                                                                                                       gCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGa                                                    <  1:61341/143‑1 (MQ=255)
                                                                                                             cTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAg                                            >  1:606671/1‑145 (MQ=255)
                                                                                                             cTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAg                                            <  2:606671/145‑1 (MQ=255)
                                                                                                                ctctGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGc                                    <  2:350299/150‑1 (MQ=255)
                                                                                                                        aaTCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTtc                                                                 <  2:536630/113‑1 (MQ=255)
                                                                                                                        aaTCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTtc                                                                 >  1:536630/1‑113 (MQ=255)
                                                                                                                        aaTCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGc                                          <  2:560116/136‑1 (MQ=255)
                                                                                                                        aaTCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGc                                          >  1:560116/1‑136 (MQ=255)
                                                                                                                         aTCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAgg                           >  1:297396/1‑150 (MQ=255)
                                                                                                                                  ggCCCCGCCGGCAGGCAGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAAc                                                                                    <  2:571174/84‑1 (MQ=38)
                                                                                                                                  ggCCCCGCCGGCAGGCAGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAAc                                                                                    >  1:571174/1‑84 (MQ=38)
                                                                                                                                      ctgcctgcCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAATgg              >  1:361358/1‑150 (MQ=255)
                                                                                                                                       tgcctgccGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAgcgcgaccgg             <  2:249880/150‑7 (MQ=255)
                                                                                                                                       tgcctgccGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAATGGt             <  2:361358/150‑1 (MQ=255)
                                                                                                                                          ctgccGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAATGGTCCa          >  2:44685/1‑150 (MQ=255)
                                                                                                                                                  ggggCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAATGGTCCACCCCGCtt  >  2:405117/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CAGGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAATGGTCCACCCCGCTT  >  NC_002947/4163407‑4163702

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: