Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,857,508:1 +T 8.8% coding (428/558 nt) PP_4272 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,857,5081.T 8.8% 211.0 / 13.5 57coding (428/558 nt)PP_4272membrane protein
Reads supporting (aligned to +/- strand):  ref base . (26/26);  new base T (2/3);  total (28/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 = 9.57e-01

CCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAAT  >  NC_002947/4857364‑4857647
                                                                                                                                                 |                                                                                                                                           
ccAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCgctg‑gctg                                                                                                                                         >  2:157557/1‑149 (MQ=255)
  aaCTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCa                                                                                                                                       >  5:271665/1‑149 (MQ=255)
     tCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATg                                                                                                                                     >  7:145234/1‑148 (MQ=255)
       ccTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGaa                                                                                                                                   <  2:107165/148‑1 (MQ=255)
       ccTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAg                                                                                                                                  <  5:243005/149‑1 (MQ=255)
             gTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCTAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAgc                                                                                                                            <  5:294755/149‑1 (MQ=255)
             gTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTTGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAgc                                                                                                                            <  2:170079/149‑1 (MQ=255)
                acacTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGTCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAgcggc                                                                                                                         <  8:71500/149‑1 (MQ=255)
                  acTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTAGGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgatgg                                                                                                                       <  8:101883/149‑5 (MQ=255)
                  acTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGgcgc                                                                                                                       >  5:169733/1‑149 (MQ=255)
                            gTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGt                                                                                                              <  7:198051/148‑1 (MQ=255)
                            gTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGt                                                                                                              >  8:198051/1‑148 (MQ=255)
                                  ccGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑CCTGCATGAAGCACAGCGCCGCCACGGCTGTacccaca                                                                                                       >  6:89392/1‑149 (MQ=255)
                                   cGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTacccacac                                                                                                      <  6:169733/149‑1 (MQ=255)
                                        tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgatgg                                                                                                                       <  2:66906/127‑5 (MQ=255)
                                        tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgatgg                                                                                                                       >  1:66906/1‑123 (MQ=255)
                                        tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgat                                                                                                                         <  8:191381/125‑3 (MQ=255)
                                        tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgat                                                                                                                         >  7:191381/1‑123 (MQ=255)
                                        tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCAt                                                                                                                                      >  3:340204/1‑112 (MQ=255)
                                        tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCAt                                                                                                                                      <  4:340204/112‑1 (MQ=255)
                                        tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTacccacacccac                                                                                                  >  8:12830/1‑148 (MQ=255)
                                        tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTacccacacccac                                                                                                  >  5:146449/1‑148 (MQ=255)
                                            gAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCAc                                                                                                                    <  6:115831/126‑1 (MQ=255)
                                            gAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCAc                                                                                                                    >  5:115831/1‑126 (MQ=255)
                                            gAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAacta                                                                                             <  1:157557/149‑1 (MQ=255)
                                                  ccccACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAg                                                                                       <  8:145234/149‑1 (MQ=255)
                                                      acacCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTaa                                                                                   <  7:266652/149‑1 (MQ=255)
                                                             acacAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAagag                                                                              >  7:167119/1‑147 (MQ=255)
                                                             acacAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTa                                                                            <  8:167119/149‑1 (MQ=255)
                                                                  aaaTGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGa                                                                       >  7:83555/1‑149 (MQ=255)
                                                                   aaTGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACgg                                                                                                                  >  4:185113/1‑105 (MQ=255)
                                                                         tGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAAc                                                                <  6:146449/149‑1 (MQ=255)
                                                                         tGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAAc                                                                <  6:271665/149‑1 (MQ=255)
                                                                         tGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAAc                                                                >  2:161266/1‑149 (MQ=255)
                                                                           gCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAAcc                                                               <  2:291858/148‑1 (MQ=255)
                                                                           gCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGAAGGTAAAAGAGTAGGCGCCCCCAACCa                                                              >  2:305861/1‑149 (MQ=255)
                                                                            ccAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACgg                                                                                                                  <  8:254132/96‑1 (MQ=255)
                                                                            ccAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACgg                                                                                                                  >  7:254132/1‑96 (MQ=255)
                                                                                   ctTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTccc                                                      >  7:83840/1‑149 (MQ=255)
                                                                                       aGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTAc                                                  <  8:83555/149‑1 (MQ=255)
                                                                                       aGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTAc                                                  >  2:90099/1‑149 (MQ=255)
                                                                                           aGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATc                                              <  7:49865/149‑1 (MQ=255)
                                                                                            gTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATca                                             <  8:83840/149‑1 (MQ=255)
                                                                                                cgcgCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACaa                                          <  1:161266/148‑1 (MQ=255)
                                                                                                            ggcgtcCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATcgc                             >  2:74877/1‑149 (MQ=255)
                                                                                                                tcCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGaaa                         <  7:12830/149‑1 (MQ=255)
                                                                                                                   tGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGCCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCAACGATCACAACCAGCAGCATCGCGaaaacc                      >  5:245851/1‑147 (MQ=255)
                                                                                                                        gAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGcc                 >  1:339110/1‑149 (MQ=255)
                                                                                                                           aaGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACCACTACTGCAGGTAAAAGCGTAGGCGACCGCACCCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGt              >  6:367853/1‑149 (MQ=255)
                                                                                                                            aGAACCCTGCCCCCATCGCTG‑GCTGCATGAACCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAATGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTg             <  5:89392/149‑1 (MQ=255)
                                                                                                                                cccTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTaaaa                                                                                 >  4:4080/1‑77 (MQ=255)
                                                                                                                                cccTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTaaaa                                                                                 <  3:4080/77‑1 (MQ=255)
                                                                                                                                 ccTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTccc                                                      >  7:234932/1‑103 (MQ=255)
                                                                                                                                 ccTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTccc                                                      <  8:234932/103‑1 (MQ=255)
                                                                                                                                  cTGCCCCCATCGCTG‑CCTGCATGAAGCACAGCGGCGCCACGGCTTTACCCACACCCACAAGTACTGCAGGTAAAAGAGTAGGCTACCGCAACAAGCCCTCCCCTCCGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAg       <  1:74877/149‑1 (MQ=255)
                                                                                                                                   tGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCACCCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGt      >  7:165582/1‑149 (MQ=255)
                                                                                                                                       cccATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAAt  >  6:38575/1‑149 (MQ=255)
                                                                                                                                       cccATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAAt  >  6:284228/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
CCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAAT  >  NC_002947/4857364‑4857647

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