Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,560,671 T→C 9.9% intergenic (‑124/+178) hisZ ← / ← hflC ATP phosphoribosyltransferase regulatory subunit/integral membrane ATP‑dependent zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,560,6710TC 9.9% 110.5 / 3.0 44intergenic (‑124/+178)hisZ/hflCATP phosphoribosyltransferase regulatory subunit/integral membrane ATP‑dependent zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base C (2/2);  total (22/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 = 3.42e-01

CGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCACCCCGCCTGGTGCATTAGCTGCCAGGCGG  >  NC_002947/5560531‑5560816
                                                                                                                                            |                                                                                                                                                 
cGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGa                                                                                                                                          <  2:1049823/150‑1 (MQ=255)
   tCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCt                                                                                                                                       <  1:307325/150‑1 (MQ=255)
      cccGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGAc                                                                                                                                    <  2:407917/149‑1 (MQ=255)
         gTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGc                                                                                                                                                 >  2:507333/1‑134 (MQ=255)
         gTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGc                                                                                                                                                 <  1:507333/134‑1 (MQ=255)
                ttACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCa                                                                                                                          >  1:622212/1‑150 (MQ=255)
                                ctGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCaa                                                                                                          <  2:1075684/150‑1 (MQ=255)
                                       gggCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCt                                                                                                   >  2:769554/1‑150 (MQ=255)
                                             cTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGa                                                                                             >  1:321638/1‑150 (MQ=255)
                                              ttGCCAGCAGGCAAACCGTCAGAAAAGCCGATATCGCGCACAGCAGCGGGCAGCACTGATACCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGAc                                                                                            >  1:281410/1‑150 (MQ=255)
                                                       ggCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGCCACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAg                                                                                   <  1:625195/150‑1 (MQ=255)
                                                        gCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCAAAGATCCAAGCGCAGACACACAACCCGGCACCAGAAAAGAAGCGTGCGGTGTGAGCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGc                                                                                  >  2:678965/1‑150 (MQ=255)
                                                          aaaCCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCg                                                                                >  2:845640/1‑150 (MQ=255)
                                                          aaaCCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCg                                                                                >  2:303780/1‑150 (MQ=255)
                                                                        gCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCTAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGc                                                                  >  1:463508/1‑150 (MQ=255)
                                                                               tCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGGGCGGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACACAAAAAAGCCGGGAAATTTCCCGGCTGCCtca                                                           >  1:356284/1‑150 (MQ=255)
                                                                                            cGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCtt                                              <  1:845640/150‑1 (MQ=255)
                                                                                               gcagcaCTGATCCCAGAGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTAttc                                           >  2:651747/1‑150 (MQ=255)
                                                                                                cagcaCTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTttca                                          <  2:836191/150‑1 (MQ=255)
                                                                                                cagcaCTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTttca                                          >  2:344136/1‑150 (MQ=255)
                                                                                                  gcaCTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTttcatt                                        <  1:940571/150‑1 (MQ=255)
                                                                                                   caCTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTttcattc                                       <  2:946323/150‑1 (MQ=255)
                                                                                                   caCTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTttcattc                                       >  1:571182/1‑150 (MQ=255)
                                                                                                           ccAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCa                               <  1:651747/150‑1 (MQ=255)
                                                                                                            cAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCAc                              >  2:343272/1‑150 (MQ=255)
                                                                                                            cAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCAc                              >  1:883200/1‑150 (MQ=255)
                                                                                                            cAGCGCAGCCACCGTACCCGACACCCTCAAATCGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCtt                                                                                                  >  2:929882/1‑82 (MQ=255)
                                                                                                            cAGCGCAGCCACCGTACCCGACACCCTCAAATCGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCtt                                                                                                  <  1:929882/82‑1 (MQ=255)
                                                                                                            cAGCGCAGCCACCGTACCCGACACCCTCAAATCGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCtt                                                                                                  <  1:929874/82‑1 (MQ=255)
                                                                                                            cAGCGCAGCCACCGTACCCGACACCCTCAAATCGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCtt                                                                                                  >  2:929874/1‑82 (MQ=255)
                                                                                                               cgcAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCAcccc                           >  2:881911/1‑150 (MQ=255)
                                                                                                                   gCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCACCCCGCCt                       <  1:344136/150‑1 (MQ=255)
                                                                                                                    ccACCGTACCCGACACCCTCAACTTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCACCCCGCCTg                      <  1:654490/150‑1 (MQ=255)
                                                                                                                    ccACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGc                                                               >  2:459434/1‑109 (MQ=255)
                                                                                                                    ccACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGc                                                               <  1:459434/109‑1 (MQ=255)
                                                                                                                    ccACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGC‑ACAAAGCCTGCAAATGTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCACCCCGCCTgg                     <  2:427547/150‑1 (MQ=255)
                                                                                                                      aCCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAATCGGCGCTGAGCCAATCAACCTTGCAGACGCTAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCACCCCGCCTGGt                    <  2:571182/150‑1 (MQ=255)
                                                                                                                      aCCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCACCCCGCCTggg                    <  2:321638/150‑2 (MQ=255)
                                                                                                                        cGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCACCCCGCCTGGTGc                  >  2:192558/1‑150 (MQ=255)
                                                                                                                             ccGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATaca                                                    >  1:419684/1‑111 (MQ=255)
                                                                                                                             ccGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATaca                                                    <  2:419684/111‑1 (MQ=255)
                                                                                                                                    ccTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCACCCCGCCTGGTGCATTAGCTGCCAg      >  2:944064/1‑150 (MQ=255)
                                                                                                                                        aaaTTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCTAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCACCCCGCCTGGTGCATTAGCTGCCAGGCgg  <  2:463508/150‑1 (MQ=255)
                                                                                                                                        aaaTTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCACCCCGCCTGGTGCATTAGCTGCCAGGCgg  <  2:281410/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAGCCGGGAAATTTCCCGGCTGCCTCATCATACACCCCTTTTCATTCAGGATGCACCCCGCCTGGTGCATTAGCTGCCAGGCGG  >  NC_002947/5560531‑5560816

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