Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,527,151 T→C Q87Q (CAA→CAG nanS ← 9‑O‑acetyl‑N‑acetylneuraminic acid deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,527,1510TC90.4% 137.0 / 8.2 53Q87Q (CAA→CAGnanS9‑O‑acetyl‑N‑acetylneuraminic acid deacetylase
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (22/25);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

ACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAA  >  NZ_CP009273/4527065‑4527236
                                                                                      |                                                                                     
aCGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCtga                                                                                    <  2:630713/90‑1 (MQ=255)
             cTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCa                                                                       >  2:557400/1‑90 (MQ=255)
              tGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCag                                                                      <  2:699193/90‑1 (MQ=255)
              tGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCag                                                                      <  2:343893/90‑1 (MQ=255)
              tGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCag                                                                      <  2:240670/90‑1 (MQ=255)
              tGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCag                                                                      <  2:16331/90‑1 (MQ=255)
              tGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCag                                                                      <  2:157397/90‑1 (MQ=255)
                 aTTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAgg                                                                   >  2:643138/1‑90 (MQ=255)
                  ttATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGa                                                                  <  2:827445/90‑1 (MQ=255)
                   tATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGAt                                                                 >  2:28156/1‑90 (MQ=255)
                    aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                >  1:338389/1‑90 (MQ=255)
                         gAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGAt                                                           <  1:180370/90‑1 (MQ=255)
                             gAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       <  1:452467/90‑1 (MQ=255)
                             gAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       <  1:602794/90‑1 (MQ=255)
                                  gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  1:554395/90‑1 (MQ=255)
                                  gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  2:602140/90‑1 (MQ=255)
                                  gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  1:648170/90‑1 (MQ=255)
                                  gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  1:723961/90‑1 (MQ=255)
                                    aGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                <  1:143561/90‑1 (MQ=255)
                                     gTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCatat                                               >  1:454476/1‑90 (MQ=255)
                                               tGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATc                                     <  2:414270/90‑1 (MQ=255)
                                                   atatGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCgtgt                                 <  1:255964/90‑1 (MQ=255)
                                                   atatGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCgtgt                                 <  1:26489/90‑1 (MQ=255)
                                                         aGTGCCTGGCCAACGGTGCAGTACTGTGTATGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAg                           <  2:377384/90‑1 (MQ=255)
                                                          gTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATc                                     >  1:330483/1‑79 (MQ=255)
                                                          gTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATc                                     <  2:330483/79‑1 (MQ=255)
                                                          gTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGt                          >  2:813894/1‑90 (MQ=255)
                                                           tGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                <  1:723943/51‑1 (MQ=255)
                                                           tGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                >  2:723943/1‑51 (MQ=255)
                                                           tGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTg                         >  2:379657/1‑90 (MQ=255)
                                                             ccTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                >  1:688794/1‑65 (MQ=255)
                                                             ccTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                <  2:688794/65‑1 (MQ=255)
                                                             ccTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAg                       <  1:628575/90‑1 (MQ=255)
                                                               tGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTc                     >  2:185983/1‑90 (MQ=255)
                                                               tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTc                     <  1:813894/90‑1 (MQ=255)
                                                               tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTc                     >  2:541262/1‑90 (MQ=255)
                                                                    aaCGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTgg                >  2:305781/1‑90 (MQ=255)
                                                                    aaCGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTgg                >  1:457545/1‑90 (MQ=255)
                                                                    aaCGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTgg                >  2:701556/1‑90 (MQ=255)
                                                                          gCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGaataat          >  1:347261/1‑90 (MQ=255)
                                                                          gCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGaataat          >  2:661134/1‑90 (MQ=255)
                                                                                cTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGtt    <  1:503804/90‑1 (MQ=255)
                                                                                 tgtgTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGGGAGTCAGTGGAATAATGTCGTTa   >  2:659376/1‑90 (MQ=255)
                                                                                  gtgtCTGATGATTCGTTGCCAGAGGGTGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  <  2:246070/90‑1 (MQ=255)
                                                                                  gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACa                                       >  2:588957/1‑53 (MQ=255)
                                                                                  gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACa                                       <  1:588957/53‑1 (MQ=255)
                                                                                  gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTc                     >  2:149131/1‑71 (MQ=255)
                                                                                  gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTc                     <  1:149131/71‑1 (MQ=255)
                                                                                  gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  2:602947/1‑90 (MQ=255)
                                                                                  gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  <  1:643138/90‑1 (MQ=255)
                                                                                  gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  1:689419/1‑90 (MQ=255)
                                                                                  gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  2:346647/1‑90 (MQ=255)
                                                                                  gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  1:309810/1‑90 (MQ=255)
                                                                                      |                                                                                     
ACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAA  >  NZ_CP009273/4527065‑4527236

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: