Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,342,597 T→C 9.8% intergenic (+43/‑164) panC → / → pgi‑II pantothenate synthetase/glucose 6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,342,5970TC 9.8% 135.0 / 5.6 52intergenic (+43/‑164)panC/pgi‑IIpantothenate synthetase/glucose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base C (2/3);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-01

AATGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACG  >  NC_002947/5342459‑5342740
                                                                                                                                          |                                                                                                                                               
aaTGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTc                                                                                                                                                >  4:91131/1‑140 (MQ=255)
aaTGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTc                                                                                                                                                <  3:91131/140‑1 (MQ=255)
aaTGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGtt                                                                                                                                       >  6:9108/1‑149 (MQ=255)
aaTGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGtt                                                                                                                                       >  4:221350/1‑149 (MQ=255)
 aTGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGttt                                                                                                                                      <  2:291282/149‑1 (MQ=255)
    aTCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTAc                                                                                                                                    <  5:213714/148‑1 (MQ=255)
        agaCCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTa                                                                                                                               <  5:146816/149‑1 (MQ=255)
         gaCCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAg                                                                                                                              <  1:259560/149‑1 (MQ=255)
             aggtggtGATTGCGGCGGCTTACCTCGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAACAGAAGACCGCATAACGCTCTTGCGCAGCCTTGGAGCCCATTCGCGGGCAAGCCCCCTCCTTCACCGTTACCTGAACgggc                                                                                                                          >  4:210773/2‑149 (MQ=255)
               gtggtgATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCCCCTACAGGTAAACGTGCAGgagcgg                                                                                                                        <  8:208534/149‑1 (MQ=38)
                  gtgATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCACTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCATGCGGCCCAATCCCGGGCACACCCGCTCCTGGACGGGTTCCAGGAGgggcgggcg                                                                                                                     >  8:194623/1‑148 (MQ=255)
                    gATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGgcttg                                                                                                                   >  5:283870/1‑149 (MQ=255)
                     aTTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGgcttgc                                                                                                                  <  6:187099/149‑1 (MQ=255)
                        gcggcgGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGgcttgctt                                                                                                                <  5:9108/148‑1 (MQ=255)
                           gcggcTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTg                                                                                                                                 <  7:160785/128‑1 (MQ=255)
                           gcggcTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCAGCACGTTTACCTg                                                                                                                                 >  8:160785/1‑128 (MQ=255)
                                     tGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCCCTACTGCACCTTTACCTGTAGGGGCGGGCTTGCTTGCAATGaccccca                                                                                                   >  6:6927/1‑146 (MQ=255)
                                       ggTACCACGCGTTTTATCGACATCCTCTACCTGGATTTTGAAGAGAAGTCGGCATAACGGTCTTTGGCTGCCTTGCAGCCCATTCGCGGGCAACCCCGCTCCTTCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCaaaa                                                                                                 <  3:210773/148‑1 (MQ=255)
                                       ggTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCaaaa                                                                                                 <  2:110241/148‑1 (MQ=255)
                                        gTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATc                                                                                               >  1:40251/1‑149 (MQ=255)
                                          aaCACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATccc                                                                                             >  4:203905/1‑149 (MQ=255)
                                          aaCACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGAGCAAGCCCGCCCCTACAGGTAAACGTGCAGgagcggg                                                                                                                       <  2:98946/123‑1 (MQ=37)
                                          aaCACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGAGCAAGCCCGCCCCTACAGGTAAACGTGCAGgagcggg                                                                                                                       >  1:98946/1‑123 (MQ=37)
                                             acGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCt                                                                                          <  3:368047/149‑1 (MQ=255)
                                                     gATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATAc                                                                                   <  7:194623/148‑1 (MQ=255)
                                                     gATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCCCCTACAGGTAAACGTGCAGGAGCGGgctt                                                                                                                    <  6:268894/115‑1 (MQ=25)
                                                     gATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCCCCTACAGGTAAACGTGCAGGAGCGGgctt                                                                                                                    >  5:268894/1‑115 (MQ=25)
                                                       tCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCa                                                                                >  3:245315/1‑149 (MQ=255)
                                                         gACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCAtt                                                                              <  3:221350/149‑1 (MQ=255)
                                                                    ccTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGGCCCCAAAATAACCCTGGCCTACCCCTTACCAGCACCTg                                                                   >  5:6485/1‑149 (MQ=38)
                                                                       gCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGGAGCCCATTACCGGGGAACCCCGCGCCGGCACGGTTACCTGTTTGGGCGGGCATGGTTGCAATGCCCCCCAAATACCCCGGCCAGACCCAATACCACCACCAGGcc                                                                >  4:219590/1‑149 (MQ=255)
                                                                          tttGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGcttgcttg                                                                                                               >  5:152386/1‑99 (MQ=255)
                                                                          tttGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGcttgcttg                                                                                                               <  6:152386/99‑1 (MQ=255)
                                                                          tttGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCttt                                                             <  2:40251/149‑1 (MQ=255)
                                                                                agagaagaCCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTAGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCAGCCAAACCCATTACCAGCCCGTTGGCTTTGCCta                                                        >  4:149233/1‑148 (MQ=255)
                                                                                 gagaagaCCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCGCCTGCCCGTTTACCAGTATGGGCGGGCTTGCTTGCAATGCCCCCCAAATCCCCCAGCCATACCCATTCCCACCCCGTGTCCTTTGCCgat                                                       >  4:244120/1‑148 (MQ=255)
                                                                                  agaagaCCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATaa                                                     <  3:219590/149‑1 (MQ=255)
                                                                                      gaCCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCa                                                                                <  3:74961/118‑1 (MQ=255)
                                                                                                  gTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAAc                                     <  4:245315/149‑1 (MQ=255)
                                                                                                       gggCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCgg                                 <  6:283870/148‑1 (MQ=255)
                                                                                                                 gCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGcc                      >  7:190089/1‑149 (MQ=255)
                                                                                                                  cAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGcc                      <  3:149233/148‑1 (MQ=255)
                                                                                                                       cATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGCCTg                        >  5:7567/1‑141 (MQ=255)
                                                                                                                       cATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTg                        <  6:7567/141‑1 (MQ=255)
                                                                                                                         ttCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATAACCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCaa              <  3:203905/149‑1 (MQ=255)
                                                                                                                           cgcgGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGc                                                           <  7:138258/102‑1 (MQ=255)
                                                                                                                           cgcgGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGc                                                           >  8:138258/1‑102 (MQ=255)
                                                                                                                            gcgGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGc            <  8:190089/148‑1 (MQ=255)
                                                                                                                                gCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCTCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGcc           >  1:201521/1‑145 (MQ=255)
                                                                                                                                gCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCTCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGcc           <  2:201521/145‑1 (MQ=255)
                                                                                                                                    gCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTcaccac   >  6:38125/1‑149 (MQ=255)
                                                                                                                                     cccGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACg  >  8:179673/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                               
AATGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACG  >  NC_002947/5342459‑5342740

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