Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,559,994 T→A 11.8% Q185L (CAG→CTG)  hisZ ← ATP phosphoribosyltransferase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,559,9940TA11.8% 89.0 / 4.7 34Q185L (CAG→CTG) hisZATP phosphoribosyltransferase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base T (15/15);  new base A (2/2);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

GGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTA  >  NC_002947/5559860‑5560137
                                                                                                                                      |                                                                                                                                               
gggCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCgcgc                                                                                                                                  >  2:198991/1‑150 (MQ=255)
  gCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCg                                                                                                                                <  2:158205/150‑1 (MQ=255)
    ttCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATAGCGTCGAACAGCTGCTGTTCGACCGCGCCGGa                                                                                                                              >  2:516537/1‑150 (MQ=255)
     tCAGCCAGTACTTCCCGACCACCACCCAGCTCGACCAGGGCTCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGAc                                                                                                                             <  1:274214/150‑1 (MQ=255)
                ttCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCAAGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCggcggc                                                                                                                  >  2:474323/1‑150 (MQ=255)
                         accacACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAg                                                                                                         >  1:363964/1‑150 (MQ=255)
                          ccacACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAgg                                                                                                        >  1:260142/1‑150 (MQ=255)
                           cacacaGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGc                                                                                                       >  1:420870/1‑150 (MQ=255)
                                              gcgcAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATg                                                                                    >  1:543572/1‑150 (MQ=255)
                                                 cagcaTGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGACTGCACTTCATCGACAGCCTTGCGTTGAATGGAGTGGAACAGCTGCTGTTAAACAGCACTTGCTGGACCGGCGTGACTGGCCAGGCCGCGGTAGATGCCGACATGGcc                                                                                 <  1:322995/150‑1 (MQ=255)
                                                    caTGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAg                                                                              <  2:272273/150‑1 (MQ=255)
                                                    caTGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAg                                                                              <  1:821089/150‑1 (MQ=255)
                                                         tgcCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCa                                                                                                          >  1:492492/1‑117 (MQ=255)
                                                         tgcCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCa                                                                                                          <  2:492492/117‑1 (MQ=255)
                                                          gcCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAt                                                                        <  2:363964/150‑1 (MQ=255)
                                                               ggTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGaa                                                                   <  2:673576/150‑1 (MQ=255)
                                                                      cGGCAGGTCGGCCGGCAGCGCCTGCACTTCATCGACAGCCGTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCGGACATGGCCCAGGTCCATGTGAACATCCgg                                                            <  1:111561/150‑1 (MQ=255)
                                                                        gCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtgt                                                                      >  1:373991/1‑138 (MQ=255)
                                                                        gCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtgt                                                                      <  2:373991/138‑1 (MQ=255)
                                                                                        cgcCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAg                                          <  1:590963/150‑1 (MQ=255)
                                                                                          ccTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCg                                        >  1:423467/1‑150 (MQ=255)
                                                                                                 ttCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAg                                 <  2:260142/150‑1 (MQ=255)
                                                                                                   cATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCAAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCa                               <  1:198991/150‑1 (MQ=255)
                                                                                                       gACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAg                           >  2:434340/1‑150 (MQ=255)
                                                                                                                  gTTGCACGGCGTTGCACAGCTGGTGTTTGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACtt                <  2:266249/150‑1 (MQ=255)
                                                                                                                     gCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGa             <  1:434340/150‑1 (MQ=255)
                                                                                                                       aTGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGCCAGACCGGCGGCCCGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACa           >  2:223126/1‑150 (MQ=255)
                                                                                                                         ggCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATc         >  2:111915/1‑150 (MQ=255)
                                                                                                                            gTCGAACAGCAGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCtga                      >  2:847864/1‑134 (MQ=255)
                                                                                                                            gTCGAACAGCAGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCtga                      <  1:847864/134‑1 (MQ=255)
                                                                                                                             tCGAACAGCAGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtgt                                                                      >  2:535201/1‑85 (MQ=255)
                                                                                                                             tCGAACAGCAGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtgt                                                                      <  1:535201/85‑1 (MQ=255)
                                                                                                                                aaCAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTa  >  1:815589/1‑150 (MQ=255)
                                                                                                                                aaCAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGACAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTa  >  1:57333/1‑150 (MQ=255)
                                                                                                                                      |                                                                                                                                               
GGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTA  >  NC_002947/5559860‑5560137

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