Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,406,173 G→A 7.4% G43G (GGC→GGT)  PP_3018 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,406,1730GA 7.4% 153.6 / 2.6 54G43G (GGC→GGT) PP_3018hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/23);  new base A (2/2);  total (29/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 = 7.44e-01

CAAAGCCTGGCGCTGGACCTGTAAGGCGTCAACCGAACCCGCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGCTAGAAACAGATAGATCATGCAAGGTACGGGTTCAAGGG  >  NC_002947/3406031‑3406321
                                                                                                                                              |                                                                                                                                                    
cAAAGCCTGGCGCTGGACCTGTAAGGCGTCAACCGAACCCGCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGc                                                                                                                                                <  3:82053/149‑1 (MQ=255)
 aaaGCCTGGCGCTGGACCTGTAAGGCGTCAACCGAACCCGCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCg                                                                                                                                               >  5:307632/1‑149 (MQ=255)
               caCCTGTAAGCCGCCAACCGACCCCGCGGCCCTTCAGGCTGGCCAGAACGGCCAGCACCGAACCGGCGAGCAGGGGGAATTGCCGCTGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCgg                                                                                                                                   <  1:241694/146‑1 (MQ=255)
                aCCTGTAAGGCGTCAACCGAACCCGCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGcc                                                                                                                                >  8:190683/1‑149 (MQ=255)
                    gTAAGGCGTCAACCGAACCCGCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGgc                                                                                                                                          <  4:165086/135‑1 (MQ=255)
                    gTAAGGCGTCAACCGAACCCGCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGgc                                                                                                                                          >  3:165086/1‑135 (MQ=255)
                     tAAGGCGTCAACCGAACCCGCTTGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATg                                                                                                                            <  8:44201/148‑1 (MQ=255)
                              aaCCGAACCCGCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAAccc                                                                                                                  >  5:9085/1‑149 (MQ=255)
                                ccGAACCCGCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTc                                                                                                                >  5:202833/1‑149 (MQ=255)
                                  gtaCCCGCTGGCGATGAGGCAGGCCAGCACGGCCTGCATCGAACCGGCGATCAGGGTCTATTGCCGCCGGGCCTGACGTGCGCCGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGGAAACCCTCgg                                                                                                              <  2:162740/147‑1 (MQ=255)
                                  gAACCCGCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCgg                                                                                                              <  7:190683/149‑1 (MQ=255)
                                   aaCCCGCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGc                                                                                                             <  7:135510/149‑1 (MQ=255)
                                     cccGCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCca                                                                                                           >  5:2004/1‑149 (MQ=255)
                                        gCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCCAACCCTCGGCCACAt                                                                                                        >  1:184489/1‑149 (MQ=255)
                                            gCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCgg                                                                                                     >  4:154023/1‑148 (MQ=255)
                                            gCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGc                                                                                                    <  6:2004/149‑1 (MQ=255)
                                              gATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGcc                                                                                                                                >  1:280588/1‑119 (MQ=255)
                                              gATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGcc                                                                                                                                <  2:280588/119‑1 (MQ=255)
                                              gATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATctgtag                                                                                                  >  1:263451/1‑144 (MQ=255)
                                              gATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGc                                                                                                    >  5:73340/1‑147 (MQ=255)
                                                       ggccagcacggccagcacCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCt                                                                                         <  3:154023/149‑1 (MQ=255)
                                                        gccagcacggccagcacCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCACCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCt                                                                                         >  3:212418/1‑148 (MQ=255)
                                                        gccagcacggccagcacCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCACCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCt                                                                                         >  3:216453/1‑148 (MQ=255)
                                                          cagcacggccagcacCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGaa                                                                                      >  8:204808/1‑149 (MQ=255)
                                                               cggccagcacCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGc                                                                                 <  1:12755/149‑1 (MQ=255)
                                                                         cGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGa                                                                                       <  7:235415/133‑1 (MQ=255)
                                                                         cGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGa                                                                                       >  8:235415/1‑133 (MQ=255)
                                                                            aCCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGc                                                                    <  6:73340/149‑1 (MQ=255)
                                                                                 cGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAg                                                               <  7:204808/149‑1 (MQ=255)
                                                                                 cGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAg                                                               <  6:9085/149‑1 (MQ=255)
                                                                                       gggTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTaccccgc                                                          >  7:18081/1‑142 (MQ=255)
                                                                                       gggTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCACCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAg                                                         <  4:216453/149‑1 (MQ=255)
                                                                                       gggTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCACCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAg                                                         <  4:212418/149‑1 (MQ=255)
                                                                                              ttGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTaa                                                   <  7:232803/148‑1 (MQ=255)
                                                                                                gccgccGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACt                                                 >  2:34173/1‑148 (MQ=255)
                                                                                                 ccgccgGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGTTGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAg                                                         <  3:12469/139‑1 (MQ=255)
                                                                                                 ccgccgGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGCTGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAg                                                         >  4:12469/1‑139 (MQ=255)
                                                                                                 ccgccgGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCa                                               <  6:307632/149‑1 (MQ=255)
                                                                                                  cgccgGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAg                                              >  2:298275/1‑149 (MQ=255)
                                                                                                           tGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGc                                        >  3:109334/1‑146 (MQ=255)
                                                                                                             aCGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGCTAGa                                    >  7:176578/1‑148 (MQ=255)
                                                                                                               gTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGAcc                                                                        >  8:203265/1‑110 (MQ=255)
                                                                                                               gTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGAcc                                                                        <  7:203265/110‑1 (MQ=255)
                                                                                                                  cgcgGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGCTAGAAACaga                              <  6:202833/149‑1 (MQ=255)
                                                                                                                  cgcgGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGCTAGAAACaga                              <  1:298275/149‑1 (MQ=255)
                                                                                                                          ccGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCTAACCCTCGTCGACATGGGCAAGTTGAACCTGAATCAGCATACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGCTAGAAACagatagttgat                      <  2:260593/149‑6 (MQ=255)
                                                                                                                                  aCGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGCTAGAAACAGATAGATCATGCAAGGt               >  2:211292/1‑148 (MQ=255)
                                                                                                                                  aCGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGCTAGAAACAGATAGATCATGCAAGGTa              >  5:65424/1‑149 (MQ=255)
                                                                                                                                      aGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGCTAGAAACAGATAGATCATGCAAGGTACggg          >  8:50256/1‑149 (MQ=255)
                                                                                                                                        gTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACg                                         >  7:42372/1‑116 (MQ=255)
                                                                                                                                        gTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACg                                         <  8:42372/116‑1 (MQ=255)
                                                                                                                                          aCTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGACAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGc                                        >  2:123252/1‑115 (MQ=255)
                                                                                                                                          aCTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGACAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGCTAGAAACAGAGAGATCATGCAAGGTACGGGTTc       >  2:127666/1‑148 (MQ=255)
                                                                                                                                            tCGCCGGGCGTAGGCGCATCCGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGCTAGAAACAGATAGATCATGCAAGGTACGGCTTCAAg    >  4:100018/1‑149 (MQ=255)
                                                                                                                                              gCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGCCCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGCTAGAAACAGATAGATCATGCAAGGTACGGGTTCAAggg  >  7:246476/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CAAAGCCTGGCGCTGGACCTGTAAGGCGTCAACCGAACCCGCTGGCGATCAGGCAGGCCAGCACGGCCAGCACCGAACCGGCGATCAGGGTGTATTGCCGCCGGGCCTGACGTGCGCGGGGCCCGACCTCACGCAGGTACTCGCCGGGCGTAGGCGCATCGGGCCGATGGCGCAAACCCTCGGCCACATCGGCAAGTTGACCCTGAATCAGCAAACCGACCAGGCAGTAGGTGCAGGCGTAACTCAGCAACGCTAGAAACAGATAGATCATGCAAGGTACGGGTTCAAGGG  >  NC_002947/3406031‑3406321

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