Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,227,422 T→A 8.9% intergenic (‑21/‑426) PP_1964 ← / → holB dGMP/dTMP/dhydroxyCMP kinase/DNA polymerase III subunit delta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,227,4220TA 8.9% 120.9 / 3.9 45intergenic (‑21/‑426)PP_1964/holBdGMP/dTMP/dhydroxyCMP kinase/DNA polymerase III subunit delta'
Reads supporting (aligned to +/- strand):  ref base T (17/23);  new base A (2/2);  total (20/25)
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.49e-01

CTTGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGA  >  NC_002947/2227275‑2227563
                                                                                                                                                   |                                                                                                                                             
cTTGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAtt                                                                                                                                              >  1:244427/1‑149 (MQ=255)
 ttGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAttt                                                                                                                                             >  7:104738/1‑149 (MQ=255)
      aaCCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCggg                                                                                                                                         >  4:152509/1‑148 (MQ=255)
      aaCCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCggg                                                                                                                                         <  3:152509/148‑1 (MQ=255)
           gCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGt                                                                                                                                      >  2:126735/1‑146 (MQ=255)
              ttGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCggg                                                                                                                                         >  8:104790/1‑140 (MQ=255)
              ttGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCggg                                                                                                                                         <  7:104790/140‑1 (MQ=255)
                gAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAggg                                                                                                                               <  1:248885/148‑1 (MQ=255)
                    ggATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTa                                                                                                                          <  8:104738/149‑1 (MQ=255)
                     gATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATg                                                                                                                            >  3:212777/1‑146 (MQ=255)
                                   tatGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACa                                                                                                           <  6:109995/149‑1 (MQ=255)
                                     tGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACAtt                                                                                                         <  6:145873/149‑1 (MQ=255)
                                     tGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACAtt                                                                                                         >  5:53983/1‑149 (MQ=255)
                                      gTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTc                                                                                                        >  7:58204/1‑149 (MQ=255)
                                          cgACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAGTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAg                                                                                                    >  8:233520/1‑149 (MQ=255)
                                             ccTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTaa                                                                                                 <  1:126735/149‑1 (MQ=255)
                                              cTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCt                                                                                                       >  7:213347/1‑142 (MQ=255)
                                              cTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCt                                                                                                       <  8:213347/142‑1 (MQ=255)
                                                        tCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTataa                                                                                      >  4:175136/1‑149 (MQ=255)
                                                          cATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTataata                                                                                    <  3:175136/149‑1 (MQ=255)
                                                              aGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGc                                                                                 <  4:124609/148‑1 (MQ=255)
                                                                  tGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATg                                                                            >  6:185000/1‑149 (MQ=255)
                                                                        cGAAGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACa                                                                      <  4:212777/149‑1 (MQ=255)
                                                                          aCGGTGTCATTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACAtt                                                                    <  6:40188/149‑1 (MQ=255)
                                                                            ggTGTCTTTACCGGATCGAGCTCATGCGGCGAGGCCGATGATGGTTGGGATTTTTGGCTCCAGGCCCAGAATTTGGGGGGTTagag                                                                                                                                 >  7:169108/1‑86 (MQ=255)
                                                                            ggTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTagag                                                                                                                                 <  8:169108/86‑1 (MQ=255)
                                                                             gtgtCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAgg                                                                                               >  2:219618/1‑119 (MQ=255)
                                                                             gtgtCTTTACCGGATCGAGCTACAGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAgg                                                                                               <  1:219618/119‑1 (MQ=255)
                                                                                tCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCa                                                              <  2:244427/149‑1 (MQ=255)
                                                                                              aGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCt                                                                                                       <  5:231911/94‑1 (MQ=255)
                                                                                              aGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCt                                                                                                       >  6:231911/1‑94 (MQ=255)
                                                                                                      cggcgAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCAc                                         <  6:53983/148‑1 (MQ=255)
                                                                                                                       gTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGggcgg                       >  3:211316/1‑149 (MQ=255)
                                                                                                                          ggCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATATGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTa                                <  2:160468/137‑1 (MQ=255)
                                                                                                                          ggCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTa                                >  1:160468/1‑137 (MQ=255)
                                                                                                                             aTTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCaa                                                             >  6:132785/1‑105 (MQ=255)
                                                                                                                             aTTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCaa                                                             <  5:132785/105‑1 (MQ=255)
                                                                                                                                       ctaacccccGAAATTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTg               <  7:156535/136‑1 (MQ=255)
                                                                                                                                       ctaacccccGAAATTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTg               >  8:156535/6‑141 (MQ=255)
                                                                                                                                        cAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtc        <  7:233520/147‑1 (MQ=255)
                                                                                                                                         aGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtc        <  7:87229/146‑1 (MQ=255)
                                                                                                                                         aGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtc        >  8:87229/1‑146 (MQ=255)
                                                                                                                                          acccccGAAATTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGc                                   >  4:95151/3‑118 (MQ=255)
                                                                                                                                          acccccGAAATTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGc                                   <  3:95151/116‑1 (MQ=255)
                                                                                                                                            ccccGAATTTCGTGGGTTAGAGCGATCTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGTCCTATGTTTACATATACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGa  <  2:73862/149‑1 (MQ=255)
                                                                                                                                              ccGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGAc                                              >  7:247660/1‑103 (MQ=255)
                                                                                                                                              ccGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGAc                                              <  8:247660/103‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CTTGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGA  >  NC_002947/2227275‑2227563

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: