Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,440,133 G→A 10.2% G352D (GGC→GAC)  PP_3059 → tail sheath protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,440,1330GA10.2% 103.9 / 2.7 39G352D (GGC→GAC) PP_3059tail sheath protein
Reads supporting (aligned to +/- strand):  ref base G (17/18);  new base A (2/2);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.67e-01

GCTCGATCACCGCGACCTACGTCAAGGACGTGACCGAGGGCCTGCAGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACGCTGCCGCTTAAGGAGGCCCCGACATGGCATTCATTCCCCAAATTCT  >  NC_002947/3439992‑3440280
                                                                                                                                             |                                                                                                                                                   
gCTCGATCACCGCGACCTACGTCAAGGACGTGACCGAGGGCCTGCAGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGgtg                                                                                                                                              <  5:56227/149‑1 (MQ=255)
    gATCACCGCGACCTACGTCAAGGACGTGACCGAGGGCCTGCAGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGc                                                                                                                                                    >  1:233202/1‑139 (MQ=255)
    gATCACCGCGACCTACGTCAAGGACGTGACCGAGGGCCTGCAGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGc                                                                                                                                                    <  2:233202/139‑1 (MQ=255)
         ccGCGACCTACGTCAAGGACGTGACCGAGGGCCTGCAGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAAc                                                                                                                                     <  1:71415/149‑1 (MQ=255)
                  aCGTCAAGGACGTGACCGAGGGCCTGCAGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGttt                                                                                                                             >  6:127226/1‑148 (MQ=255)
                             gTGACCGAGGGCCTGCAGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTgcc                                                                                                                 >  5:171838/1‑149 (MQ=255)
                             gTGACCGAGGGCCTGCAGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTgcc                                                                                                                 <  5:127226/149‑1 (MQ=255)
                                aCCGAGGGCCTGCAGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTgccgcc                                                                                                              <  8:42789/149‑1 (MQ=255)
                                  cGAGGGCCTGCAGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCgccgg                                                                                                            >  6:265368/1‑149 (MQ=255)
                                           gCAGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAAcc                                                                                                   >  8:69312/1‑149 (MQ=255)
                                            cAGGCTTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAAcc                                                                                                   <  5:265368/148‑1 (MQ=255)
                                             aGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAAcccc                                                                                                 <  7:69312/149‑1 (MQ=255)
                                             aGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAAccc                                                                                                  >  2:50843/1‑148 (MQ=255)
                                               gCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCaa                                                                                               >  6:3054/1‑149 (MQ=255)
                                                       gcgcGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCcgcg                                                                                       >  2:143748/1‑149 (MQ=255)
                                                         gcgACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTc                                                                                     >  3:158762/1‑149 (MQ=255)
                                                                       cttggCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGACAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAAc                                                                                                    <  2:106742/117‑1 (MQ=255)
                                                                       cAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCaccaac                                                                         >  5:11370/1‑147 (MQ=255)
                                                                       cAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGACAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAAc                                                                                                    >  1:109955/1‑120 (MQ=255)
                                                                       cAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGACAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAAc                                                                                                    <  2:109955/120‑1 (MQ=255)
                                                                       cAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGACAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAAc                                                                                                    >  1:106742/1‑120 (MQ=255)
                                                                            cgcAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGt                                                                  <  4:158762/149‑1 (MQ=255)
                                                                               aaTCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGa                                                               <  6:11370/149‑1 (MQ=255)
                                                                                 tcatcaACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGAcc                                                             >  3:40854/1‑149 (MQ=255)
                                                                                    tcaACTTCGAGGTCTACGCCGACCCGGTGCACAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCg                                                            >  1:24196/1‑147 (MQ=255)
                                                                                     caACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCACTGGTTggtgac                                                           >  6:95265/1‑142 (MQ=255)
                                                                                     caACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCACTGGGTggtgac                                                           <  5:95265/147‑6 (MQ=255)
                                                                                              ggTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACg                                                <  1:143748/149‑1 (MQ=255)
                                                                                                 cTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACGCTg                                             >  5:49638/1‑149 (MQ=255)
                                                                                                  tACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGt                                                                                                                               <  5:62106/66‑1 (MQ=255)
                                                                                                  tACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGt                                                                                                                               >  6:62106/1‑66 (MQ=255)
                                                                                                    cGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACGCTg                                             <  5:3054/146‑1 (MQ=255)
                                                                                                    cGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACGCTGcc                                           <  6:171838/148‑1 (MQ=255)
                                                                                                          cgcGGTGATCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCGAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACGCTGCCGCTTAag                                    <  2:24196/147‑1 (MQ=255)
                                                                                                                  tCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACGCTGCCGCTTAAGGAGGcccc                            >  4:233691/1‑149 (MQ=255)
                                                                                                                    aaCACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACGCTGCCGCTTAAGGAGGCCCCg                           <  4:21760/148‑1 (MQ=255)
                                                                                                                         ggCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACGCTGCCGCTTAAGGAGGCCCCGACATgg                     >  5:267318/1‑149 (MQ=255)
                                                                                                                             agccagcTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACGCTGCCGCTTAAGGAGGCCCCGACATGGcatt                 <  4:62707/149‑1 (MQ=255)
                                                                                                                                  gcTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACGCTGCCGCTTAAGGAGGCCCCGACATGGCattcattc            <  6:49638/149‑1 (MQ=255)
                                                                                                                                            ggCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGaa                                                          <  7:34247/93‑1 (MQ=255)
                                                                                                                                            ggCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGaa                                                          >  8:34247/1‑93 (MQ=255)
                                                                                                                                            ggCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACGCTGCCGCTTAAGGAGGCCCCGACTTGGCATTCATTCCCCAAATtct  >  8:109093/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GCTCGATCACCGCGACCTACGTCAAGGACGTGACCGAGGGCCTGCAGGCGTTCATGCGCGACCTGAAGAACCAAGGCGCAATCATCAACTTCGAGGTCTACGCCGACCCGGTGCTCAACACGGCCAGCCAGCTGGAGCAGGGCAAGGTGTACTGGAACATTCGTTTCACCGACGTGCCGCCGGCCGAAAACCCCAACTTCCGCGTCGAGGTCACCAACCAGTGGTTGACCGAAGTGCTTAACGCTGCCGCTTAAGGAGGCCCCGACATGGCATTCATTCCCCAAATTCT  >  NC_002947/3439992‑3440280

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