Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,904,244 Δ1 bp 7.5% coding (499/1593 nt) ilvA‑I ← threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,904,2440T. 7.5% 158.5 / 7.2 53coding (499/1593 nt)ilvA‑Ithreonine deaminase
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base . (2/2);  total (27/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 = 9.23e-01

TAGGCGGCGATCCCGGCGATCAAGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGCCTTCA  >  NC_002947/3904099‑3904390
                                                                                                                                                  |                                                                                                                                                   
tAGGCGGCGATCCCGGCGATCAAGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTcgc                                                                                                                                                  <  8:49813/149‑1 (MQ=255)
        gATCCCGGCGATCAAGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGc                                                                                                                                         >  5:101530/1‑149 (MQ=255)
         aTCCCGGCGATCAAGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCt                                                                                                                                        >  4:292801/1‑149 (MQ=255)
                   tCAAGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGc                                                                                                                              <  2:136142/149‑1 (MQ=255)
                     aaGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGCGCGACCTTCgtg                                                                                                                                        >  1:100609/1‑135 (MQ=255)
                     aaGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGCGCGACCTTCgtg                                                                                                                                        <  2:100609/137‑3 (MQ=255)
                     aaGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGCGCGACCTTCgt                                                                                                                                         <  5:110806/136‑2 (MQ=255)
                     aaGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGCGCGACCTTCgt                                                                                                                                         >  6:110806/1‑135 (MQ=255)
                         cgccACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGt                                                                                                                        >  1:276412/1‑149 (MQ=255)
                             aCCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGt                                                                                                                    <  6:140836/149‑1 (MQ=255)
                             aCCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGt                                                                                                                    >  6:256380/1‑149 (MQ=255)
                              ccGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTg                                                                                                                   <  3:79395/149‑1 (MQ=255)
                               cGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGTGCGCTGACGGAGGATTTCCATGGCGACGGTTCCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCGTTCGCTGTCGGCCAGCTTCAGTGCGTgg                                                                                                                  <  7:7458/149‑1 (MQ=255)
                               cGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTgg                                                                                                                  <  6:101530/149‑1 (MQ=255)
                                     cGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAgc                                                                                                            <  3:292801/149‑1 (MQ=255)
                                      gggCACAAAGCTCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAgcg                                                                                                           <  8:164713/149‑1 (MQ=255)
                                        gCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGc                                                                                                                                         <  5:102596/117‑1 (MQ=255)
                                        gCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGc                                                                                                                                         >  6:102596/1‑117 (MQ=255)
                                          acaAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATg                                                                                                       >  3:118232/1‑149 (MQ=255)
                                           caAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGc                                                                                                      <  8:41447/149‑1 (MQ=255)
                                            aaaGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCg                                                                                                     >  6:100230/1‑149 (MQ=255)
                                                 tCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAAGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAg                                                                                                <  4:6165/149‑1 (MQ=255)
                                                     gTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGctct                                                                                            >  1:57730/1‑149 (MQ=255)
                                                         aaCGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGACGCGCTCggc                                                                                        >  6:251140/1‑147 (MQ=255)
                                                         aaCGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCtctc                                                                                           >  3:169682/1‑146 (MQ=255)
                                                         aaCGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCtctc                                                                                           <  4:169682/146‑1 (MQ=255)
                                                          aCGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCa                                                                                       >  1:89888/1‑149 (MQ=255)
                                                                cAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATgcagc                                                                                 <  2:276412/149‑1 (MQ=255)
                                                                            ggaggaTTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCccg                                                                     >  2:253035/1‑149 (MQ=255)
                                                                             gaggaTTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTAGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCccg                                                                     >  1:255978/1‑149 (MQ=255)
                                                                              aggaTTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCccgcc                                                                   >  7:95894/1‑149 (MQ=255)
                                                                                       aTGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAacg                                                          >  2:57488/1‑149 (MQ=255)
                                                                                        tGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAacg                                                          >  7:302885/1‑148 (MQ=255)
                                                                                           ccACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTgg                                                                                                                  >  8:136724/1‑89 (MQ=255)
                                                                                           ccACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTgg                                                                                                                  <  7:136724/89‑1 (MQ=255)
                                                                                            cACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAacgcacgc                                                     >  1:60268/1‑149 (MQ=255)
                                                                                               ggTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAacgcac                                                       <  1:24827/144‑1 (MQ=255)
                                                                                               ggTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAacgcac                                                       >  2:24827/1‑144 (MQ=255)
                                                                                                    ccTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACt                                             <  2:57730/149‑1 (MQ=255)
                                                                                                     cTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACtt                                            >  4:98111/1‑149 (MQ=255)
                                                                                                        cccGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTc                                          <  1:57488/148‑1 (MQ=255)
                                                                                                            gCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGa                                     >  6:345606/1‑149 (MQ=255)
                                                                                                              gATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATg                                   <  2:89888/149‑1 (MQ=255)
                                                                                                              gATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATg                                   <  2:60268/149‑1 (MQ=255)
                                                                                                                   caTCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCgtg                              <  8:302885/149‑1 (MQ=255)
                                                                                                                           gTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGtggtggtgg                       <  4:77467/149‑1 (MQ=255)
                                                                                                                           gTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGtggtggtgg                       <  1:253035/149‑1 (MQ=255)
                                                                                                                                 gAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGa                 >  4:264904/1‑149 (MQ=255)
                                                                                                                                      gTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATg            >  2:130102/1‑149 (MQ=255)
                                                                                                                                           aCGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGc       >  6:271333/1‑149 (MQ=255)
                                                                                                                                             gAAGGGCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGCCt     <  8:126594/149‑1 (MQ=255)
                                                                                                                                               aGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGCCtt    >  4:180635/1‑148 (MQ=255)
                                                                                                                                                ggTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGCCTTCa  <  3:98111/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TAGGCGGCGATCCCGGCGATCAAGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGTCATCGAACGGTGGCACGAAGGTCGCGCCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGCCTTCA  >  NC_002947/3904099‑3904390

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: