Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,586,238 G→A 7.8% A9V (GCG→GTG)  PP_4913 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,586,2380GA 7.8% 143.9 / 2.3 52A9V (GCG→GTG) PP_4913hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/27);  new base A (2/2);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

CTCAGCAGGCGCACTTCGATAAAGCCTTGACCATGCTCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCAC  >  NC_002947/5586091‑5586385
                                                                                                                                                    |                                                                                                                                                   
cTCAGCAGGCGCACTTCGATAAAGCCTTGACCATGCTCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑GCCCcgt                                                                                                                                                    >  6:319115/1‑148 (MQ=255)
                cgatccaGCCTTGACCATGCTCTGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGgcgc                                                                                                                                     <  4:207064/142‑1 (MQ=255)
                         cTTGACCATGCTCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATg                                                                                                                           <  8:99036/149‑1 (MQ=255)
                         cTTGACCATGCTCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATg                                                                                                                           <  4:42121/149‑1 (MQ=255)
                            gACCATGCTCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCTTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGc                                                                                                                                >  5:325921/1‑141 (MQ=255)
                            gACCATGCTCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGc                                                                                                                                <  6:325921/141‑1 (MQ=255)
                             aCCATGCTCGGATTGCACCTGGTGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGttt                                                                                                                       <  1:227937/149‑1 (MQ=255)
                             aCCATGCTCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGttt                                                                                                                       <  1:225191/149‑1 (MQ=255)
                              ccATGCTCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTg                                                                                                                      <  8:122212/149‑1 (MQ=255)
                               cATGCTCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTgg                                                                                                                     <  4:19425/149‑1 (MQ=255)
                               cATGCTCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTgg                                                                                                                     >  3:128007/1‑149 (MQ=255)
                                    tCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTc                                                                                                                >  6:173564/1‑149 (MQ=255)
                                    tCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTc                                                                                                                >  4:143102/1‑149 (MQ=255)
                                      ggATTGCACCTGGCGGTTGGGGTGACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCa                                                                                                              <  4:49575/149‑1 (MQ=255)
                                             aCCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCcgcgc                                                                                                                                                  <  7:158611/106‑1 (MQ=255)
                                             aCCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCcgcgc                                                                                                                                                  >  8:158611/1‑106 (MQ=255)
                                             aCCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCACGGAGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCt                                                                                                       >  1:155930/1‑149 (MQ=255)
                                               cTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTc                                                                                                     >  3:160016/1‑149 (MQ=255)
                                               cTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTc                                                                                                     <  8:125988/149‑1 (MQ=255)
                                               cTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTc                                                                                                     <  1:27902/149‑1 (MQ=255)
                                               cTGGCGGTTGGGGTTACGCTGGCTGTCTACCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGACGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTc                                                                                                     >  8:175458/1‑149 (MQ=255)
                                                        ggggTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGt                                                                                             <  2:231157/148‑1 (MQ=255)
                                                           gTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGgc                                                                                                                                       >  1:189286/1‑103 (MQ=255)
                                                           gTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGgc                                                                                                                                       <  2:189286/103‑1 (MQ=255)
                                                              aCGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAgcg                                                                                      >  6:312717/1‑149 (MQ=255)
                                                                   ggctgTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAga                                                                                 >  6:175487/1‑149 (MQ=255)
                                                                      tgTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAgacga                                                                              >  3:54005/1‑149 (MQ=255)
                                                                            cccACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCGgccgc                                                                                                         <  6:214718/116‑1 (MQ=255)
                                                                            cccACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCGgccgc                                                                                                         >  5:214718/1‑116 (MQ=255)
                                                                             ccACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCACGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTcc                                                                                                    <  6:13711/120‑1 (MQ=255)
                                                                             ccACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCACGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTcc                                                                                                    <  6:18449/120‑1 (MQ=255)
                                                                             ccACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCACGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTcc                                                                                                    >  5:18449/1‑120 (MQ=255)
                                                                             ccACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCACGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTcc                                                                                                    >  5:13711/1‑120 (MQ=255)
                                                                              ctcTGGCCGAAAAAGCTCGTGGCAAGAAACAGCTCCGCCGCCCAGGCGACGATCATCAGTATCT‑TGCCCGCGCGCTTGCCCGGCGCCTGTCTCTTGGTTTGGCGGTCGAGCCGCCTTCGGGTGGGTCGAAGCTCAAGACGAccggccgg                                                                      <  2:155930/147‑1 (MQ=255)
                                                                               aCTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCggg                                                                     >  2:270875/1‑149 (MQ=255)
                                                                                cTGGCCGAAAAAGCGCGGGGCAACATCCAGCGCCGCCGCCCAGGCGACGATCATCCGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCATCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGt                                                                    <  4:281917/149‑1 (MQ=255)
                                                                                   gCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACTAAGCCCCCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGttt                                                                  >  1:259778/1‑149 (MQ=255)
                                                                                   gCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGttt                                                                  >  5:215853/1‑148 (MQ=255)
                                                                                      gAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAg                                                                                        <  8:292683/123‑1 (MQ=255)
                                                                                      gAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAg                                                                                        >  7:292683/1‑123 (MQ=255)
                                                                                                gTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACgg                                                    <  5:173564/149‑1 (MQ=255)
                                                                                                  ggCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGgc                                                  >  1:220883/1‑149 (MQ=255)
                                                                                                    cAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGgcgc                                                >  5:89179/1‑149 (MQ=255)
                                                                                                         aaCAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCgtt                                            <  4:160016/148‑1 (MQ=255)
                                                                                                           cAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCgttgtt                                         <  6:22341/149‑1 (MQ=255)
                                                                                                            aGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCgttgttg                                        <  5:175487/149‑1 (MQ=255)
                                                                                                                ccgccgccCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTtgttgttgt                                    <  2:29887/149‑1 (MQ=255)
                                                                                                                             gacgaTCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATc                         >  7:266573/1‑147 (MQ=255)
                                                                                                                              acgaTCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAg                      <  6:89179/149‑1 (MQ=255)
                                                                                                                                  tcatcaGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGgcgc                  <  5:312717/149‑1 (MQ=255)
                                                                                                                                     tcaGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCt               >  5:285947/1‑149 (MQ=255)
                                                                                                                                       aGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCtt              <  2:259778/148‑1 (MQ=255)
                                                                                                                                           aCT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGt                                                                                                 >  2:122257/1‑61 (MQ=255)
                                                                                                                                           aCT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGt                                                                                                 <  1:122257/61‑1 (MQ=255)
                                                                                                                                                gCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGa      <  6:215853/148‑1 (MQ=255)
                                                                                                                                                    gcgcgcTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCAc  >  3:245685/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTCAGCAGGCGCACTTCGATAAAGCCTTGACCATGCTCGGATTGCACCTGGCGGTTGGGGTTACGCTGGCTGTCTTCCCACTGGCCGAAAAAGCGCGTGGCAAGAAACAGCGCCGCCGCCCAGGCGACGATCATCAGTAACT‑TGCCCGCGCGCTTGCCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCAC  >  NC_002947/5586091‑5586385

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