Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,609,520 C→T 7.5% G238D (GGC→GAC)  PP_3183 ← SCO1/SenC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,609,5200CT 7.5% 151.1 / 2.5 53G238D (GGC→GAC) PP_3183SCO1/SenC family protein
Reads supporting (aligned to +/- strand):  ref base C (22/27);  new base T (2/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

TTCGAACAGCAGCATGGCCAAGGGGTCTTTTTCGGCCAGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCTGCTTCCAGTTGTGCAGGCTGTTGCCCAGGCGGTC  >  NC_002947/3609375‑3609659
                                                                                                                                                 |                                                                                                                                            
ttCGAACAGCAGCATGGCCAAGGGGTCTTTTTCGGCCAGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGTCCa                                                                                                                                           >  3:85363/1‑149 (MQ=255)
  cGAACAGCAGCATGGCCAAGGGGTCTTTTTCGGCCAGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACg                                                                                                                                         >  7:151854/1‑149 (MQ=255)
      cagcagCATGGCCAAGGGGTCTTTTTCGGCCAGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGTCCACGGTAt                                                                                                                                     <  4:85363/149‑1 (MQ=255)
      cagcagCATGGCCAAGGGGTCTTTTTCGGCCAGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTAt                                                                                                                                     >  5:384775/1‑149 (MQ=255)
               ggCCAAGGGGTCTTTTTCGGCCAGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGAc                                                                                                                             >  7:64540/1‑148 (MQ=255)
                       ggTCTTTTTCGGCCAGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACgg                                                                                                                    <  6:384775/149‑1 (MQ=255)
                        gtgtttttCGGCCAGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACggg                                                                                                                   >  3:270614/4‑149 (MQ=255)
                        gTCTTTTTCGGCCAGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACggg                                                                                                                   <  3:263753/149‑1 (MQ=255)
                                 ggCCAGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAat                                                                                                           <  4:261880/148‑1 (MQ=255)
                                   cctgCATCTGGGCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCTGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATc                                                                                                        <  8:146905/146‑1 (MQ=255)
                                   ccAGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATc                                                                                                        <  3:154792/149‑1 (MQ=255)
                                    cAGCATCTGGTCGGCCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGACGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGCTATGGCAGCACGAGCAACGGGTACGGAAAATCt                                                                                                       >  6:383775/1‑149 (MQ=255)
                                     aGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTg                                                                                                      <  8:151854/149‑1 (MQ=255)
                                           tGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGAAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTcgcc                                                                                                >  1:206868/1‑149 (MQ=255)
                                            ggTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTcgccg                                                                                               >  6:151043/1‑149 (MQ=255)
                                             gTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTcgccgc                                                                                              <  8:64540/149‑1 (MQ=255)
                                                  gCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATACCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCaa                                                                                                                       <  3:165460/119‑1 (MQ=255)
                                                  gCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATACCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCaa                                                                                                                       >  4:165460/1‑119 (MQ=255)
                                                  gCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGATCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGctgctg                                                                                         <  2:49425/149‑1 (MQ=255)
                                                   cACCTTCAGCCAGCGGACCAGGCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTTTGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTgg                                                                                        <  6:137738/149‑1 (MQ=255)
                                                     ccTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGt                                                                                      <  3:168111/149‑1 (MQ=255)
                                                       ttCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGt                                                                                      >  1:408885/1‑147 (MQ=255)
                                                           gccagcGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGt                                                                                                                                       <  1:313310/94‑1 (MQ=255)
                                                           gccagcGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGt                                                                                                                                       >  2:313310/1‑94 (MQ=255)
                                                            ccagcGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCGACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGa                                                                              <  2:206868/149‑1 (MQ=255)
                                                            ccagcGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGa                                                                              >  1:299848/1‑149 (MQ=255)
                                                                  gACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCg                                                                        <  2:408885/149‑1 (MQ=255)
                                                                     cagccagTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAg                                                                     >  1:239517/1‑149 (MQ=255)
                                                                       gccagTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCgctgct                                                                                          <  8:123639/127‑1 (MQ=255)
                                                                       gccagTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCgctgct                                                                                          >  7:123639/1‑127 (MQ=255)
                                                                       gccagTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGcc                                                                   >  4:178062/1‑149 (MQ=255)
                                                                            ttGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGc                                                              <  5:151043/149‑1 (MQ=255)
                                                                                cGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAg                                                          <  3:178062/149‑1 (MQ=255)
                                                                                  tCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGt                                                         >  4:138001/1‑148 (MQ=255)
                                                                                    aCGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCAt                                                      >  8:314598/1‑149 (MQ=255)
                                                                                    aCGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCAt                                                      <  3:294135/149‑1 (MQ=255)
                                                                                      gCTGCCGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCGCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGTCGAGCGGATTTGCGGAGCCTGAGCGTAGTCAtt                                                     <  1:48660/149‑1 (MQ=255)
                                                                                      gCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTAGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTg                                                    >  5:408862/1‑149 (MQ=255)
                                                                                                   gCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTgg                                       >  1:234914/1‑149 (MQ=255)
                                                                                                   gCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTgg                                       >  3:220736/1‑149 (MQ=255)
                                                                                                    cccAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGc                                      >  8:11961/1‑149 (MQ=255)
                                                                                                     ccAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGcc                                     >  7:271385/1‑149 (MQ=255)
                                                                                                       agcagGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCTg                                   <  7:11791/149‑1 (MQ=255)
                                                                                                           ggTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCTGCTTc                               <  2:299848/149‑1 (MQ=255)
                                                                                                            gTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCTGCTTc                               <  4:220736/148‑1 (MQ=255)
                                                                                                             tCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCTGCTTCCa                             <  7:314598/149‑1 (MQ=255)
                                                                                                                gggCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCTGCTTCCAGtt                          <  8:271385/149‑1 (MQ=255)
                                                                                                                          ctggctggCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCTGCTTCCAGTTGTGCAGGCTg                <  4:270614/149‑1 (MQ=255)
                                                                                                                             gctggcCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCTGCTTCCAGTTGTGCAGGCTGTTg             >  1:15432/1‑149 (MQ=255)
                                                                                                                                  cctggcTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCTGCTTCCAGTTGTGCAGGCTGTTGCCCAg        <  2:234914/149‑1 (MQ=255)
                                                                                                                                        tCGGTGTTGTCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCt                                    <  4:289518/115‑1 (MQ=255)
                                                                                                                                        tCGGTGTTGTCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCt                                    >  3:289518/1‑115 (MQ=255)
                                                                                                                                          ggTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCTGCTTCCAGTTGTGCAGGCTGTTGCCCAGGCGGTc  <  5:383775/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
TTCGAACAGCAGCATGGCCAAGGGGTCTTTTTCGGCCAGCATCTGGTCGGGCACCTTCAGCCAGCGGACCAGCCAGTTGGCGTCACGCTGACGGGTCACGCCCAGCAGGTCGGGGCCGATGCCTGGCTGGCCTGGCTCGGTGTTGCCCACGGTATGGCAGGACGAGCAACGGGTACGGAATATCTGCTCGCCGCTGCTGGGT‑GAGCGGATTTGCGGAGCCTGAGCGTAGTCATTGCTCATGGCGCTGGCCTGCTTCCAGTTGTGCAGGCTGTTGCCCAGGCGGTC  >  NC_002947/3609375‑3609659

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