Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,497,329 Δ1 bp 9.1% coding (268/765 nt) PP_5562 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,497,3290C. 9.1% 134.9 / 7.2 44coding (268/765 nt)PP_5562hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/22);  new base . (2/2);  total (20/24)
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.64e-01

AAGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGAT  >  NC_002947/3497193‑3497467
                                                                                                                                        |                                                                                                                                           
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                                                                               >  2:91027/1‑139 (MQ=255)
                                                           gtgtTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTgg                                                                               <  1:91027/139‑1 (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                                                          >  2:162709/1‑129 (MQ=255)
                                                                                          aaCAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCgcgc                                                          <  1:162709/129‑1 (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                                                                          >  3:143716/1‑87 (MQ=255)
                                                                                                                    gTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTg                                                                          <  4:143716/87‑1 (MQ=255)
                                                                                                                         gCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcgga       <  7:22106/149‑1 (MQ=255)
                                                                                                                          cGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggac      <  4:44647/149‑1 (MQ=255)
                                                                                                                            gTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTc                          >  1:208639/1‑127 (MQ=255)
                                                                                                                            gTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTc                          <  2:208639/127‑1 (MQ=255)
                                                                                                                            gTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggatgg    <  4:252768/149‑1 (MQ=255)
                                                                                                                             tACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggacgga   >  6:291051/1‑149 (MQ=255)
                                                                                                                             tACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGc          >  6:281826/1‑142 (MQ=255)
                                                                                                                             tACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGc          <  5:281826/142‑1 (MQ=255)
                                                                                                                               gttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggacgg    <  4:139333/145‑1 (MQ=255)
                                                                                                                               gttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggacgg    >  3:139333/2‑146 (MQ=255)
                                                                                                                                ttCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGAt  >  6:286445/1‑147 (MQ=255)
                                                                                                                                ttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCt                                                                                 <  4:196757/68‑1 (MQ=37)
                                                                                                                                ttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCt                                                                                 >  3:196757/1‑68 (MQ=37)
                                                                                                                                        |                                                                                                                                           
AAGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGAT  >  NC_002947/3497193‑3497467

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