Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,560,540 A→G 9.2% intergenic (‑44/+258) 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,5400AG 9.2% 148.3 / 5.3 54intergenic (‑44/+258)hisZ/hflCATP phosphoribosyltransferase regulatory subunit/integral membrane ATP‑dependent zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base A (28/21);  new base G (2/3);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

AGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAG  >  NC_002947/5560396‑5560684
                                                                                                                                                |                                                                                                                                                 
aGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGt                                                                                                                                               >  6:211006/1‑149 (MQ=255)
   tCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGACGGTTTGCCTGCTGGCAAGCCgcccg                                                                                                                                            <  1:157421/149‑5 (MQ=255)
                      cgATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGca                                                                                                                          >  2:126964/1‑148 (MQ=255)
                          ctcGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCg                                                                                                                                    <  5:97269/134‑1 (MQ=255)
                          ctcGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCg                                                                                                                                    >  6:97269/1‑134 (MQ=255)
                           tcGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCAcagcag                                                                                                                    >  3:73003/1‑149 (MQ=255)
                           tcGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGACAGCGGTCTACCGTTGCCATTACGCCATTTCCACTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCAcagcag                                                                                                                    >  5:210564/1‑149 (MQ=255)
                              aTGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATcgcg                                                                                                                            <  8:230249/138‑1 (MQ=255)
                              aTGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATcgcg                                                                                                                            >  7:230249/1‑138 (MQ=255)
                                  gcgcAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGgcag                                                                                                             <  1:126964/149‑1 (MQ=255)
                                         cTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGAt                                                                                                      >  5:141225/1‑149 (MQ=255)
                                           cAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGaaaa                                                                                                                                        <  8:196406/113‑1 (MQ=255)
                                           cAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGaaaa                                                                                                                                        >  7:196406/1‑113 (MQ=255)
                                              gTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAg                                                                                                 >  1:122881/1‑149 (MQ=255)
                                                 gCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAgcg                                                                                               <  4:165036/148‑1 (MQ=255)
                                                        ttCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCAcc                                                                                       >  6:47636/1‑149 (MQ=255)
                                                         tCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGgca                                                                                                              <  7:157673/125‑1 (MQ=255)
                                                         tCCTCGATGCAATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGgca                                                                                                              >  8:157673/1‑125 (MQ=255)
                                                           cTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCa                                                                                                                                             >  6:38206/1‑92 (MQ=255)
                                                           cTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCa                                                                                                                                             <  5:38206/92‑1 (MQ=255)
                                                           cTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTa                                                                                    >  8:261284/1‑149 (MQ=255)
                                                              gATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTAccc                                                                                 <  7:261284/149‑1 (MQ=255)
                                                              gATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTAccc                                                                                 <  5:211006/149‑1 (MQ=255)
                                                                    aTCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACAcc                                                                          >  8:133165/1‑149 (MQ=255)
                                                                    aTCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACAcc                                                                          >  1:88859/1‑149 (MQ=255)
                                                                           agcagcCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAAt                                                                   <  7:133165/149‑1 (MQ=255)
                                                                              agccagcGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGc                                                                <  1:93195/149‑1 (MQ=255)
                                                                              agccagcGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGc                                                                <  6:141225/149‑1 (MQ=255)
                                                                              agccagcGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGc                                                                <  5:47636/149‑1 (MQ=255)
                                                                               gccagcGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGcc                                                               >  1:9259/1‑149 (MQ=255)
                                                                                  agcGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAAtt                                                                  <  6:192459/143‑1 (MQ=255)
                                                                                  agcGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAAtt                                                                  >  5:192459/1‑143 (MQ=255)
                                                                                   gcGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGc                                                           <  5:286362/149‑1 (MQ=255)
                                                                                   gcGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGc                                                           <  5:96984/149‑1 (MQ=255)
                                                                                          aCCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGa                                                    >  5:194194/1‑149 (MQ=255)
                                                                                            cGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGAc                                                   >  2:167929/1‑148 (MQ=255)
                                                                                                   aTTACGCCATTTCCCCTCTGGTCCGGACGGTTTGCCTGCTGGCAAGCCgc                                                                                                                                               >  5:99079/1‑49 (MQ=17)
                                                                                                   aTTACGCCATTTCCCCTCTGGTCCGGACGGTTTGCCTGCTGGCAAGCCgc                                                                                                                                               <  6:99079/50‑2 (MQ=17)
                                                                                                       cGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACAGTCAAc                                       >  5:190626/1‑149 (MQ=255)
                                                                                                          cATTTCCCCTCTTGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGc                                                       <  1:289850/130‑1 (MQ=255)
                                                                                                          cATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGCTGAGCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCg                                     >  1:257021/1‑149 (MQ=255)
                                                                                                          cATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGc                                                       >  2:289850/1‑130 (MQ=255)
                                                                                                                cccTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCTCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGa                              <  6:210564/149‑1 (MQ=255)
                                                                                                                 ccTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACCGCAGCGGGCAGCACTGCTCCCAGCGCCGCCGCCGTAGCCGACACCCTCAAATTGCCTTGCGAGCTGCCCAAACCGTCACCCGGCGCTGac                             >  7:280854/1‑148 (MQ=255)
                                                                                                                   tctGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGc                            >  3:214656/1‑148 (MQ=255)
                                                                                                                        tCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATc                       >  4:248797/1‑148 (MQ=255)
                                                                                                                          cGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGa                                                         >  8:261480/1‑112 (MQ=255)
                                                                                                                          cGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGa                                                         <  7:261480/112‑1 (MQ=255)
                                                                                                                            ggcggcTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCt                  >  4:54030/1‑149 (MQ=255)
                                                                                                                                gcTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCATCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCa              <  3:154731/149‑1 (MQ=255)
                                                                                                                                     cctgctggCAAGCCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGAc           >  3:74959/7‑147 (MQ=255)
                                                                                                                                     cctgctggCAAGCCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGAc           <  4:74959/141‑1 (MQ=255)
                                                                                                                                     ccAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGc         <  4:73003/149‑1 (MQ=255)
                                                                                                                                            gCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAg  >  2:25532/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGACGCAAAAAAG  >  NC_002947/5560396‑5560684

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: