Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,497,333 G→C 8.7% V88V (GTC→GTG)  PP_5562 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,497,3330GC 8.7% 144.1 / 4.9 46V88V (GTC→GTG) PP_5562hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/23);  new base C (2/2);  total (21/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 = 4.44e-01

AAGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGCGGTGCCACGG  >  NC_002947/3497193‑3497481
                                                                                                                                            |                                                                                                                                                     
aaGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCTGTAACAGTTCTGGGGTGGTTCTGATCATGTTCAGCGGTAGTTCCAGCTCGGTGACC‑AGCTg                                                                                                                                              <  4:290269/149‑1 (MQ=255)
    cGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAAc                                                                                                                                          >  6:79756/1‑149 (MQ=255)
     gTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACg                                                                                                                                         <  5:79756/149‑1 (MQ=255)
        gAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCCGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGgcg                                                                                                                                      >  7:331532/1‑149 (MQ=255)
             ccaccaCCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgcc                                                                                                                                 <  4:271248/149‑1 (MQ=255)
                              gggCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTg                                                                                                                                              <  6:105409/119‑1 (MQ=255)
                              gggCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTg                                                                                                                                              >  5:105409/1‑119 (MQ=255)
                                                 ccTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTg                                                                                              >  1:331440/1‑148 (MQ=255)
                                                  cTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgccgagccga                                                                                                                          <  6:234539/119‑1 (MQ=255)
                                                  cTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgccgagccga                                                                                                                          >  5:234539/1‑119 (MQ=255)
                                                       cGCTGTGTTCCAGGGTCTTGTCGTATTCCTGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTTGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTgg                                                                                       <  1:250802/149‑1 (MQ=255)
                                                        gCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGc                                                                                      <  5:251415/149‑1 (MQ=255)
                                                           gtgtTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTgg                                                                                             <  1:91027/139‑1 (MQ=255)
                                                           gtgtTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTgg                                                                                             >  2:91027/1‑139 (MQ=255)
                                                            tgtTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTg                                                                                  <  5:137036/149‑1 (MQ=255)
                                                              ttCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTgg                                                                                 <  2:331440/148‑1 (MQ=255)
                                                                    ggTCTTGTCGTATTCCACCGGGCACCGTTCGGGGGTGGGGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTcgc                                                                          <  7:218207/149‑1 (MQ=255)
                                                                        ttGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGccc                                                                      >  8:68372/1‑149 (MQ=255)
                                                                                   cAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCt                                                           >  4:150827/1‑149 (MQ=255)
                                                                                       gggAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTgc                                                       <  3:150827/149‑1 (MQ=255)
                                                                                        ggAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTgcg                                                      >  3:252768/1‑149 (MQ=255)
                                                                                          aaCAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCgcgc                                                                        <  1:162709/129‑1 (MQ=255)
                                                                                          aaCAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCgcgc                                                                        >  2:162709/1‑129 (MQ=255)
                                                                                             aGTTCGGGGATGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTg                                                 >  3:253311/1‑149 (MQ=255)
                                                                                               ttCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCa                                               <  2:221902/149‑1 (MQ=255)
                                                                                                tCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAg                                              <  7:68372/149‑1 (MQ=255)
                                                                                                   ggggTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTa                                           >  4:213504/1‑149 (MQ=255)
                                                                                                        ggtgCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGc                                      >  8:22106/1‑149 (MQ=255)
                                                                                                        ggtgCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGc                                      <  3:213504/149‑1 (MQ=255)
                                                                                                                    gTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTg                                                                                        <  4:143716/87‑1 (MQ=255)
                                                                                                                    gTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTg                                                                                        >  3:143716/1‑87 (MQ=255)
                                                                                                                         gCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcgga                     <  7:22106/149‑1 (MQ=255)
                                                                                                                          cGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggac                    <  4:44647/149‑1 (MQ=255)
                                                                                                                            gTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTc                                        <  2:208639/127‑1 (MQ=255)
                                                                                                                            gTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTc                                        >  1:208639/1‑127 (MQ=255)
                                                                                                                            gTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggatgg                  <  4:252768/149‑1 (MQ=255)
                                                                                                                             tACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggacgga                 >  6:291051/1‑149 (MQ=255)
                                                                                                                             tACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGc                        <  5:281826/142‑1 (MQ=255)
                                                                                                                             tACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGc                        >  6:281826/1‑142 (MQ=255)
                                                                                                                               gttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggacgg                  >  3:139333/2‑146 (MQ=255)
                                                                                                                               gttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggacgg                  <  4:139333/145‑1 (MQ=255)
                                                                                                                                ttCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGAt                >  6:286445/1‑147 (MQ=255)
                                                                                                                                ttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCt                                                                                               >  3:196757/1‑68 (MQ=37)
                                                                                                                                ttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCt                                                                                               <  4:196757/68‑1 (MQ=37)
                                                                                                                                         ggTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGCgg          >  6:217701/1‑144 (MQ=255)
                                                                                                                                         ggTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGCgg          <  5:217701/144‑1 (MQ=255)
                                                                                                                                            gACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGCGGTGCCACgg  >  5:299109/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
AAGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGCGGTGCCACGG  >  NC_002947/3497193‑3497481

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