Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,714,353 G→A 9.7% R36H (CGT→CAT)  hisI → phosphoribosyl‑AMP cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,714,3530GA 9.7% 107.7 / 2.7 41R36H (CGT→CAT) hisIphosphoribosyl‑AMP cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base G (17/20);  new base A (2/2);  total (19/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.68e-01

AAAGGGCCCGGCACAGGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGAGTTGCGCCTGGACTGTGACGCCGACGTGATCATC  >  NC_002947/5714206‑5714498
                                                                                                                                                   |                                                                                                                                                 
aaaGGGCCCGGCACAGGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGt                                                                                                                                                  <  4:75643/149‑1 (MQ=255)
  aGGGCCCGGCACAGGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGa                                                                                                                                                <  6:189689/149‑1 (MQ=255)
   gggCCCGGCACAGGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGTTGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGaa                                                                                                                                               <  2:38766/149‑1 (MQ=255)
        cGGCACAGGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCt                                                                                                                                          >  4:137548/1‑149 (MQ=255)
          gCACAGGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTgg                                                                                                                                        >  2:225994/1‑149 (MQ=255)
                      gCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGAccgccg                                                                                                                            >  3:71421/1‑149 (MQ=255)
                           cTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGCCCACAAGACCGGCCGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGag                                                                                                                       >  4:144144/1‑149 (MQ=255)
                                 aGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGATCCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGt                                                                                                                 <  1:225994/149‑1 (MQ=255)
                                 aGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGt                                                                                                                 <  5:127161/149‑1 (MQ=255)
                                     gacgaTGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCATGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCa                                                                                                             >  6:97875/1‑149 (MQ=255)
                                      acgaTGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTgg                                                                                                                                        <  5:262798/121‑1 (MQ=255)
                                      acgaTGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTgg                                                                                                                                        >  6:262798/1‑121 (MQ=255)
                                        gaTGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCt                                                                                                          <  3:78729/149‑1 (MQ=255)
                                          tGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTa                                                                                                         >  4:11654/1‑148 (MQ=255)
                                                aCTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGtcg                                                                                                  >  5:99211/1‑149 (MQ=255)
                                                       ggACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGtcgcgttcgc                                                                                           >  6:14497/1‑149 (MQ=255)
                                                         aCGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGtcgcgttcgcgt                                                                                         <  3:11654/149‑1 (MQ=255)
                                                              aTCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCATGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCaa                                                                                    <  5:97875/149‑1 (MQ=255)
                                                                    tGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCtgtg                                                                              <  5:14497/149‑1 (MQ=255)
                                                                      gAACAGCGATGGCCTGGTACCGGCGATCGCCCAGTACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCATGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGGTCGCGTGGCAAGCTGTGGc                                                                            >  5:176082/1‑149 (MQ=255)
                                                                      gAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGAccgccgc                                                                                                                           <  2:307531/102‑1 (MQ=255)
                                                                      gAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGAccgccgc                                                                                                                           >  1:307531/1‑102 (MQ=255)
                                                                            cGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAg                                                                      >  3:84737/1‑149 (MQ=255)
                                                                                      gTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTc                                                            >  3:278215/1‑149 (MQ=255)
                                                                                           ggCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCggg                                                        >  1:15581/1‑148 (MQ=255)
                                                                                             cGATCGCCCAGGACCACTAGTCCGGACGCGGGCTGCTGTTGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCAt                                                     <  6:52163/149‑1 (MQ=255)
                                                                                                    ccAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCAt                                                     >  7:282705/1‑142 (MQ=255)
                                                                                                    ccAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCAt                                                     <  8:282705/142‑1 (MQ=255)
                                                                                                      aGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCaaaaa                                            >  2:178924/1‑149 (MQ=255)
                                                                                                            acaAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCAt                                      <  4:71421/149‑1 (MQ=255)
                                                                                                              aaGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGa                                    <  3:144144/149‑1 (MQ=255)
                                                                                                                gACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGAGt                                  >  2:2044/1‑149 (MQ=255)
                                                                                                                   cGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGAGTTgc                               <  1:178924/149‑1 (MQ=255)
                                                                                                                          gTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGAGTTGCGCCTGGa                        <  4:278215/149‑1 (MQ=255)
                                                                                                                               gatgatGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGAGTTGCGCCTGGACtgtg                   <  2:15581/149‑1 (MQ=255)
                                                                                                                                atgatgGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGAGTTGCGCCTGGACtgtg                   <  6:99211/148‑1 (MQ=255)
                                                                                                                                       ccgggATGAACCATGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGGTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGAGTTGCGCCTGGACTGTGACGCCGAc           <  6:176082/146‑1 (MQ=255)
                                                                                                                                       ccTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGAGTTGCGCCTGGACTGTGACGCCGAc           >  1:128838/1‑149 (MQ=255)
                                                                                                                                          ggATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGAGTTGCGCCTGGACTGTGACGCCGACGt         >  6:306348/1‑148 (MQ=255)
                                                                                                                                           gATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGAGTTGCGCCTGGACTGTGACGCCGACGt         <  1:2044/147‑1 (MQ=255)
                                                                                                                                                aCCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGAGTTGCGCCTGGACTGTGACGCCGACGTGatcatc  <  5:306348/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AAAGGGCCCGGCACAGGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACAGCGATGGCCTGGTACCGGCGATCGCCCAGGACCACAAGACCGGACGCGTGCTGATGATGGCCTGGATGAACCGTGAATCGCTGGCCCTGACCGCCGCCGAGCAGCGTGCCATCTACTGGTCGCGTTCGCGTGGCAAGCTGTGGCGTAAAGGCGAGGAGTCCGGGCATGTGCAAAAACTGCATGAGTTGCGCCTGGACTGTGACGCCGACGTGATCATC  >  NC_002947/5714206‑5714498

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