Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,592,806 A→T 10.2% Y334F (TAC→TTC)  PP_2270 → DNA primase/helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,592,8060AT10.2% 128.5 / 3.7 47Y334F (TAC→TTC) PP_2270DNA primase/helicase
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base T (2/3);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGGAAACGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTGGGAATCATCATCACGT  >  NC_002947/2592659‑2592949
                                                                                                                                                   |                                                                                                                                                
ccGGGAAACGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTa                                                                                                                                                  >  2:96252/1‑148 (MQ=255)
    gAAACGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCACCTTGCGCTTCacaa                                                                                                                                             <  1:155097/149‑4 (MQ=255)
     aaaCGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGc                                                                                                                                            <  8:123700/149‑1 (MQ=255)
       aCGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTg                                                                                                                                          >  6:122673/1‑149 (MQ=255)
        cGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCACCTTGCGCTTCa                                                                                                                                                >  8:211721/1‑142 (MQ=255)
        cGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCACCTTGCGCTTCa                                                                                                                                                <  7:211721/142‑1 (MQ=255)
        cGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCACCTTGCGCTTCa                                                                                                                                                >  4:245377/1‑142 (MQ=255)
        cGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCACCTTGCGCTTCa                                                                                                                                                <  3:245377/142‑1 (MQ=255)
         gCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGtc                                                                                                                                        >  8:320596/1‑149 (MQ=255)
               aCGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATgg                                                                                                                                             >  5:986/1‑138 (MQ=255)
                    cccAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCAt                                                                                                                            >  8:170805/1‑149 (MQ=255)
                     ccAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCAt                                                                                                                            >  8:327952/1‑148 (MQ=255)
                      cAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCa                                                                                                                          <  8:290646/149‑1 (MQ=255)
                       aaGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATc                                                                                                                           >  4:263898/1‑147 (MQ=255)
                         ggACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGc                                                                                                                       <  3:263898/149‑1 (MQ=255)
                          gACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATgg                                                                                                                                             >  8:311799/1‑127 (MQ=255)
                          gACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATgg                                                                                                                                             <  7:311799/127‑1 (MQ=255)
                           acacCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCtt                                                                                                                     <  5:122673/149‑1 (MQ=255)
                             acCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTAt                                                                                                                   <  7:274294/149‑1 (MQ=255)
                                        gACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGac                                                                                                         <  7:170805/148‑1 (MQ=255)
                                               gaggagCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACa                                                                                                              >  1:13094/1‑136 (MQ=255)
                                               gaggagCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACa                                                                                                              <  2:13094/136‑1 (MQ=255)
                                                            cAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGac                                                                                    <  7:327952/149‑1 (MQ=255)
                                                                      tGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCTCCCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGa                                                                            <  1:96252/148‑1 (MQ=255)
                                                                         gcTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATgag                                                                        >  7:245792/1‑148 (MQ=255)
                                                                               cGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCt                                                                 >  8:171656/1‑149 (MQ=255)
                                                                                gAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTc                                                                >  5:198425/1‑149 (MQ=255)
                                                                                gAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTc                                                                >  1:155289/1‑149 (MQ=255)
                                                                                         gTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGc                                                        <  2:155289/148‑1 (MQ=255)
                                                                                            aTGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTa                                                                                                                    >  8:19530/1‑85 (MQ=255)
                                                                                            aTGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTa                                                                                                                    <  7:19530/85‑1 (MQ=255)
                                                                                                   acgCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCTATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGaaa                                             <  7:320596/149‑1 (MQ=255)
                                                                                                    cgCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGa                                               >  2:26013/1‑146 (MQ=255)
                                                                                                      cATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAagag                                          <  7:171656/149‑1 (MQ=255)
                                                                                                       aTTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGCTGAGAAGGGCTCGCTTGAGCAGATCATGaaa                                             >  2:87908/1‑145 (MQ=255)
                                                                                                       aTTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGaaa                                             >  1:62866/1‑145 (MQ=255)
                                                                                                                   cGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTgg                             <  6:198425/149‑1 (MQ=255)
                                                                                                                              ttGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATga                                                                         <  7:257950/94‑1 (MQ=255)
                                                                                                                              ttGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATga                                                                         >  8:257950/1‑94 (MQ=255)
                                                                                                                               tgtgAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTgg                 >  7:101730/1‑149 (MQ=255)
                                                                                                                               tgtgAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTgg                 <  8:245792/149‑1 (MQ=255)
                                                                                                                                   aaGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGCTGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTGGGaa              <  1:87908/148‑1 (MQ=255)
                                                                                                                                   aaGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAg                                                         >  4:266216/1‑105 (MQ=255)
                                                                                                                                   aaGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAg                                                         <  3:266216/105‑1 (MQ=255)
                                                                                                                                   aaGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTGGGAat             >  2:115262/1‑149 (MQ=255)
                                                                                                                                      cgcAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGCACGAGCTGGGAatcat          >  4:253819/1‑149 (MQ=255)
                                                                                                                                       gcAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTGGGAatcatc         <  1:115262/149‑1 (MQ=255)
                                                                                                                                              gcgcTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTGGGAATCATCATCACGt  >  4:45439/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCGGGAAACGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTGGGAATCATCATCACGT  >  NC_002947/2592659‑2592949

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