Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,537,643 C→A 6.9% intergenic (‑49/+16) nanS ← / ← nanM 9‑O‑acetyl N‑acetylneuraminic acid esterase/N‑acetylneuraminic acid mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,537,6430CA 6.9% 200.7 / 5.2 58intergenic (‑49/+16)nanS/nanM9‑O‑acetyl N‑acetylneuraminic acid esterase/N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base C (22/32);  new base A (2/2);  total (24/34)
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.12e-01

TGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTA  >  NC_000913/4537498‑4537784
                                                                                                                                                 |                                                                                                                                             
tGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCgg                                                                                                                                            <  2:216519/149‑1 (MQ=255)
tGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCgg                                                                                                                                            <  1:424020/149‑1 (MQ=255)
   ccGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGc                                                                                                                                         <  3:427684/149‑1 (MQ=255)
     ggTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGccc                                                                                                                                       >  2:402101/1‑149 (MQ=255)
     ggTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGccc                                                                                                                                       <  3:256029/149‑1 (MQ=255)
      gTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCt                                                                                                                                      <  1:402101/149‑1 (MQ=255)
         aTGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGaa                                                                                                                                   <  4:515241/149‑1 (MQ=255)
                  cTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGtt                                                                                                                          >  1:28163/1‑149 (MQ=255)
                  cTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGtt                                                                                                                          <  2:28163/149‑1 (MQ=255)
                  cTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGtt                                                                                                                          <  1:98702/149‑1 (MQ=255)
                  cTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGtt                                                                                                                          <  1:396322/149‑1 (MQ=255)
                       ccATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTa                                                                                                                     <  3:545468/149‑1 (MQ=255)
                         aTACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACt                                                                                                                   >  3:56903/1‑149 (MQ=255)
                           aCGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACtgt                                                                                                                 <  1:564351/149‑1 (MQ=255)
                              ccATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGAc                                                                                                              >  2:123433/1‑149 (MQ=255)
                              ccATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGAc                                                                                                              >  4:234994/1‑149 (MQ=255)
                                 tGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACttt                                                                                                           >  3:480667/1‑149 (MQ=255)
                                  ggCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTtta                                                                                                          <  3:471328/149‑1 (MQ=255)
                                     aTTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTttatta                                                                                                       >  1:458531/1‑149 (MQ=255)
                                     aTTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGccc                                                                                                                                       >  2:195323/1‑117 (MQ=255)
                                     aTTAGACTGACCAGCAACGGTAAGAACAGAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGccc                                                                                                                                       <  1:195323/117‑1 (MQ=255)
                                      ttAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTtcttat                                                                                                      >  4:382348/1‑149 (MQ=255)
                                       tAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATc                                                                                                     <  1:339114/149‑1 (MQ=255)
                                            tGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACttt                                                                                                           <  4:71576/138‑1 (MQ=255)
                                            tGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACttt                                                                                                           >  3:71576/1‑138 (MQ=255)
                                            tGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTca                                                                                                <  2:356877/149‑1 (MQ=255)
                                            tGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTca                                                                                                <  3:373253/149‑1 (MQ=255)
                                              aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACtgt                                                                                                                 <  4:242283/130‑1 (MQ=255)
                                              aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACtgt                                                                                                                 >  3:242283/1‑130 (MQ=255)
                                              aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTcaca                                                                                              >  1:190438/1‑149 (MQ=255)
                                              aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTcaca                                                                                              <  2:458531/149‑1 (MQ=255)
                                                  gCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGa                                                                                          >  2:308930/1‑149 (MQ=255)
                                                    aaCGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATc                                                                                        >  3:357868/1‑149 (MQ=255)
                                                        gTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCaaaa                                                                                    <  4:357868/149‑1 (MQ=255)
                                                        gTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGCCCTgct                                                                                                                                   >  1:235121/1‑100 (MQ=38)
                                                        gTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGCCCTgct                                                                                                                                   <  2:235121/102‑3 (MQ=38)
                                                              aCATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAAt                                                                              <  4:129520/149‑1 (MQ=255)
                                                                                gATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTgccgc                                                            <  1:123433/149‑1 (MQ=255)
                                                                                gATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTgccgc                                                            <  3:234994/149‑1 (MQ=255)
                                                                                  tattatTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCt                                                          >  3:101621/1‑149 (MQ=255)
                                                                                       ttGCGTTCATATAGATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGt                                                     <  1:32804/149‑1 (MQ=255)
                                                                                           gTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCa                                                 >  1:299320/1‑149 (MQ=255)
                                                                                            ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTg                                                                                                                       <  3:62795/78‑1 (MQ=255)
                                                                                            ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTg                                                                                                                       >  4:62795/1‑78 (MQ=255)
                                                                                                      tttCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGCGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCaataat                                      >  4:583938/1‑149 (MQ=255)
                                                                                                       ttCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATc                                     <  3:382348/149‑1 (MQ=255)
                                                                                                         cgcgTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCaa                                   <  1:517744/149‑1 (MQ=255)
                                                                                                          gcgTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAAt                                  <  3:583938/149‑1 (MQ=255)
                                                                                                                   aaaTGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatta                         <  4:101621/149‑1 (MQ=255)
                                                                                                                   aaaTGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatta                         >  4:473313/1‑149 (MQ=255)
                                                                                                                         cAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCaataat                                      <  4:104625/130‑1 (MQ=255)
                                                                                                                         cAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCaataat                                      >  3:104625/1‑130 (MQ=255)
                                                                                                                          aaaTTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAgg                  >  1:394156/1‑149 (MQ=255)
                                                                                                                          aaaTTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCCGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAgg                  <  4:339733/149‑1 (MQ=255)
                                                                                                                               gATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAAt             <  1:308930/149‑1 (MQ=255)
                                                                                                                                       ggCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTcc     >  3:216030/1‑149 (MQ=255)
                                                                                                                                         ccccTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGt   >  3:562979/1‑149 (MQ=255)
                                                                                                                                          cccTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTa  <  2:528967/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
TGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTA  >  NC_000913/4537498‑4537784

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: