Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,048,866 C→T 7.9% H473Y (CAC→TAC)  quiA → quinate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,048,8660CT 7.9% 144.5 / 2.8 51H473Y (CAC→TAC) quiAquinate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (24/22);  new base T (2/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

AACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGACCGCGCCACTGGCAAGCCGCT  >  NC_002947/4048725‑4049011
                                                                                                                                             |                                                                                                                                                 
aaCACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGa                                                                                                                                            >  2:312923/1‑149 (MQ=255)
  cacCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTATTCCAGACCGTGCACAACGAcc                                                                                                                                          >  4:29302/1‑149 (MQ=255)
     gctggtCCTGCCGATGCGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGCACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACctct                                                                                                                                       <  8:144672/144‑1 (MQ=255)
      gtgtTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTg                                                                                                                                      <  1:312923/149‑1 (MQ=255)
            cTGCCGATGGGCTCGTCCTCCACCGACATATAAGGCGTTGAGCGCACCAAACTGGACCACACTTACGGCGCTTCGGTGCTCGCTTTGGACCCCACCACCGGCACCCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGAAtt                                                                                                                                >  8:32607/1‑149 (MQ=255)
              gCCGATGGGCTCGTCCTCCACCGCCATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGTACAACGACCTCTGGGACTTc                                                                                                                               >  8:157941/1‑148 (MQ=255)
              gCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGAGTACAACGACCTCTGGCACtac                                                                                                                               >  8:162477/1‑146 (MQ=255)
                         cGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATg                                                                                                                   >  3:135661/1‑149 (MQ=255)
                           tcctccACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCa                                                                                                                 <  7:62226/149‑1 (MQ=255)
                                        tCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATc                                                                                                    >  6:172613/1‑149 (MQ=255)
                                            cGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGACTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACt                                                                                                >  1:243885/1‑149 (MQ=255)
                                             ggCGTTGAGCGCAGCAAACTGGACCACTCTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGTACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACtt                                                                                               <  7:162477/149‑1 (MQ=255)
                                             ggCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGTACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACtt                                                                                               <  7:157941/149‑1 (MQ=255)
                                              gCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACtt                                                                                               >  1:31493/1‑148 (MQ=255)
                                                         agcaAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGgacga                                                                                   >  8:1424/1‑149 (MQ=255)
                                                         agcaAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGCACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGtacga                                                                                   >  4:197687/1‑149 (MQ=255)
                                                           caAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGacgacg                                                                                 <  4:135661/149‑1 (MQ=255)
                                                             aaCTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGc                                                                               >  3:210781/1‑149 (MQ=255)
                                                                     cacaCTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGTCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGt                                                                       <  2:311339/149‑1 (MQ=255)
                                                                      acacTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCaa                                                                                                            <  8:87962/111‑1 (MQ=255)
                                                                      acacTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCaa                                                                                                            >  7:87962/1‑111 (MQ=255)
                                                                        acTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTg                                                                                                                         <  8:173083/96‑1 (MQ=255)
                                                                        acTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTg                                                                                                                         >  7:173083/1‑96 (MQ=255)
                                                                         cTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCt                                                                   <  4:210781/149‑1 (MQ=255)
                                                                          ttACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCAcc                                                                                           <  8:126830/124‑1 (MQ=255)
                                                                          ttACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCAcc                                                                                           >  7:126830/1‑124 (MQ=255)
                                                                             cGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCg                                                                <  2:31493/148‑1 (MQ=255)
                                                                               gcgcTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCggtg                                                             <  5:172613/149‑1 (MQ=255)
                                                                                 gcTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCggtggt                                                           <  4:299940/149‑1 (MQ=255)
                                                                                      ggTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATc                                                       <  6:51231/148‑1 (MQ=255)
                                                                                          cTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCAc                                                  >  3:265902/1‑149 (MQ=255)
                                                                                              cTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAg                                              <  8:276627/149‑1 (MQ=255)
                                                                                               tttGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAgg                                             >  1:122541/1‑149 (MQ=255)
                                                                                                     cGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGc                                       <  7:1424/149‑1 (MQ=255)
                                                                                                        caccacCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCggg                                        >  3:44669/1‑145 (MQ=255)
                                                                                                             ccGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTAc                               <  2:243885/149‑1 (MQ=255)
                                                                                                              cGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACg                              >  1:41586/1‑149 (MQ=255)
                                                                                                                gCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTg                            >  2:152172/1‑149 (MQ=255)
                                                                                                                   aCCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTc                         <  2:122541/149‑1 (MQ=255)
                                                                                                                       gAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGAcc                     >  2:287822/1‑149 (MQ=255)
                                                                                                                          aTGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGcc                                           <  7:42844/124‑1 (MQ=255)
                                                                                                                          aTGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGcc                                           >  8:42844/1‑124 (MQ=255)
                                                                                                                           tGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGACCgcgc                 <  1:287822/149‑1 (MQ=255)
                                                                                                                             ggTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCAACCTACCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGAGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGAACGCGCCa               >  7:49418/1‑149 (MQ=255)
                                                                                                                                    cAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGACCGCGCCACTGGc          >  4:278260/1‑147 (MQ=255)
                                                                                                                                    cAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGACCGCGCCACTGGCaa        <  4:282730/149‑1 (MQ=255)
                                                                                                                                     aGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATc                                                       <  4:161445/101‑1 (MQ=255)
                                                                                                                                      gACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCACTCGGTGCCTGCggtgg                                                            >  7:212513/1‑95 (MQ=255)
                                                                                                                                      gACCGTGAACAACGACCTCTGTGACTTCGACCTGCCCATGCAGGCAAGCCTTATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCggtgg                                                            <  8:212513/95‑1 (MQ=255)
                                                                                                                                        ccGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGACCGCGCCACTGGCAAGCCg    <  3:197687/149‑1 (MQ=255)
                                                                                                                                           tGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGACCGCGCCACTGGCAAGCCGCt  >  4:281422/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
AACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCCAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGACCGCGCCACTGGCAAGCCGCT  >  NC_002947/4048725‑4049011

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