Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,282,294:1 +C 12.1% coding (72/300 nt) PP_2883 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,282,2941.C12.1% 108.0 / 12.4 33coding (72/300 nt)PP_2883hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (15/14);  new base C (2/2);  total (17/16)
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.28e-01

AGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACTGTGCGGCAG  >  NC_002947/3282146‑3282442
                                                                                                                                                     |                                                                                                                                                    
aGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtgg                                                                                                                                                       >  4:121760/1‑149 (MQ=255)
        aGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCACCACCTTGAGGTGGCTGCACGcc                                                                                                                                              <  8:226145/149‑1 (MQ=255)
                          tCGGGGGCTTGGCAACGAACGCAATCCGCCAGGCCTGCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAg                                                                                                                            <  6:318850/149‑1 (MQ=255)
                                       aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCCCGCCAAGGGCCACATTGCGCAGGCACTGGTAGACa                                                                                                               >  2:139854/1‑149 (MQ=255)
                                       aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACa                                                                                                               >  6:241192/1‑149 (MQ=255)
                                                aaTGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCg                                                                                                      >  8:163189/1‑149 (MQ=255)
                                                     gcCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGtt                                                                                                    >  6:260580/1‑146 (MQ=255)
                                                     gcCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGtt                                                                                                    <  5:260580/146‑1 (MQ=255)
                                                      ccAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTgcgc                                                                                                                              <  7:34840/119‑1 (MQ=255)
                                                      ccAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTgcgc                                                                                                                              >  8:34840/1‑119 (MQ=255)
                                                          gCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAAc                                                                                            <  3:121760/149‑1 (MQ=255)
                                                            cTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAg                                                                                                                            >  8:5750/1‑115 (MQ=255)
                                                                 gTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCtgtg                                                                                     <  7:163189/149‑1 (MQ=255)
                                                                  tATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTgg                                                                                    <  5:275461/149‑1 (MQ=255)
                                                                      cAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTc                                                                                 <  5:241192/148‑1 (MQ=255)
                                                                        gggCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTc                                                                                 >  2:62150/1‑146 (MQ=255)
                                                                           ccAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGAt                                                                            >  3:295302/1‑148 (MQ=255)
                                                                           ccAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGAt                                                                            >  5:275747/1‑148 (MQ=255)
                                                                                         cAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGc                                                             >  5:208869/1‑149 (MQ=255)
                                                                                              cTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGTAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGtt                                                        >  6:294322/1‑149 (MQ=255)
                                                                                                 gATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTg                                                       <  4:295302/147‑1 (MQ=255)
                                                                                                   tGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGgcgc                                                   <  6:275747/149‑1 (MQ=255)
                                                                                                      cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtga                                                >  1:51040/1‑149 (MQ=255)
                                                                                                      cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtga                                                >  5:220823/1‑149 (MQ=255)
                                                                                                       gACATCGATCAGGCCGCCATCGCGGCGTGTAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtgat                                               <  5:294322/149‑1 (MQ=255)
                                                                                                             gATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGAtctc                                         >  3:39104/1‑149 (MQ=255)
                                                                                                                 aGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGc                                      <  1:62150/148‑1 (MQ=255)
                                                                                                                      gccATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATc                                <  1:139854/149‑1 (MQ=255)
                                                                                                                             cggcgTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTAGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCg                         <  8:74127/149‑1 (MQ=255)
                                                                                                                              ggcgTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGa                        <  2:51040/149‑1 (MQ=255)
                                                                                                                              ggcgTGCACCACCTTGAGGTGGCTGCACGCCAAGGGCCACATTGAGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCg                         >  6:38533/1‑148 (MQ=255)
                                                                                                                               gcgTGCACCACCTTGAGGTGGCTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGa                        <  5:38533/148‑1 (MQ=255)
                                                                                                                                       ccACCTTGAGGTGGCTGCACGCCGAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCtt                 >  4:309399/1‑148 (MQ=255)
                                                                                                                                         aCCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGaaa             >  8:195421/1‑149 (MQ=255)
                                                                                                                                                    gtgCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACTGTGCGGCAg  >  8:102384/1‑149 (MQ=255)
                                                                                                                                                    gtgCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTAAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACTGTGCGGCa   <  4:39104/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACTGTGCGGCAG  >  NC_002947/3282146‑3282442

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: