Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,675,431 G→A 8.7% A200T (GCT→ACT)  hom → homoserine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,675,4310GA 8.7% 128.9 / 3.0 46A200T (GCT→ACT) homhomoserine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base A (2/2);  total (24/22)
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.84e-01

CTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGAACTACGCCGAAGCCTTGGGTTACCGCATCAAGCACCTGGGCGT  >  NC_002947/1675283‑1675579
                                                                                                                                                    |                                                                                                                                                    
ctggctggCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACg                                                                                                                                                      <  4:240655/149‑1 (MQ=255)
   gctggcCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTg                                                                                                                                                   >  6:193613/1‑149 (MQ=255)
       gcCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCTACGTGGAAGGCATCGACGCTGCGca                                                                                                                                               >  6:281452/1‑149 (MQ=255)
       gcCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTgcgc                                                                                                                                                <  3:201042/148‑1 (MQ=255)
         cGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGcac                                                                                                                                              >  8:225219/1‑148 (MQ=255)
         cGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACACTGCGcac                                                                                                                                              >  2:280046/1‑148 (MQ=255)
         cGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACACTGCGcac                                                                                                                                              >  6:103543/1‑148 (MQ=255)
             atcatcAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCt                                                                                                                                         <  5:193613/149‑1 (MQ=255)
                atcaACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGAc                                                                                                                                      >  6:150627/1‑149 (MQ=255)
                   aaCGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGAt                                                                                                                                   <  2:270443/149‑1 (MQ=255)
                                     aTCCTCACCTAAATTGGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGAAGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGtt                                                                                                                 <  8:231556/149‑1 (MQ=255)
                                     aTCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTATGGCGTCCATCGCGt                                                                                                                  >  7:231556/1‑148 (MQ=255)
                                      tCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTc                                                                                                                <  7:225219/149‑1 (MQ=255)
                                        cTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCg                                                                                                               <  4:88083/148‑1 (MQ=255)
                                                         gaAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTTGGTATCCCGCTGCAGTTc                                                                                              >  2:96811/1‑148 (MQ=255)
                                                         gaAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTc                                                                                              >  8:176175/1‑148 (MQ=255)
                                                          aaGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGa                                                                                            <  5:150627/149‑1 (MQ=255)
                                                                  tGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGcc                                                                                     >  5:204467/1‑148 (MQ=255)
                                                                           ggACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAg                                                                           <  7:176175/149‑1 (MQ=255)
                                                                           ggACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGCTGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAg                                                                           <  1:75353/149‑1 (MQ=255)
                                                                                tGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATc                                                                      <  7:42839/149‑1 (MQ=255)
                                                                                   tggctgAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCAcc                                                                   >  8:269113/1‑149 (MQ=255)
                                                                                   tggctgAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCAcc                                                                   >  7:203505/1‑149 (MQ=255)
                                                                                         aaGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCt                                                              >  6:137575/1‑148 (MQ=255)
                                                                                          aGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGa                                                            <  5:137575/149‑1 (MQ=255)
                                                                                              cAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAAc                                                        <  7:269113/149‑1 (MQ=255)
                                                                                                ggCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACgg                                                      >  7:175888/1‑149 (MQ=255)
                                                                                                 gcgcTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGt                                                                         <  1:4796/129‑1 (MQ=255)
                                                                                                 gcgcTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGt                                                                         >  2:4796/1‑129 (MQ=255)
                                                                                                         cTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACgg                                                      <  4:166470/140‑1 (MQ=255)
                                                                                                         cTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGatcg                                               >  3:166470/1‑144 (MQ=255)
                                                                                                          tACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCt                                                              <  5:300659/131‑1 (MQ=255)
                                                                                                          tACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCt                                                              >  6:300659/1‑131 (MQ=255)
                                                                                                            cGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACACTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGaa                                            <  1:280046/147‑1 (MQ=255)
                                                                                                            cGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACACTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGaa                                            <  5:103543/147‑1 (MQ=255)
                                                                                                                   gCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGAACTACGCCg                                    >  4:222841/1‑148 (MQ=255)
                                                                                                                    ccgacccgacCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCTGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGAc                                                 <  4:262291/134‑1 (MQ=255)
                                                                                                                    ccgacccgacCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGAc                                                 >  3:262291/1‑134 (MQ=255)
                                                                                                                         ccgaccTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGAACTACGCCGAAGCCtt                             >  5:86769/1‑149 (MQ=255)
                                                                                                                           gaccTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTTGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGAACTACGCCGAAGCCTTgg                           <  1:96811/149‑1 (MQ=255)
                                                                                                                           gaccTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGAACTACGCCGAAGCCTTgg                           >  7:132729/1‑149 (MQ=255)
                                                                                                                           gaccTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGAACTACGCCGAAGCCTTgg                           >  7:271857/1‑149 (MQ=255)
                                                                                                                             ccTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGAACTACGCCGAAGCCTTGGGt                         <  3:222841/149‑1 (MQ=255)
                                                                                                                                tCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGAACTACGCCGAAGCCTTGGGTTAc                      >  1:282904/1‑149 (MQ=255)
                                                                                                                                          aGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGAACTACGCCGAAGCCTTGGGTTACCGCATCAAGc            <  6:86769/149‑1 (MQ=255)
                                                                                                                                              aTCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGa                                             <  4:88819/112‑1 (MQ=255)
                                                                                                                                              aTCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGa                                             >  3:88819/1‑112 (MQ=255)
                                                                                                                                                    gCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGAACTACGCCGAAGCCTTGGGTTACCGCATCAAGCACCTGGGCGt  <  8:271857/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCAGACGTGAACTACGCCGAAGCCTTGGGTTACCGCATCAAGCACCTGGGCGT  >  NC_002947/1675283‑1675579

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