Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,051,975 C→T 8.5% Q464Q (CAG→CAA)  PP_0910 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,051,9750CT 8.5% 133.7 / 2.4 47Q464Q (CAG→CAA) PP_0910hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/24);  new base T (2/2);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

GGCGGAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTCAGGCCGATGCCCC  >  NC_002947/1051830‑1052120
                                                                                                                                                 |                                                                                                                                                  
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)
                                                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                   >  6:50671/1‑145 (MQ=255)
                                                                                                                                 cgAGCAGCGCGCTGCGTTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTc                   <  5:50671/145‑1 (MQ=255)
                                                                                                                                              gcgcTGATCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTCAGGCCGATGcccc  <  4:8644/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGCGGAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAC‑GGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTCAGGCCGATGCCCC  >  NC_002947/1051830‑1052120

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