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,1510TC96.7% 196.0 / ‑1.1 60Q87Q (CAA→CAGnanS9‑O‑acetyl‑N‑acetylneuraminic acid deacetylase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (28/30);  total (29/31)
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.72e-01

GAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAA  >  NZ_CP009273/4527072‑4527236
                                                                               |                                                                                     
gAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCg                                                                             <  1:548316/90‑1 (MQ=255)
      cTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCa                                                                       >  2:193971/1‑90 (MQ=255)
        gCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCaga                                                                     >  2:155938/1‑90 (MQ=255)
           ttATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGa                                                                  <  2:34355/90‑1 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATg                                                                <  1:193971/90‑1 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                >  2:94394/1‑90 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                <  2:920005/90‑1 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                <  2:699903/90‑1 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                >  2:62702/1‑90 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                <  1:210414/90‑1 (MQ=255)
              tCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATgg                                                               >  2:406409/1‑90 (MQ=255)
               cAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGt                                                              >  1:46655/1‑90 (MQ=255)
                  gAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGAt                                                           >  1:383458/1‑90 (MQ=255)
                  gAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGAt                                                           <  2:352135/90‑1 (MQ=255)
                      gAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       <  1:433686/90‑1 (MQ=255)
                      gAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       >  1:296877/1‑90 (MQ=255)
                      gAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       <  2:383458/90‑1 (MQ=255)
                       aaGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAccc                                                      >  1:424718/1‑90 (MQ=255)
                           gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  2:179257/90‑1 (MQ=255)
                           gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  2:27758/90‑1 (MQ=255)
                           gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  2:46655/90‑1 (MQ=255)
                           gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  2:664979/90‑1 (MQ=255)
                             aGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                <  2:604468/90‑1 (MQ=255)
                             aGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                <  1:754535/90‑1 (MQ=255)
                             aGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                <  1:505019/90‑1 (MQ=255)
                                 aTTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCt                                            <  1:444099/90‑1 (MQ=255)
                                   ttCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCag                                                                      <  1:617121/62‑1 (MQ=255)
                                   ttCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCag                                                                      >  2:617121/1‑62 (MQ=255)
                                   ttCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGa                                          >  1:607922/1‑90 (MQ=255)
                                      cGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACa                                       <  1:406409/90‑1 (MQ=255)
                                        tGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATc                                     >  1:453368/1‑90 (MQ=255)
                                            atatGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACa                                       <  1:255420/84‑1 (MQ=255)
                                            atatGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACa                                       >  2:255420/1‑84 (MQ=255)
                                            atatGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCgtgt                                 <  1:8312/90‑1 (MQ=255)
                                            atatGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCgtgt                                 >  2:206312/1‑90 (MQ=255)
                                                   gTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGt                          >  1:783411/1‑90 (MQ=255)
                                                   gTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGt                          >  1:244826/1‑90 (MQ=255)
                                                   gTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGt                          <  1:301123/90‑1 (MQ=255)
                                                    tGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCNAGAGGATGGTGATAACNCTTNATANCCTTAACATCGTGTGGGCAGTg                         >  2:224631/1‑90 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGa                                                                  >  2:700424/1‑45 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGa                                                                  <  1:700424/45‑1 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       >  1:430781/1‑56 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       <  2:430781/56‑1 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTc                     >  1:283441/1‑90 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTc                     <  1:323105/90‑1 (MQ=255)
                                                                 gTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  1:665345/52‑1 (MQ=255)
                                                                 gTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  >  2:665345/1‑52 (MQ=255)
                                                                   gCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGaataat          >  2:468991/1‑90 (MQ=255)
                                                                   gCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGaataat          >  2:756420/1‑90 (MQ=255)
                                                                     cGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGt        <  2:424718/90‑1 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                <  2:854943/44‑1 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                >  2:864935/1‑44 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                >  1:854943/1‑44 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                <  1:864935/44‑1 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  1:902688/1‑90 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  1:117313/1‑90 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  2:847137/1‑90 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  <  1:468991/90‑1 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  1:778736/1‑90 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  <  2:117313/90‑1 (MQ=255)
                                                                               |                                                                                     
GAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAA  >  NZ_CP009273/4527072‑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: