Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,857,517 A→C 9.4% F140C (TTC→TGC)  PP_4272 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,857,5170AC 9.4% 153.5 / 8.0 53F140C (TTC→TGC) PP_4272membrane protein
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base C (2/3);  total (26/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 = 8.68e-01

CCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAATCACCCATCCCCACGC  >  NC_002947/4857371‑4857662
                                                                                                                                                   |                                                                                                                                                 
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                                                                                                                                      >  1:66906/1‑123 (MQ=255)
                                 tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgatgg                                                                                                                                      <  2:66906/127‑5 (MQ=255)
                                 tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgat                                                                                                                                        >  7:191381/1‑123 (MQ=255)
                                 tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgat                                                                                                                                        <  8:191381/125‑3 (MQ=255)
                                 tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTacccacacccac                                                                                                                 >  8:12830/1‑148 (MQ=255)
                                 tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTacccacacccac                                                                                                                 >  5:146449/1‑148 (MQ=255)
                                     gAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCAc                                                                                                                                   >  5:115831/1‑126 (MQ=255)
                                     gAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCAc                                                                                                                                   <  6:115831/126‑1 (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                                                                 >  2:90099/1‑149 (MQ=255)
                                                                                aGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTAc                                                                 <  8:83555/149‑1 (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)
                                                                                                                                                 tGAAGCACAGCCGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGAAAAAGAGTAGGCGACCGCAACCAGCCCACCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACCCCGCTGGCGTAGTGAATCACCCATACCCACGc  >  5:203489/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAATCACCCATCCCCACGC  >  NC_002947/4857371‑4857662

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