Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,537,406 A→G 8.1% I63I (ATT→ATC)  nanS ← 9‑O‑acetyl N‑acetylneuraminic acid esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,537,4060AG 8.1% 145.5 / 3.6 49I63I (ATT→ATC) nanS9‑O‑acetyl N‑acetylneuraminic acid esterase
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base G (2/2);  total (23/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 = 9.21e-01

CTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGA  >  NC_000913/4537261‑4537544
                                                                                                                                                 |                                                                                                                                          
cTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGaat                                                                                                                                          <  5:151351/148‑1 (MQ=255)
 tGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGaat                                                                                                                                          >  5:219525/1‑147 (MQ=255)
   cATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGaata                                                                                                                                         <  8:256765/146‑1 (MQ=255)
                gAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGt                                                                                                                                     <  7:245448/137‑1 (MQ=255)
                gAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGt                                                                                                                                     >  8:245448/1‑137 (MQ=255)
                   gggCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTAGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGGATAATGTCGTTAAAGTGACAt                                                                                                                       >  5:121398/1‑148 (MQ=255)
                          aaTTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGAc                                                                                                                         <  8:277120/139‑1 (MQ=255)
                          aaTTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGAc                                                                                                                         >  7:277120/1‑139 (MQ=255)
                          aaTTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGACGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAAACCTGCATATACTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGc                                                                                                                >  2:16338/1‑148 (MQ=255)
                             ttCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGGACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTc                                                                                                             >  5:215443/1‑148 (MQ=255)
                                cGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTTGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGAAATGGCGGGCCTCCg                                                                                                           <  4:264440/147‑1 (MQ=255)
                                cGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGGATAATGTCGTTAAAGTGACATGGCGGGCCTCCgg                                                                                                          >  3:288341/1‑148 (MQ=255)
                                             gTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATTGCGGGCCTCCGGGATGCGTATgcg                                                                                              >  1:247797/1‑147 (MQ=255)
                                                ccTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGGCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCaa                                                                                           <  2:222743/147‑1 (MQ=255)
                                                   ggCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGCGATGCGTATGCGCAAAtct                                                                                       >  2:152316/1‑148 (MQ=255)
                                                        aCGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGa                                                                                                                          >  8:23202/1‑108 (MQ=255)
                                                        aCGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGa                                                                                                                          <  7:23202/108‑1 (MQ=255)
                                                                gTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGAAGCGTATGCGCAAATCTCGCTAATTGttt                                                                           >  6:35964/1‑147 (MQ=255)
                                                                  aCTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAAt                                                                        >  3:219933/1‑148 (MQ=255)
                                                                  aCTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAAt                                                                        >  8:272909/1‑148 (MQ=255)
                                                                      tgtTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGGATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTg                                                                    <  4:288341/148‑1 (MQ=255)
                                                                      tgtTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTg                                                                    <  2:247797/148‑1 (MQ=255)
                                                                        tttGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGAAGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTggg                                                                  <  5:35964/148‑1 (MQ=255)
                                                                               atTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGt                                                                                                   >  4:214229/1‑108 (MQ=255)
                                                                               atTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGt                                                                                                   <  3:214229/108‑1 (MQ=255)
                                                                               atTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAgg                                                            <  4:106845/147‑1 (MQ=255)
                                                                                ttCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGgcg                                                          >  4:129171/1‑148 (MQ=255)
                                                                                  cGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTggg                                                                  >  3:89120/1‑138 (MQ=255)
                                                                                  cGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTggg                                                                  <  4:89120/138‑1 (MQ=255)
                                                                                   gTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCa                                                        <  7:272909/147‑1 (MQ=255)
                                                                                         agagGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTccc                                                 >  4:163412/1‑148 (MQ=255)
                                                                                          gagGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAAtctc                                                                                      >  1:60186/1‑110 (MQ=255)
                                                                                          gagGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAAtctc                                                                                      <  2:60186/110‑1 (MQ=255)
                                                                                           atgATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCAt                                                       <  6:10956/138‑1 (MQ=255)
                                                                                           agGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCAt                                                       >  5:10956/1‑140 (MQ=255)
                                                                                              aaggTGATAAACCTGCATATCCTTTAAATCGTGTGGGCAGTGAGGAAGTGGAATAATGTCCTTAAAGTGACAAGGCGGGCCACCGTGATGCGTATGCGCAAATCTCGCAAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTc                                             <  1:16338/145‑1 (MQ=255)
                                                                                                  tGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTa                                        <  3:129171/148‑1 (MQ=255)
                                                                                                         ccTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAg                                 <  1:152316/148‑1 (MQ=255)
                                                                                                                      aaCATCGTGTGGGCAGTGAGGCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCGCCGGGATGCGGATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCAt                      <  6:219525/146‑1 (MQ=255)
                                                                                                                           cGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTTCCTGTCCGGTa                                        <  7:214335/123‑1 (MQ=255)
                                                                                                                           cGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTTCCTGTCCGGTa                                        >  8:214335/1‑123 (MQ=255)
                                                                                                                            gtgtGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATg              >  2:193462/1‑148 (MQ=255)
                                                                                                                               tgGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGTCGGGCCTCCGGGATGCGTATGGGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCGTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCa           <  3:217343/148‑1 (MQ=255)
                                                                                                                                gggCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGc                                   <  7:153464/123‑1 (MQ=255)
                                                                                                                                gggCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGc                                   >  8:153464/1‑123 (MQ=255)
                                                                                                                                 ggCTGTGAGTCAGTGGAATAATGTCTTTAAAGTGTGATGGCGGGCCTCCGGGATGCGTATGCGCATATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCAtt         <  6:215443/148‑1 (MQ=255)
                                                                                                                                 ggCAGTGAGTCAGTGGGATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCTCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCAtt         <  6:121398/148‑1 (MQ=255)
                                                                                                                                 ggCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCAtt         <  4:219933/148‑1 (MQ=255)
                                                                                                                                       gAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGCGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACt    >  1:200277/1‑147 (MQ=255)
                                                                                                                                       gAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTa                                        >  5:65893/1‑111 (MQ=255)
                                                                                                                                       gAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTa                                        <  6:65893/111‑1 (MQ=255)
                                                                                                                                        aGTCAGTGGAATAATGTCGTTAAAGTGACAAGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGTATGAGGCGCAGCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGa  <  4:33075/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
CTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGA  >  NC_000913/4537261‑4537544

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