Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,163,553 G→A 15.6% 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,5530GA15.6% 82.4 / 8.9 39intergenic (‑48/‑134)pssA/PP_3665CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (14/18);  new base A (3/3);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

GGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAAT  >  NC_002947/4163409‑4163688
                                                                                                                                                |                                                                                                                                       
ggcaggcGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGcc                                                                                                                                    >  2:428567/1‑150 (MQ=255)
 gcaggcGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCCTAGGTTCGGGGTTCACTGACGCGCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCa                                                                                                                                   >  1:656112/1‑150 (MQ=255)
  caggcGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCCGCCGGCAGGCAGGGccga                                                                                                                                  <  1:661209/150‑3 (MQ=255)
   aggcGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCAtt                                                                                                                                 >  2:516502/1‑150 (MQ=255)
    ggcGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTc                                                                                                                                <  2:534128/150‑1 (MQ=255)
    ggcGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCCGCCGGCAGGCAGGGccg                                                                                                                                   <  1:359460/147‑2 (MQ=255)
    ggcGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCCGCCGGCAGGCAGGGccg                                                                                                                                   >  2:359460/1‑146 (MQ=255)
       gcgGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTg                                                                                                                             <  2:667676/150‑1 (MQ=255)
               cccTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGc                                                                                                                     >  2:664114/1‑150 (MQ=255)
                      aTCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCggg                                                                                                              <  2:705299/150‑1 (MQ=255)
                          ggCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGgtgt                                                                                                          <  2:116299/150‑1 (MQ=255)
                              gcaCCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGa                                                                                                      <  1:428567/150‑1 (MQ=255)
                               caCCTTGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGAt                                                                                                     <  1:288378/150‑1 (MQ=255)
                                aCCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCTGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATc                                                                                                    <  1:186591/150‑1 (MQ=255)
                                aCCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATc                                                                                                    >  1:799573/1‑150 (MQ=255)
                                 ccTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCATGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTTGCCGCTGCATCCGGGGTGTATGATCg                                                                                                   >  2:270967/1‑150 (MQ=255)
                                 ccTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCg                                                                                                   >  2:585104/1‑150 (MQ=255)
                                  cTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGa                                                                                                  >  2:734448/1‑150 (MQ=255)
                                      aTCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGGGTATGATCGACCCg                                                                                              >  2:186113/1‑150 (MQ=255)
                                          cTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCaa                                                                                          <  2:656112/150‑1 (MQ=255)
                                              gggggCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCGTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGAc                                                                                      <  1:188876/150‑1 (MQ=255)
                                                   ccATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATa                                                                                 >  2:466863/1‑150 (MQ=255)
                                                                                       cGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTg                                             <  1:32222/150‑1 (MQ=255)
                                                                                              cACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGa                                      <  1:568734/150‑1 (MQ=255)
                                                                                                        cACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGc                            <  1:466863/150‑1 (MQ=255)
                                                                                                            ttCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCa                        <  1:516502/150‑1 (MQ=255)
                                                                                                                   tgCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGca                 <  2:799573/150‑1 (MQ=255)
                                                                                                                     cAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGcaca               <  1:186113/150‑1 (MQ=255)
                                                                                                                          aggaatggCCCCGCCGGCAGGCAGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAAc          >  2:762751/7‑150 (MQ=255)
                                                                                                                           cTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACa         <  1:385228/150‑1 (MQ=255)
                                                                                                                            tGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCa                        <  2:798960/134‑1 (MQ=255)
                                                                                                                            tGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCa                        >  1:798960/1‑134 (MQ=255)
                                                                                                                                ggCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCa    >  2:311143/1‑150 (MQ=255)
                                                                                                                                 gCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCaa   <  1:585104/150‑1 (MQ=255)
                                                                                                                                  ccccgccggcaGGCAGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGa            <  2:548089/136‑1 (MQ=255)
                                                                                                                                  ccccGCCGGCAGGCAGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGa            >  1:548089/1‑140 (MQ=255)
                                                                                                                                  cccTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAgcg      <  2:819121/146‑1 (MQ=255)
                                                                                                                                  cccTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAgcg      >  1:819121/1‑146 (MQ=255)
                                                                                                                                  cccTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAAt  >  1:694618/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                       
GGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAAT  >  NC_002947/4163409‑4163688

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: ‑