Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 852,159 A→T 6.9% Q313L (CAG→CTG)  hemA → glutamyl‑tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947852,1590AT 6.9% 175.2 / 2.7 58Q313L (CAG→CTG) hemAglutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base A (25/29);  new base T (2/2);  total (27/31)
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.38e-01

AAACAGCGCCGACACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAA  >  NC_002947/852014‑852300
                                                                                                                                                 |                                                                                                                                             
aaaCAGCGCCGACACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAggg                                                                                                                                            >  2:115555/1‑149 (MQ=255)
  aCAGCGCCGACACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGgcg                                                                                                                                          >  8:230492/1‑149 (MQ=255)
    aGCGCCGACACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGcgcg                                                                                                                                        <  7:267483/149‑1 (MQ=255)
    aGCGCCGACACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGcgcg                                                                                                                                        <  1:331104/149‑1 (MQ=255)
     gcgcCGACACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCgg                                                                                                                                       <  1:318153/149‑1 (MQ=255)
                   ccATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGaa                                                                                                                         <  3:327438/149‑1 (MQ=255)
                   ccATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGaa                                                                                                                         <  1:115555/149‑1 (MQ=255)
                    cATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCg                                                                                                                           >  5:89262/1‑146 (MQ=255)
                    cATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGaa                                                                                                                         <  6:89262/148‑1 (MQ=255)
                     aTCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGaa                                                                                                                         >  2:328285/1‑147 (MQ=255)
                          cATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTg                                                                                                                   <  8:331868/148‑1 (MQ=255)
                          cATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTg                                                                                                                   <  1:328285/148‑1 (MQ=255)
                          cATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTg                                                                                                                   <  7:276118/148‑1 (MQ=255)
                           aTGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCTGGGCCGCGCCCTggcggc                                                                                                                             <  2:248916/137‑1 (MQ=255)
                           aTGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCCAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCTGGGCCGCGCCCTggcggc                                                                                                                             >  1:248916/1‑137 (MQ=255)
                                  aCATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGATCTGGTATACGTg                                                                                                          >  3:221603/1‑149 (MQ=255)
                                  aCATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTg                                                                                                          >  5:17077/1‑149 (MQ=255)
                                        cTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGcc                                                                                                    >  5:63284/1‑149 (MQ=255)
                                          gTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCg                                                                                                   >  2:3196/1‑148 (MQ=255)
                                           tACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGgc                                                                                                                                           >  4:297167/1‑107 (MQ=255)
                                           tACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGgc                                                                                                                                           <  3:297167/107‑1 (MQ=255)
                                                    atgtCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCAtg                                                                                        <  4:230484/146‑1 (MQ=255)
                                                     tatCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGcc                                                                                                                            >  7:308120/1‑112 (MQ=255)
                                                     tatCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGcc                                                                                                                            <  8:308120/112‑1 (MQ=255)
                                                       tCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCAtgctg                                                                                     >  3:55254/1‑149 (MQ=255)
                                                         gAAACAGAAGTCGGCGAGCTGGGTGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGctgcg                                                                                   >  2:306269/1‑149 (MQ=255)
                                                            aCCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCa                                                                                          >  6:82796/1‑139 (MQ=255)
                                                            aCCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCa                                                                                          <  5:82796/139‑1 (MQ=255)
                                                                       cGAGCTGGATGATGTCACCCTGTACCCCGTCGACGACCTCCACGATGTAGTGGCGTAAAACCTCAAGCGCCGCCAGGGCGCGGCCCAGGCGGCCTAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTgg                                                                     <  7:292168/149‑1 (MQ=255)
                                                                        gAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGc                                                                    <  6:17077/149‑1 (MQ=255)
                                                                          gCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTg                                                                  >  5:146689/1‑149 (MQ=255)
                                                                              ggtgatgTCTACCTGTACACCGTCGACGAACTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCa                                                                <  1:306269/145‑1 (MQ=255)
                                                                              gatgatGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCTGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACGGGCTGCa                                                                <  7:230492/147‑1 (MQ=255)
                                                                              gatgatGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGt                                                              >  4:36840/1‑149 (MQ=255)
                                                                                 gatgTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTg                                                                  >  5:62770/1‑142 (MQ=255)
                                                                                 gatgTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTg                                                                  <  6:62770/142‑1 (MQ=255)
                                                                                   tgTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGAc                                                          <  6:63284/148‑1 (MQ=255)
                                                                                       tACCTGTACACCGTCGACGACCCGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCt                                                     >  2:96254/1‑149 (MQ=255)
                                                                                             tACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAgg                                                 <  4:55254/147‑1 (MQ=255)
                                                                                                      gacgacCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAg                                                               <  2:226218/124‑1 (MQ=255)
                                                                                                      gacgacCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAg                                                               >  1:226218/1‑124 (MQ=255)
                                                                                                         gacCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGc                                    >  2:208417/1‑148 (MQ=255)
                                                                                                            cTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAg                                >  1:313359/1‑149 (MQ=255)
                                                                                                               cACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGt                               >  1:155407/1‑147 (MQ=255)
                                                                                                                aCGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAg                            <  1:208417/149‑1 (MQ=255)
                                                                                                                  gATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCT‑AGGGCCTATCGCCAGCAAAGTGAgcgc                         >  8:125305/1‑149 (MQ=255)
                                                                                                                        gTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTcg                    <  1:3196/149‑1 (MQ=255)
                                                                                                                           gCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGcga                 >  5:103845/1‑149 (MQ=255)
                                                                                                                           gCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGcga                 <  2:155407/149‑1 (MQ=255)
                                                                                                                             ggAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAggg                                                <  6:282260/116‑1 (MQ=255)
                                                                                                                             ggAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAggg                                                >  5:282260/1‑116 (MQ=255)
                                                                                                                               aaaaCCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACgaa             <  3:36840/149‑1 (MQ=255)
                                                                                                                                aaaCCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACgaag            <  6:103845/149‑1 (MQ=255)
                                                                                                                                     tCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCa                                                   <  7:173783/105‑1 (MQ=255)
                                                                                                                                     tCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCa                                                   >  8:173783/1‑105 (MQ=255)
                                                                                                                                         gagCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCT‑AGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCaaaa  <  7:125305/149‑1 (MQ=255)
                                                                                                                                            ccgccTGGGCCGCGCCCTGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGcgac                <  7:62327/133‑1 (MQ=255)
                                                                                                                                            ccgccTGGGCCGCGCCCTGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGcgac                >  8:62327/1‑133 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
AAACAGCGCCGACACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAA  >  NC_002947/852014‑852300

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