Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,644,838 A→C 9.5% intergenic (‑97/+107) lon‑I ← / ← gcd DNA‑binding ATP‑dependent protease/quinoprotein glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,644,8380AC 9.5% 139.6 / 7.5 53intergenic (‑97/+107)lon‑I/gcdDNA‑binding ATP‑dependent protease/quinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (29/19);  new base C (2/3);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

TTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAATGATGTAATCGCCCATCTT  >  NC_002947/1644694‑1644986
                                                                                                                                                |                                                                                                                                                    
ttCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGcagcag                                                                                                                                                  >  3:273295/1‑149 (MQ=255)
        gtgcggcAGGGGTTTCGGGGAACACGTGCTGGGCCCTCATGGGCCACCTGGGGTAATGACTAGGGGGGTTCGTTGGGGCGGGCAGGGAGGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGaa                                                                                                                                          <  3:131869/143‑1 (MQ=19)
                  gttttggggtAATTCCTGCTGTTCGGTCATGGGGAACTTGCGTTAATTAATTTGGGTGTTTGATGGGGGGGGTTGGTTTGGTTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGc                                                                                                                                <  4:273295/143‑1 (MQ=19)
                        ggggAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGCGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTc                                                                                                                            >  5:176788/1‑147 (MQ=255)
                          ggAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTc                                                                                                                            >  8:196883/1‑145 (MQ=255)
                            aaaTCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTc                                                                                                                            >  7:193050/1‑143 (MQ=255)
                             accccTGCGGGTCGCGCATGGGGCACCTGCGTGAATGACTAGTGTGGTTTTATGGGGCGGGGAGGTATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATc                                                                                                                      <  5:103971/145‑1 (MQ=38)
                                cgtcgtgGTCCCTCTTGGGGCACCGGGGTGAAGGCATTTGTGGCTTCTAGGGGGGGGGGTTGGTGGGGTTTCAATGGTTCGGTAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAg                                                                                                                         <  6:35965/137‑1 (MQ=19)
                                ccTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAg                                                                                                                         >  5:35965/1‑142 (MQ=255)
                                                   gCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTTGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGa                                                                                                >  5:269393/1‑148 (MQ=255)
                                                          tgTGAATTACTTTGGTGGTTCGTTGGGCGGGGGAGGGTTGGGTTTCAATGGTGCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGa                                                                                         <  8:193050/147‑1 (MQ=38)
                                                               aTGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAg                                                                                   >  1:159897/1‑149 (MQ=255)
                                                                tGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAgg                                                                                  >  8:103728/1‑149 (MQ=255)
                                                                 gACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGc                                                                                 >  7:176454/1‑149 (MQ=255)
                                                                    tATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCg                                                                               >  8:213024/1‑148 (MQ=255)
                                                                    tATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCg                                                                               >  1:175360/1‑148 (MQ=255)
                                                                       ggTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCTGCGCAGcccaa                                                                                                                                             >  5:342500/1‑81 (MQ=38)
                                                                       ggTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCTGCGCAGcccaa                                                                                                                                             <  6:342500/83‑3 (MQ=38)
                                                                            ttAGATGGGGCGGGGAGGGATGGGTTACAATGGTTCGGTATCAGCCTTGTGTCGCGATTGGGCTGAGAAGCAGCCCCAGAATCCATACAGCATTCTAGATCGTTTGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAg                                                                      >  1:6139/1‑149 (MQ=255)
                                                                                   gggCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAg                                                               >  8:127661/1‑149 (MQ=255)
                                                                                        gggagggaTGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGc                                                          >  8:169475/1‑149 (MQ=255)
                                                                                                 gggTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGgt                                                 <  7:213024/149‑1 (MQ=255)
                                                                                                    tttCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGGATTCTAGATCGTTGGGGCCGCTGCGCGACTGAAACGCGACACAAGCCCGCTCCCGCAGGGGTGAGCCTGCGGGAGCGGTGTc                                              >  8:33285/1‑149 (MQ=39)
                                                                                                    tttCAATGGTTCGGGAGCAGCCTTGTGTCGCGAGTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTc                                              <  3:206236/149‑1 (MQ=255)
                                                                                                     ttCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCg                                             >  8:134469/1‑149 (MQ=255)
                                                                                                      tCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGc                                            >  3:55176/1‑149 (MQ=255)
                                                                                                      tCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGc                                            >  6:95814/1‑149 (MQ=255)
                                                                                                        aaTGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAgg                                                                     <  7:358123/122‑1 (MQ=255)
                                                                                                        aaTGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAgg                                                                     >  8:358123/1‑122 (MQ=255)
                                                                                                         aTGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCttt                                         <  7:103728/149‑2 (MQ=255)
                                                                                                         aTGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCtt                                          >  8:344199/1‑148 (MQ=255)
                                                                                                           ggTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCAGTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGGCACAAGGCCGCTCCGTCAGGGTTGAGCCTGCGGGAGCGGGGTCGcggacg                                       >  8:180453/1‑144 (MQ=39)
                                                                                                                gggAGCAGCCTTGTGTCGCGATTTGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGc                                   <  6:269393/148‑1 (MQ=255)
                                                                                                                 ggAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCt                                  >  7:225950/1‑148 (MQ=255)
                                                                                                                  gAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTaa                                >  5:35953/1‑149 (MQ=255)
                                                                                                                    gcagcCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAAtt                              >  1:55295/1‑149 (MQ=255)
                                                                                                                     cagcCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAAttt                             <  7:134469/149‑1 (MQ=255)
                                                                                                                        ccTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGCCCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCCCTTACACGGCTAATTTGt                           >  5:279148/1‑148 (MQ=255)
                                                                                                                        ccTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGt                           >  6:301746/1‑148 (MQ=255)
                                                                                                                         cTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTaa                         <  8:225950/149‑1 (MQ=255)
                                                                                                                            gtgtCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCAGCGGGAGCGGTGTCGCTTACTCGGCTAATATGTAAGc                       >  8:159321/1‑148 (MQ=255)
                                                                                                                                cgcgATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAatg                   <  7:174551/148‑1 (MQ=255)
                                                                                                                                cgcgATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAatg                   >  8:174551/1‑148 (MQ=255)
                                                                                                                                 gcgATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAatga                  <  7:127661/148‑1 (MQ=255)
                                                                                                                                 attctggggCTGCTGCGCAGCCCCAGAATCCATCCAGCATTCTTGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAatgat                 <  5:316682/143‑1 (MQ=255)
                                                                                                                                  cgATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAatgatg                <  8:176454/149‑1 (MQ=255)
                                                                                                                                  cgATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAatgatg                <  7:344199/149‑1 (MQ=255)
                                                                                                                                     tggggCTGCTGCGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGa                                                      >  4:125139/3‑108 (MQ=39)
                                                                                                                                     tggggCTGCTGCGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGa                                                      <  3:125139/106‑1 (MQ=39)
                                                                                                                                      tGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAATGATGTAAt            >  8:337728/1‑149 (MQ=255)
                                                                                                                                        ggCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAATGATGTAATc           <  4:55176/148‑1 (MQ=255)
                                                                                                                                         gCTGCGCAGCAGCCCCAGAATCCATCCAACATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAATGATGTAATCGc         >  1:216608/1‑149 (MQ=255)
                                                                                                                                             cgcAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAATGATGTAATCGccc       <  5:95814/147‑1 (MQ=255)
                                                                                                                                             cgcAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAATGATGTAATCGCCCa      <  7:196883/148‑1 (MQ=255)
                                                                                                                                                agcagcCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAATGATGTAATCGCCCATCtt  >  7:4311/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAATCCATCCAGCATTCTAGATCGTTGGGGCCGCTGCGCGACCGAAACGCGACACAAGGCCGCTCCCTCAGGGTTGAGCCTGCGGGAGCGGTGTCGCTTACTCGGCTAATTTGTAAGCAATGATGTAATCGCCCATCTT  >  NC_002947/1644694‑1644986

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