Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,618,706 A→T 6.9% Y203N (TAC→AAC)  PP_2292 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,618,7060AT 6.9% 166.5 / 3.3 58Y203N (TAC→AAC) PP_2292hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/24);  new base T (2/2);  total (32/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 = 8.66e-01

ATTTGTTTCTAGCCTCGCTTTCGTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACGCAGCATCGATTCCAC  >  NC_002947/2618568‑2618853
                                                                                                                                          |                                                                                                                                                    
aTTTGTTTCTAGCCTCGCTTTCGTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAgc                                                                                                                                            >  8:4409/1‑149 (MQ=255)
      ttCTAGCCTCGCTTTCGTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTg                                                                                                                                      >  5:251567/1‑149 (MQ=255)
           gCCTCGCTTTCGTTTTTGAACGTTAAGTAACACATCGCAAGCGCTTTGCTGACTATTTACTCTCGCGGAGGAAATGGCTGGAATGTATAAATGAGTCACAATATATCGACTGCGGCCAATTTCTTGTAGACGTCCTGCGCACTGCCGa                                                                                                                                  >  4:179834/1‑148 (MQ=255)
            ccTCGCTTTCGTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAAt                                                                                                                                >  3:265396/1‑149 (MQ=255)
              tCGCTTTCGTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAAtt                                                                                                                               >  3:232765/1‑148 (MQ=255)
                  tttCGTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCa                                                                                                                                         <  1:24813/134‑1 (MQ=255)
                  tttCGTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCa                                                                                                                                         >  2:24813/1‑134 (MQ=255)
                   ttCGTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTTGAAGTCCAGCGCACTGCCGAATTTAgtgt                                                                                                                         >  5:10799/1‑149 (MQ=255)
                   ttCGTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTg                                                                                                                                      >  7:203394/1‑136 (MQ=255)
                   ttCGTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTg                                                                                                                                      <  8:203394/136‑1 (MQ=255)
                    tCGTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAgtg                                                                                                                          <  6:251567/147‑1 (MQ=255)
                      gTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTTGAAGTCCAGCGCACTGCCGAATTTAGTGTCAg                                                                                                                      >  5:86455/1‑149 (MQ=255)
                                gTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGcc                                                                                                           <  6:126247/149‑1 (MQ=255)
                                   aaGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCAt                                                                                                         <  7:4409/148‑1 (MQ=255)
                                            aTCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGc                                                                                                >  6:55551/1‑148 (MQ=255)
                                                cAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTg                                                                                           <  6:28765/149‑1 (MQ=255)
                                                     gcTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCAt                                                                                                         <  3:138093/130‑1 (MQ=255)
                                                     gcTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCAt                                                                                                         >  4:138093/1‑130 (MQ=255)
                                                       tttGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAg                                                                                    >  2:229008/1‑149 (MQ=255)
                                                           ctctctATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGtgc                                                                                 <  4:232765/148‑1 (MQ=255)
                                                           ctatctATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGCGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGtgct                                                                                <  4:7147/146‑1 (MQ=255)
                                                                tATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGt                                                                           >  7:297336/1‑149 (MQ=255)
                                                                tATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGt                                                                           >  5:101312/1‑149 (MQ=255)
                                                                tATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGt                                                                           >  2:41346/1‑149 (MQ=255)
                                                                  ttGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGATTGCTGCCAATGTATTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCAATGAGTTACGCCATTGCGATTAGTGCTTCCGTTg                                                                         >  1:94018/1‑149 (MQ=255)
                                                                    gACTCTCGCGGATGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTATAACTCCAGCGCACTGGCGAATTTAGTGTCAGTAAAATCGCCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTg                                                                         <  1:229008/148‑1 (MQ=255)
                                                                     aCTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTc                                                                      <  5:55551/149‑1 (MQ=255)
                                                                       tctcGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTcc                                                                     >  4:188799/1‑147 (MQ=255)
                                                                           gcgGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTTGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACatat                                                                <  6:86455/149‑1 (MQ=255)
                                                                           gcgGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACata                                                                 <  1:41346/148‑1 (MQ=255)
                                                                             ggaggaAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATAttt                                                              <  8:297336/149‑1 (MQ=255)
                                                                                ggaAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTTGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTc                                                           <  6:10799/149‑1 (MQ=255)
                                                                                  aaaTGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTc                                                           >  5:267134/1‑147 (MQ=255)
                                                                                  aaaTGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTc                                                           <  6:267134/147‑1 (MQ=255)
                                                                                        cTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGc                                                   <  4:265396/149‑1 (MQ=255)
                                                                                                tAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCtata                                               <  4:209414/145‑1 (MQ=255)
                                                                                                tAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCAGTTGGTCACATATTTGTCGGAGTCGCtata                                               >  3:209414/1‑145 (MQ=255)
                                                                                                  gAAATGAGTCACAATATATCGACTGCTGCCAATGTCATGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAAttt                                         <  2:94018/149‑1 (MQ=255)
                                                                                                               atatatCGTCTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGATTTACGCCATTTCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAAtttt                            <  3:179834/149‑1 (MQ=255)
                                                                                                                tatatCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAAttttt                           <  8:169324/149‑1 (MQ=255)
                                                                                                                tatatCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAAttttt                           >  7:169324/1‑149 (MQ=255)
                                                                                                                    tCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTgg                         >  3:175801/1‑147 (MQ=255)
                                                                                                                    cgatctgctgCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTgga                        <  3:38364/144‑1 (MQ=255)
                                                                                                                           ctgcCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTgga                        >  4:38364/1‑141 (MQ=255)
                                                                                                                           ctgcCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACgc                <  2:339369/149‑1 (MQ=255)
                                                                                                                            tgcCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCg                                                          >  4:112633/1‑106 (MQ=255)
                                                                                                                            tgcCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCg                                                          <  3:112633/106‑1 (MQ=255)
                                                                                                                            tgcCAAGGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGGCATGAGTTACGCCATTGCGATTAGTGCGGCCGTTGGTCACATATTAGGCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTggaaggaagg                 >  6:255324/1‑145 (MQ=255)
                                                                                                                              ccAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACgcagc             <  4:175801/149‑1 (MQ=255)
                                                                                                                              ccAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACgcagc             <  7:203817/149‑1 (MQ=255)
                                                                                                                               cAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACgcagca            >  1:3319/1‑149 (MQ=255)
                                                                                                                               aaaTGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACgcagca            >  1:3322/2‑149 (MQ=255)
                                                                                                                                 aTGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACGCAGCATc          >  2:84421/1‑149 (MQ=255)
                                                                                                                                     cTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACGCAGCATCGAtt      >  3:293904/1‑149 (MQ=255)
                                                                                                                                      ttGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACGCAGCATCGATTc     >  3:206805/1‑149 (MQ=255)
                                                                                                                                        gTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACGCAGCATCGATTcc    >  8:90454/1‑148 (MQ=255)
                                                                                                                                         tAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACGCAGCATCGATTCCAc  >  7:229501/1‑149 (MQ=255)
                                                                                                                                         tAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACGCAGCATCGATTCCAc  >  2:281826/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
ATTTGTTTCTAGCCTCGCTTTCGTTTTTGCACGTTAAGTAACACATCGCAAGCGCTTTGCTCTCTATTGACTCTCGCGGAGGAAATGGCTGGAATGTAGAAATGAGGCACAATATATCGACTGCTGCCAATGTCTTGTAGAAGTCCAGCGCACTGCCGAATTTAGTGTCAGTAAAATCGCCATGAGTTACGCCATTGCGATTAGTGCTGCCGTTGGTCACATATTTGTCGGAGTCGCTATAGAATTTTTGGAGGAATTTTTGGAAGGAACGCAGCATCGATTCCAC  >  NC_002947/2618568‑2618853

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: