Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,051,969 A→G 7.7% A466A (GCT→GCC) ‡ PP_0910 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,051,9690AG 7.7% 146.3 / 2.5 52A466A (GCT→GCC) ‡PP_0910hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/25);  new base G (2/2);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GATGGCGGCGGAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTC  >  NC_002947/1051824‑1052103
                                                                                                                                                 |                                                                                                                                       
gATGGCGGCGGAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCGCGGCTCCTTTGCGCGTTTAAAGTGATTTCAGGAACGCGAGCAACGCAgcg                                                                                                                                      >  8:161216/1‑149 (MQ=255)
gATGGCGGCGGAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAgcg                                                                                                                                      >  4:306182/1‑149 (MQ=255)
      ggcggAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGACCATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCgctgc                                                                                                                                >  5:202349/1‑149 (MQ=255)
            ggTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCa                                                                                                                          <  5:322667/149‑1 (MQ=255)
             gTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTCGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCa                                                                                                                          >  8:61010/1‑148 (MQ=255)
                    cTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGaaa                                                                                                                    <  1:285813/147‑1 (MQ=255)
                    cTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGt                                                                                                                  <  8:36985/149‑1 (MQ=255)
                    cTTGCGGGTCTTGGGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCGCCTTTGCGCGTTTAAAGTGAGTTCAGGAACCCGAGCAACGCAGCGCGCTGCTCGGCATCGaaa                                                                                                                    <  2:152620/147‑1 (MQ=255)
                       gCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTg                                                                                                                 <  6:320475/147‑1 (MQ=255)
                          ggTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTAc                                                                                                            >  6:32538/1‑149 (MQ=255)
                           gTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTAc                                                                                                            >  8:57675/1‑148 (MQ=255)
                                           gCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCt                                                                                           >  6:52992/1‑149 (MQ=255)
                                           gCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCt                                                                                           >  4:112987/1‑149 (MQ=255)
                                           gCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCt                                                                                           >  3:157974/1‑149 (MQ=255)
                                           gCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAAAAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGAGAGCAACGCAGCGAGCTGCTCGGCATCGAAAGTGAGTACGAAATAGCGtcacgaca                                                                                           >  5:203648/1‑142 (MQ=255)
                                            cACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTc                                                                                          <  4:157974/149‑1 (MQ=255)
                                            cACCGAGGGCGAGTGCGGCGAGGCTGGGGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTc                                                                                          <  6:202349/149‑1 (MQ=255)
                                               cGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGc                                                                                       <  7:61010/149‑1 (MQ=255)
                                                      gAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAATTGAGTTCAGGAACGCGAGCAACGCAgcgc                                                                                                                                     <  2:205501/96‑1 (MQ=255)
                                                      gAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAATTGAGTTCAGGAACGCGAGCAACGCAgcgc                                                                                                                                     >  1:205501/1‑96 (MQ=255)
                                                      gAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCa                                                                                <  7:265532/149‑1 (MQ=255)
                                                            ggcgAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGc                                                                          <  8:20205/149‑1 (MQ=255)
                                                              cgAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCca                                                                        <  5:52992/149‑1 (MQ=255)
                                                                  gCTGGTGAACAACAGTTTGGGTCGGAACATCCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAg                                                                     >  1:26430/1‑148 (MQ=255)
                                                                      gTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑gg                                                               >  3:108914/1‑149 (MQ=255)
                                                                         aacaacAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGcc                                                             <  3:133596/148‑1 (MQ=255)
                                                                            aacaGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTc                                                           >  6:244571/1‑147 (MQ=255)
                                                                                 tttGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAgg                                                    >  1:30738/1‑149 (MQ=255)
                                                                                       tCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCgg                                              <  5:244571/149‑1 (MQ=255)
                                                                                              aTGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGAGGGCGCgga                                       >  2:239652/1‑148 (MQ=255)
                                                                                                  ccGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTc                                   <  3:112987/149‑1 (MQ=255)
                                                                                                      cTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGc                               <  5:32538/149‑1 (MQ=255)
                                                                                                       tCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATg                                >  3:8644/1‑147 (MQ=255)
                                                                                                         cTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAg                             >  2:250500/1‑148 (MQ=255)
                                                                                                               gcgTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCTCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GTCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAAGTgg                      <  6:203648/149‑1 (MQ=255)
                                                                                                                   ttAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACTTGGGGCAAGAAGTTTGAGGGCGCGCGCGTCAAGAAGGAAATGGCTa                   >  1:263916/1‑149 (MQ=255)
                                                                                                                     aaaGTGAGTTCAGGAACGCGAGCAACGCAGCGGTCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATg                 <  3:140922/148‑1 (MQ=255)
                                                                                                                     aaaGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATg                 <  2:30738/148‑1 (MQ=255)
                                                                                                                         tGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCAc            <  3:306182/149‑1 (MQ=255)
                                                                                                                            gTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGa         >  5:82017/1‑149 (MQ=255)
                                                                                                                             ttCAGGAACGCGAGCAACGCAGTGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGAc        >  4:217856/1‑149 (MQ=255)
                                                                                                                                aGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCaa                                                         >  7:162761/1‑97 (MQ=255)
                                                                                                                                aGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCaa                                                         <  8:162761/97‑1 (MQ=255)
                                                                                                                                aGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGa         >  6:266304/1‑145 (MQ=255)
                                                                                                                                aGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGa         <  5:266304/145‑1 (MQ=255)
                                                                                                                                 ggAACGCGAGCAACGCAGTGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGtt    <  3:217856/149‑1 (MQ=255)
                                                                                                                                 ggAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGtt    <  4:108914/149‑1 (MQ=255)
                                                                                                                                  gAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGttt   <  7:57675/149‑1 (MQ=255)
                                                                                                                                     gccgAGCAGCGCGCTGCGTTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAg                             >  2:52036/3‑120 (MQ=255)
                                                                                                                                     gccgAGCAGCGCGCTGCGTTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAg                             <  1:52036/118‑1 (MQ=255)
                                                                                                                                       cgAGCAGCGCGCTGCGTTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTc  <  5:50671/145‑1 (MQ=255)
                                                                                                                                       cgAGCAGCGCGCTGCGTTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTc  >  6:50671/1‑145 (MQ=255)
                                                                                                                                                 |                                                                                                                                       
GATGGCGGCGGAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTC  >  NC_002947/1051824‑1052103

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