Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,335 G→A 10.2% P404P (CCG→CCA guaA → glutamine‑hydrolyzing GMP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,178,3350GA10.2% 101.2 / 2.7 39P404P (CCG→CCAguaAglutamine‑hydrolyzing GMP synthase
Reads supporting (aligned to +/- strand):  ref base G (14/21);  new base A (2/2);  total (16/23)
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.94e-01

GTCGCACCACAACGTAGGCGGCCTGCCGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGCGTTCGTGGTGTTCCAGC  >  NC_002947/1178191‑1178480
                                                                                                                                                |                                                                                                                                                 
gTCGCACCACAACGTAGGCGGCCTGCCGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTcccgggccc                                                                                                                                              <  2:803201/150‑1 (MQ=255)
  cGCACCACAACGTAGGCGGCCTGCCGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTcccgggcccgg                                                                                                                                            <  1:1103420/150‑1 (MQ=255)
          aaCGTAGGCGGCCTGCCGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGgc                                                                                                                                    >  2:137551/1‑150 (MQ=255)
              gagtgtgcCTGCCGGAGGAGATGAACCTCCAGCTAGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCACGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGgcgtgc                                                                                                                                <  2:790086/144‑1 (MQ=255)
                      cTGCCGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCAGGCCTGGGCGTGCGTATCCTc                                                                                                                        >  1:428518/1‑150 (MQ=255)
                         ccGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGgcgtgcgt                                                                                                                              >  1:203634/1‑141 (MQ=255)
                         ccGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGgcgtgcgt                                                                                                                              <  2:203634/141‑1 (MQ=255)
                         ccGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCAGGcccgggcccgggaacgg                                                                                                                                >  2:582452/1‑134 (MQ=255)
                         ccGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCAGGcccgggcccgggaacgg                                                                                                                                <  1:582452/139‑6 (MQ=255)
                             aggagATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGgtgaag                                                                                                                 <  2:76121/150‑1 (MQ=255)
                                     aaCCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGgtgaagtgccgttt                                                                                                         <  1:588094/150‑7 (MQ=255)
                                       ccTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGa                                                                                                       >  1:244988/1‑150 (MQ=255)
                                       ccTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGa                                                                                                       <  2:325978/150‑1 (MQ=255)
                                              cTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGcc                                                                                                >  2:392141/1‑150 (MQ=255)
                                                       ccGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTg                                                                                       <  2:729217/150‑1 (MQ=255)
                                                               tGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCgcgc                                                                               >  2:1049433/1‑150 (MQ=255)
                                                                   cTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAgg                                                                                                        >  1:1057202/1‑121 (MQ=255)
                                                                   cTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAgg                                                                                                        <  2:1057202/121‑1 (MQ=255)
                                                                      ttCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACcac                                                                        >  1:489183/1‑150 (MQ=255)
                                                                            gACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCAGGcccgggcccgggaac                                                                                                                                  <  2:666173/86‑4 (MQ=255)
                                                                            gACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCAGGcccgggcccgggaac                                                                                                                                  >  1:666173/1‑83 (MQ=255)
                                                                               gAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCg                                                                                               <  1:353403/118‑1 (MQ=255)
                                                                               gAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCg                                                                                               >  2:353403/1‑118 (MQ=255)
                                                                                     cgcAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAagag                                                         >  1:7984/1‑150 (MQ=255)
                                                                                       cAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTc                                                                                                                        >  2:1043992/1‑85 (MQ=255)
                                                                                       cAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTc                                                                                                                        <  1:1043992/85‑1 (MQ=255)
                                                                                       cAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCt                                                       <  2:179321/150‑1 (MQ=255)
                                                                                         aGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTgc                                                     <  1:650536/150‑1 (MQ=255)
                                                                                             cGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCaa                                                 <  2:7984/150‑1 (MQ=255)
                                                                                                ccTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCAGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGc                                              <  2:428518/150‑1 (MQ=255)
                                                                                                     aGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGTAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCGTCGAAGAGCTGCGCAAGGCCGAc                                          <  2:554807/149‑1 (MQ=255)
                                                                                                     aGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGAc                                          >  1:554807/1‑149 (MQ=255)
                                                                                                        tGGGCATGTGATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGt                                      <  1:326153/150‑1 (MQ=255)
                                                                                                                             tCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGgtgaagtg                                                                                                               <  1:1071259/56‑1 (MQ=255)
                                                                                                                             tCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGgtgaagtg                                                                                                               >  2:1071259/1‑56 (MQ=255)
                                                                                                                              cTACCCCCATCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGCGtt                <  1:46854/150‑1 (MQ=255)
                                                                                                                                   gCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGCGTTCgtgg           <  1:137551/150‑1 (MQ=255)
                                                                                                                                      aCCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGCGTTCGTGgtgt        <  2:489183/150‑1 (MQ=255)
                                                                                                                                       cccGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGc                   >  1:893433/1‑138 (MQ=255)
                                                                                                                                       cccGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGc                   <  2:893433/138‑1 (MQ=255)
                                                                                                                                            tCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAATGCCGACTGGTACCTCAAGACCATCCAGGCGTTCGTGGTGTTCCAGc  <  2:244988/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTCGCACCACAACGTAGGCGGCCTGCCGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGCGTTCGTGGTGTTCCAGC  >  NC_002947/1178191‑1178480

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: