Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,163,549 C→A 13.1% intergenic (‑44/‑138) 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,5490CA13.1% 105.6 / 5.8 38intergenic (‑44/‑138)pssA/PP_3665CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (14/19);  new base A (3/2);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

GCTCGAACAGGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAAT  >  NC_002947/4163400‑4163688
                                                                                                                                                     |                                                                                                                                           
gcTCGAACAGGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCtgcctgcc                                                                                                                                             >  1:116299/1‑150 (MQ=255)
  tCGAACAGGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCcgg                                                                                                                                           >  2:186591/1‑150 (MQ=255)
         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                        >  1:798960/1‑134 (MQ=255)
                                                                                                                                     tGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCa                        <  2:798960/134‑1 (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      >  1:819121/1‑146 (MQ=255)
                                                                                                                                           cccTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAgcg      <  2:819121/146‑1 (MQ=255)
                                                                                                                                           cccTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAAt  >  1:694618/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GCTCGAACAGGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAAT  >  NC_002947/4163400‑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: ‑