Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,855,539 T→G 13.2% I106M (ATT→ATG)  PP_4270 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,855,5390TG13.2% 144.6 / 12.8 53I106M (ATT→ATG) PP_4270hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/25);  new base G (4/3);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGGGCGCTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCACCGTGAG  >  NC_002947/4855392‑4855686
                                                                                                                                                   |                                                                                                                                                   
aGTGGGCGCTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTCCGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATtc                                                                                                                                                    <  2:193320/149‑1 (MQ=255)
 gTGGGCGCTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATtca                                                                                                                                                   >  5:124559/1‑149 (MQ=255)
        cTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATGAATCAgcgg                                                                                                                                            <  8:267000/149‑2 (MQ=255)
          cGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATGAATCAGCGGGa                                                                                                                                          <  5:14468/149‑1 (MQ=255)
                 aGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGACAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTg                                                                                                                                   >  5:191404/1‑149 (MQ=255)
                  gTTCCACGTCATGGCCCGCCTGCGACTGGAGGAGACCATGAGCAGCGGCGATATTGATGTGCGGTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTgg                                                                                                                                  <  4:238424/149‑1 (MQ=255)
                       aCGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGgc                                                                                                                             <  1:240642/149‑1 (MQ=255)
                               gcctgcctGCTCTTGGAGGAGGCCATGAGCAGAGGCGATTTTGTTGTGCGTTGCTAGTGCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGTAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCgg                                                                                                                     <  8:268009/149‑1 (MQ=255)
                               gCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCAACGTCAGTACCCAGTGGCGCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCgg                                                                                                                     <  1:63958/149‑1 (MQ=255)
                                cccGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGc                                                                                                                    <  1:129306/149‑1 (MQ=255)
                                 ccGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGAc                                                                                                                                         <  6:102219/127‑1 (MQ=255)
                                 ccGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGAc                                                                                                                                         >  5:102219/1‑127 (MQ=255)
                                  cgcctgcgcctgGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGa                                                                                                                  <  1:102239/149‑1 (MQ=255)
                                      tgcgcctgGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAg                                                                                                               >  6:270317/1‑148 (MQ=255)
                                                 gaggCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGcc                                                                                                   <  5:270317/149‑1 (MQ=255)
                                                    gCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATAAGCGAGACGTACTGGTTGGCCTCCCAGGAGACAGCACCCCGATGccccc                                                                                                >  7:89741/1‑149 (MQ=255)
                                                       aTGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGa                                                                                             <  6:119141/149‑1 (MQ=255)
                                                           gcagcGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGc                                                                                                                    >  7:161011/1‑122 (MQ=255)
                                                           gcagcGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGc                                                                                                                    <  8:161011/122‑1 (MQ=255)
                                                                 gcgATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCa                                                                                   >  6:222828/1‑149 (MQ=255)
                                                                   gATATTGATGTGCGTTGCTACGACAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCAcc                                                                                 <  6:191404/149‑1 (MQ=255)
                                                                          aTGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGc                                                                           >  7:143589/1‑148 (MQ=255)
                                                                                gTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGc                                                                    <  8:143589/149‑1 (MQ=255)
                                                                                      aCGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTAc                                                                                                                                     >  4:144062/1‑78 (MQ=255)
                                                                                      aCGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTAc                                                                                                                                     <  3:144062/78‑1 (MQ=255)
                                                                                         tCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTg                                                           <  6:25576/149‑1 (MQ=255)
                                                                                           aGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGa                                                                                                                  >  1:57226/1‑92 (MQ=255)
                                                                                           aGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGa                                                                                                                  <  2:57226/92‑1 (MQ=255)
                                                                                            gTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGc                                                        <  6:124559/149‑1 (MQ=255)
                                                                                               cccAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAAc                                                                                        <  6:37434/114‑1 (MQ=255)
                                                                                               cccAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTa                                                      >  7:30881/1‑148 (MQ=255)
                                                                                                ccAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTaca                                                    >  4:236821/1‑149 (MQ=255)
                                                                                                   gTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCt                                                 <  8:30881/149‑1 (MQ=255)
                                                                                                          aaCGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAAc                                          >  8:163041/1‑149 (MQ=255)
                                                                                                             gACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGATGGGCTACACCTATTCAACCCt                                       >  5:267045/1‑149 (MQ=255)
                                                                                                                cGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTa                                                      <  8:284215/131‑1 (MQ=255)
                                                                                                                cGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTa                                                      >  7:284215/1‑131 (MQ=255)
                                                                                                                    cAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCa                                >  1:275652/1‑149 (MQ=255)
                                                                                                                         gATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCActgct                           >  8:240738/1‑149 (MQ=255)
                                                                                                                           ttttGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCActgctg                          >  4:9690/1‑148 (MQ=255)
                                                                                                                             ttGTGTGTAGTATCACGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTTCCACCCCAGGCGATCTGCACGGTGCTGGTCTACACCTATTCAACCCTGCAACCACTGCTGGa                        <  8:89741/148‑1 (MQ=255)
                                                                                                                             ttGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGa                        <  7:163041/148‑1 (MQ=255)
                                                                                                                                   ggAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGc                 >  2:62764/1‑149 (MQ=255)
                                                                                                                                    gAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCActg                             >  3:304469/1‑136 (MQ=255)
                                                                                                                                    gAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCActg                             <  4:304469/136‑1 (MQ=255)
                                                                                                                                        gtcgcGCTGATGAATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGtt            >  2:131465/5‑149 (MQ=255)
                                                                                                                                         tcgcGCTGATGAATCAGCGGGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCt                                       <  7:80406/118‑1 (MQ=255)
                                                                                                                                         tcgcGCTGATGAATCAGCGGGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCt                                       >  8:80406/4‑121 (MQ=255)
                                                                                                                                         tcgcGCTGATGAATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTg           >  7:16525/4‑149 (MQ=255)
                                                                                                                                         tcgcGCTGATGAATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTg           >  7:11751/4‑149 (MQ=255)
                                                                                                                                           ccGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCa         >  1:141344/1‑149 (MQ=255)
                                                                                                                                            cGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCGCCCGGAGCAACTGCCACACCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCAc        <  5:223327/149‑1 (MQ=255)
                                                                                                                                            cGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTg           >  7:14938/1‑146 (MQ=255)
                                                                                                                                            cGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTg           <  8:14938/146‑1 (MQ=255)
                                                                                                                                                 aTTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCACCGTGa   <  3:9690/149‑1 (MQ=255)
                                                                                                                                                  ttCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCACCGTGAg  <  3:236821/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGTGGGCGCTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGGCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCACCGTGAG  >  NC_002947/4855392‑4855686

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