Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,456,994 A→C 8.9% E22A (GAG→GCG)  PP_2150 → thiamine biosynthesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,456,9940AC 8.9% 160.2 / 7.1 56E22A (GAG→GCG) PP_2150thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (27/24);  new base C (3/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATT  >  NC_002947/2456847‑2457140
                                                                                                                                                    |                                                                                                                                                  
ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAg                                                                                                                                                   <  6:157029/149‑1 (MQ=255)
ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAg                                                                                                                                                   >  2:93445/1‑149 (MQ=255)
    tCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGAGAACCGGCACCCGAGAACCCCTCCTGTGCGCACAGCCGTACTGTTCATTCTGCTCATACTCACCGCATGCAACAAAGCCCCCACCATGGAccgc                                                                                                                                                >  1:246275/1‑149 (MQ=37)
        aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                           >  8:169373/1‑149 (MQ=255)
         ccggATTATCGCTCCCCCTTTAGGAGGGGGTTAACCCGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                           <  5:195302/146‑1 (MQ=255)
         cAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                           <  5:99637/148‑1 (MQ=255)
          aGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtg                                                                                                                                         <  1:40780/149‑1 (MQ=255)
            gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                                                                                                                                            <  8:139341/144‑1 (MQ=255)
            gATTATCGATGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCACACCCTGGAGCGCCTgg                                                                                                                                            >  7:139341/1‑144 (MQ=255)
                        cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                                                                                                                                            >  4:207313/1‑132 (MQ=255)
                        cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                                                                                                                                            <  3:207313/132‑1 (MQ=255)
                                            ccGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggt                                                                                                                                          >  1:124538/1‑114 (MQ=255)
                                            ccGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggt                                                                                                                                          <  2:124538/114‑1 (MQ=255)
                                             cgcgAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg                                                                                                                            >  3:56178/1‑127 (MQ=255)
                                             cgcgAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg                                                                                                                            <  4:56178/127‑1 (MQ=255)
                                                    agGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCcgc                                                                                                >  2:147324/1‑148 (MQ=255)
                                                      gCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAg                                                                                             >  3:309423/1‑149 (MQ=255)
                                                      gCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGACCGCGAg                                                                                             >  2:243642/1‑149 (MQ=255)
                                                          aGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGACTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCt                                                                                          >  5:269653/1‑149 (MQ=255)
                                                            cACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCtgg                                                                                        <  1:31333/149‑1 (MQ=255)
                                                               cGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCTCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCtggtgg                                                                                     <  1:243642/149‑1 (MQ=255)
                                                                      cGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTg                                                                              >  4:175959/1‑149 (MQ=255)
                                                                          aaaCCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCgg                                                                          <  3:175959/149‑1 (MQ=255)
                                                                              ccTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGcc                                                                        >  3:226047/1‑147 (MQ=255)
                                                                                 aCTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTAc                                                                                                                  <  2:46994/102‑1 (MQ=255)
                                                                                 aCTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTAc                                                                                                                  >  1:46994/1‑102 (MQ=255)
                                                                                  cTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTg                                                                  <  6:19448/149‑1 (MQ=255)
                                                                                  cTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATCGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAgg                                                                    <  1:93445/147‑1 (MQ=255)
                                                                                     tgCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCa                                                                <  7:169373/148‑1 (MQ=255)
                                                                                      gcgcACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGt                                                                                                        >  6:247107/1‑107 (MQ=255)
                                                                                      gcgcACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGt                                                                                                        <  5:247107/107‑1 (MQ=255)
                                                                                           caGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCg                                                         <  4:309423/149‑1 (MQ=255)
                                                                                                cTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCCCCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGCTGGCGGTgg                                                     >  4:287453/1‑148 (MQ=255)
                                                                                                  aCTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGaaa                                                  >  1:36916/1‑149 (MQ=255)
                                                                                                   cTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGaaa                                                  <  6:146339/148‑1 (MQ=255)
                                                                                                   cTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAAc                                                 >  6:133700/1‑149 (MQ=255)
                                                                                                            ctgctgCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGc                                         <  1:147324/148‑1 (MQ=255)
                                                                                                                  ctACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACAt                                  >  7:234023/1‑149 (MQ=255)
                                                                                                                      tCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACa                                   >  4:183296/1‑144 (MQ=255)
                                                                                                                      tCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACa                                   <  3:183296/144‑1 (MQ=255)
                                                                                                                        aCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGacca                            >  7:201404/1‑149 (MQ=255)
                                                                                                                        aCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGacc                             <  3:292959/148‑1 (MQ=255)
                                                                                                                               gCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGGAGATGGCGGt                                                       <  4:298149/115‑1 (MQ=255)
                                                                                                                               gCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGt                                                       >  3:298149/1‑115 (MQ=255)
                                                                                                                                 aaCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTc                   >  4:81407/1‑149 (MQ=255)
                                                                                                                                       ggCACCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCtc             >  3:287573/1‑148 (MQ=255)
                                                                                                                                        ggaccACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTAc            >  1:16211/4‑149 (MQ=255)
                                                                                                                                          cccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACcg          <  3:81407/149‑1 (MQ=255)
                                                                                                                                          cccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACcg          >  6:300515/1‑149 (MQ=255)
                                                                                                                                          cccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGAGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACcg          >  3:235916/1‑149 (MQ=255)
                                                                                                                                          accACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCt                                           <  1:35859/115‑1 (MQ=255)
                                                                                                                                          accACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCt                                           >  2:35859/2‑116 (MQ=255)
                                                                                                                                          accACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACcg          >  4:59601/2‑149 (MQ=255)
                                                                                                                                                 tGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATc                                 <  1:61465/119‑1 (MQ=255)
                                                                                                                                                 tGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATc                                 >  2:61465/1‑119 (MQ=255)
                                                                                                                                                   gAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGAtt  <  6:269653/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCC‑AGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATT  >  NC_002947/2456847‑2457140

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