Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,496,672 T→A 6.6% intergenic (+22/+30) PP_3987 → / ← PP_5651 site‑specific recombinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,496,6720TA 6.6% 182.8 / 2.7 61intergenic (+22/+30)PP_3987/PP_5651site‑specific recombinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base A (2/2);  total (27/34)
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.44e-01

GGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCTAGT  >  NC_002947/4496533‑4496819
                                                                                                                                           |                                                                                                                                                   
ggcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGAGCGTACCCCGCttc                                                                                                                                            <  1:32451/149‑2 (MQ=255)
  cggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                          <  8:101737/149‑1 (MQ=255)
  cggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                          <  4:122539/149‑1 (MQ=255)
          gAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAAtt                                                                                                                                  <  7:127570/149‑1 (MQ=255)
            aaTCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTc                                                                                                                                 <  1:306463/148‑1 (MQ=255)
             aTCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTcac                                                                                                                               <  3:92422/149‑1 (MQ=255)
                ggCGATTTGGGTCCAGCCAGGCGCTTATTGGGGCCCACTTGCAAGGCGATCCACCATACCGTAAATATACCTGGATGTGAAATCGTTCACCCCACACGATACGTAAATGGAACCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAAc                                                                                                                            <  2:230760/149‑1 (MQ=255)
                ggCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATCTAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACaa                                                                                                                             <  8:265461/148‑1 (MQ=255)
                   gATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGa                                                                                                                         >  5:321682/1‑149 (MQ=255)
                    aTTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGa                                                                                                                         >  7:278303/1‑148 (MQ=255)
                        gggTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTa                                                                                                                     <  6:31350/148‑1 (MQ=255)
                        gggTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTAt                                                                                                                    >  7:41297/1‑149 (MQ=255)
                           tCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTAtt                                                                                                                   >  1:41935/1‑147 (MQ=255)
                           tCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGa                                                                                                                 <  8:41297/149‑1 (MQ=255)
                               gCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCta                                                                                                             >  2:304848/1‑149 (MQ=255)
                                 aaGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCtata                                                                                                           <  2:216832/149‑1 (MQ=255)
                                          aTTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGAt                                                                                                                        >  5:97086/1‑127 (MQ=255)
                                          aTTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGAt                                                                                                                        <  6:97086/127‑1 (MQ=255)
                                             gggTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATc                                                                                                <  2:259876/148‑1 (MQ=255)
                                                tCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTAtt                                                                                                                   >  6:28238/1‑126 (MQ=255)
                                                tCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTAtt                                                                                                                   <  5:28238/126‑1 (MQ=255)
                                                 cacaCTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCt                                                                                           >  5:5790/1‑149 (MQ=255)
                                                  acacTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCt                                                                                           <  6:321682/148‑1 (MQ=255)
                                                    acTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAg                                                                                        >  4:273505/1‑149 (MQ=255)
                                                      tGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCt                                                                                      >  6:336525/1‑149 (MQ=255)
                                                      tGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCt                                                                                      >  7:180092/1‑149 (MQ=255)
                                                        aCAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCt                                                                                      <  2:195478/147‑1 (MQ=255)
                                                        aCAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTtg                                                                                    >  4:312033/1‑149 (MQ=255)
                                                         cAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTtgc                                                                                   <  3:177277/149‑1 (MQ=255)
                                                          aaGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTa                                                                                         <  5:134142/142‑1 (MQ=255)
                                                          aaGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTa                                                                                         >  6:134142/1‑142 (MQ=255)
                                                            gcgcATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTtgctgc                                                                                >  6:52907/1‑149 (MQ=255)
                                                                   acacaAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCa                                                                         >  6:226842/1‑149 (MQ=255)
                                                                           ggTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAAt                                                                                             >  4:213202/1‑121 (MQ=255)
                                                                           ggTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAAt                                                                                             <  3:213202/121‑1 (MQ=255)
                                                                               tataGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAAtcatca                                                              <  5:272525/148‑1 (MQ=255)
                                                                               tataGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAAtcatca                                                              >  8:217234/1‑148 (MQ=255)
                                                                                 taGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGTACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATa                                                            <  7:49786/148‑1 (MQ=255)
                                                                                 taGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAg                                                           >  1:186713/1‑149 (MQ=255)
                                                                                 taGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAg                                                           <  5:336525/149‑1 (MQ=255)
                                                                                      gAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTc                                                      <  2:41935/149‑1 (MQ=255)
                                                                                           ggTGAAATCGTACAGCAAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCgaga                                                  <  7:217234/148‑1 (MQ=255)
                                                                                                  tCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATAcc                                           >  6:89750/1‑148 (MQ=255)
                                                                                                  tCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCt                                          >  3:98017/1‑149 (MQ=255)
                                                                                                          ccAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGcc                                  >  2:75155/1‑149 (MQ=255)
                                                                                                             acacGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGa                               <  3:312033/149‑1 (MQ=255)
                                                                                                             acacGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGa                               <  3:273505/149‑1 (MQ=255)
                                                                                                              cacGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAg                              <  8:180092/149‑1 (MQ=255)
                                                                                                                cGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTa                            <  6:5790/149‑1 (MQ=255)
                                                                                                                 gATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTTATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGctct                                                                             <  6:145481/99‑1 (MQ=255)
                                                                                                                 gATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGctct                                                                             >  5:145481/1‑99 (MQ=255)
                                                                                                                      gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCa                                                                                               <  1:281588/76‑1 (MQ=255)
                                                                                                                      gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGGATCa                                                                                               >  2:281588/1‑76 (MQ=255)
                                                                                                                         aTTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTg                    <  5:226842/148‑1 (MQ=255)
                                                                                                                           tcaaaGCGGGGAGCGTACCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                     >  4:169973/4‑145 (MQ=255)
                                                                                                                           tcaaaGCGGGGAGCGTACCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                     <  3:169973/142‑1 (MQ=255)
                                                                                                                           tcaaaGCGGGGAGCGTACCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGAc                 >  8:96067/4‑149 (MQ=255)
                                                                                                                                  ggggTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCa          <  5:52907/149‑1 (MQ=255)
                                                                                                                                   gggTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTCCCACCCAATCATCAAAGGTTTCGAGAGAATACCACGCACGCCTGAGTAGCTCTTTGGACCTTGc            >  7:215317/1‑146 (MQ=255)
                                                                                                                                         gCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCTAg   >  2:263727/1‑149 (MQ=255)
                                                                                                                                          cTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCTAGt  <  5:89750/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCTAGT  >  NC_002947/4496533‑4496819

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: