Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,597,751 C→A 10.5% Q686H (CAG→CAT) ‡ PP_4071 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,597,7510CA10.5% 108.9 / 4.5 38Q686H (CAG→CAT) ‡PP_4071membrane protein
Reads supporting (aligned to +/- strand):  ref base C (19/15);  new base A (2/2);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

CTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCA  >  NC_002947/4597612‑4597896
                                                                                                                                           |                                                                                                                                                 
cTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAgg                                                                                                                                          <  1:202651/149‑1 (MQ=255)
    ggCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGccc                                                                                                                                       <  6:76824/148‑1 (MQ=255)
      cAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCa                                                                                                                                    <  5:188521/149‑1 (MQ=255)
       aGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGAAGCAGCAGGcccagga                                                                                                                                   <  3:189640/149‑5 (MQ=255)
         tGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGCCCTGCTGCTTCAGGCCCGCATac                                                                                                                                 >  6:54126/1‑149 (MQ=255)
                            aaCACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtgtg                                                                                                                       >  1:38759/1‑140 (MQ=255)
                            aaCACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtgtg                                                                                                                       <  2:38759/140‑1 (MQ=255)
                                      tCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCa                                                                                                     <  5:171910/148‑1 (MQ=255)
                                              gCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGCCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTctgct                                                                                            <  5:54126/149‑1 (MQ=255)
                                                  gAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCAc                                                                                        >  5:21420/1‑149 (MQ=255)
                                                     ccGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGAAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTg                                                                                     >  1:75650/1‑149 (MQ=255)
                                                      cGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTctgctg                                                                                           <  1:133128/142‑1 (MQ=255)
                                                      cGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTctgctg                                                                                           >  2:133128/1‑142 (MQ=255)
                                                              aGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACg                                                                            >  8:187070/1‑149 (MQ=255)
                                                              aGCTCGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGACCAGACCGGTATGCAACTGATCCTGGGCCTGCTGCTTAAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGACTGCTGCACGTGGCTTACACg                                                                            >  2:259870/1‑149 (MQ=255)
                                                                    tAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGc                                                                       >  4:153475/1‑148 (MQ=255)
                                                                     aaaCGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGa                                                                     >  4:71722/1‑149 (MQ=255)
                                                                                     aGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACCGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCgg                                                     >  2:263750/1‑149 (MQ=255)
                                                                                                      gCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAAc                                                                                                                  >  7:145774/1‑71 (MQ=255)
                                                                                                      gCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAAc                                                                                                                  <  8:145774/71‑1 (MQ=255)
                                                                                                      gCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATAc                                                                  >  1:135766/1‑119 (MQ=255)
                                                                                                      gCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATAc                                                                  <  2:135766/119‑1 (MQ=255)
                                                                                                      gCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAgc                                    <  3:71722/149‑1 (MQ=255)
                                                                                                      gCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAgc                                    >  3:15121/1‑149 (MQ=255)
                                                                                                          gTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGt                                                    >  8:33057/1‑129 (MQ=255)
                                                                                                          gTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGt                                                    <  7:33057/129‑1 (MQ=255)
                                                                                                                   ggTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGt                                        <  8:63712/132‑1 (MQ=255)
                                                                                                                   ggTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGt                                        >  7:63712/1‑132 (MQ=255)
                                                                                                                        gCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCg                   >  8:106909/1‑148 (MQ=255)
                                                                                                                         cAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCa                 >  1:182041/1‑149 (MQ=255)
                                                                                                                          aaCTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCAt                >  6:155242/1‑149 (MQ=255)
                                                                                                                            cTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGa              <  2:253485/149‑1 (MQ=255)
                                                                                                                                    gggCCTGAAGCAGCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATAc                              >  7:279828/1‑125 (MQ=255)
                                                                                                                                    gggCCTGAAGCAGCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATAc                              <  8:279828/125‑1 (MQ=255)
                                                                                                                                    gggCCTGAAGCAGCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGt      >  5:275670/1‑149 (MQ=255)
                                                                                                                                     ggCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGtt     <  3:153475/149‑1 (MQ=255)
                                                                                                                                       ccTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTg    >  6:121187/1‑148 (MQ=255)
                                                                                                                                       ccTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGc   <  4:15121/149‑1 (MQ=255)
                                                                                                                                          gctgctTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCa  >  4:163542/1‑147 (MQ=255)
                                                                                                                                          gctgctTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCa  <  3:163542/147‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
CTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCA  >  NC_002947/4597612‑4597896

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: