Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,366,821 T→C I39V (ATC→GTC)  nanA ← N‑acetylneuraminate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,366,8210TC90.8% 182.3 / 11.2 65I39V (ATC→GTC) nanAN‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (26/33);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

CAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCT  >  NZ_CP009273/3366691‑3366956
                                                                                                                                  |                                                                                                                                       
cAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCctgct                                                                                                                                 <  1:748720/139‑1 (MQ=255)
cAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCctgct                                                                                                                                 <  2:397679/139‑1 (MQ=255)
cAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCctgct                                                                                                                                 <  1:276994/139‑1 (MQ=255)
      aCGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAatat                                                                                                                           <  1:998789/139‑1 (MQ=255)
           ggCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATtgaa                                                                                                                       <  2:707871/138‑1 (MQ=255)
           ggCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATtgaa                                                                                                                       <  2:28250/138‑1 (MQ=255)
                            tACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACg                                                                                                        >  2:1054227/1‑136 (MQ=255)
                             aCCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGa                                                                                                    >  1:157225/1‑139 (MQ=255)
                              ccTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGAc                                                                                                   >  1:134772/1‑139 (MQ=255)
                                   cGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCt                                                                                              >  2:1004154/1‑139 (MQ=255)
                                           cGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCa                                                                                       <  1:912540/138‑1 (MQ=255)
                                           cGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCa                                                                                       <  2:314515/138‑1 (MQ=255)
                                                             ccTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCaa                                                                     >  1:724907/1‑138 (MQ=255)
                                                                gTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAg                                                                  >  2:508782/1‑138 (MQ=255)
                                                                 ttCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAg                                                                  <  1:410024/137‑1 (MQ=255)
                                                                      gCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAg                                                           <  2:236885/139‑1 (MQ=255)
                                                                      gCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAg                                                           <  1:497392/139‑1 (MQ=255)
                                                                       cTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAg                                                           >  2:497392/1‑138 (MQ=255)
                                                                        tCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGcc                                                                                                                                      >  1:330905/1‑62 (MQ=255)
                                                                        tCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGcc                                                                                                                                      <  2:330905/62‑1 (MQ=255)
                                                                           gAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTg                                                      >  2:254291/1‑139 (MQ=255)
                                                                                    tGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGcc                                                                                                                                      <  1:538725/50‑1 (MQ=255)
                                                                                    tGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGcc                                                                                                                                      >  2:538725/1‑50 (MQ=255)
                                                                                    tGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCaaaa                                                                   >  1:906336/1‑117 (MQ=255)
                                                                                    tGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCaaaa                                                                   <  2:906336/117‑1 (MQ=255)
                                                                                    tGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGt                                                              <  1:618455/122‑1 (MQ=255)
                                                                                    tGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGt                                                              >  2:618455/1‑122 (MQ=255)
                                                                                    tGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGt                                                              <  1:49892/122‑1 (MQ=255)
                                                                                    tGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGt                                                              >  2:49892/1‑122 (MQ=255)
                                                                                    tGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGA‑‑ATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACg                                           >  2:1144326/1‑139 (MQ=255)
                                                                                      tACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACg                                           <  1:659914/139‑1 (MQ=255)
                                                                                        cAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGc                                          <  2:507293/138‑1 (MQ=255)
                                                                                            ggCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATcc                                                                                        >  2:583749/1‑88 (MQ=255)
                                                                                            ggCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATcc                                                                                        <  1:583749/88‑1 (MQ=255)
                                                                                                tCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAAt                                 >  1:1225369/1‑139 (MQ=255)
                                                                                                tCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGAc                                                                                                   <  1:789545/73‑1 (MQ=255)
                                                                                                tCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGAc                                                                                                   >  2:789545/1‑73 (MQ=255)
                                                                                                tCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCttgttgt                                                                            <  2:613613/96‑1 (MQ=255)
                                                                                                tCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCttgttgt                                                                            >  1:613613/1‑96 (MQ=255)
                                                                                                tCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGt                                                       <  1:891653/117‑1 (MQ=255)
                                                                                                tCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGt                                                       >  2:891653/1‑117 (MQ=255)
                                                                                                tCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAAt                                 >  1:472350/1‑139 (MQ=255)
                                                                                                tCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAAt                                 >  1:1106572/1‑139 (MQ=255)
                                                                                                   ccGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCg                              <  1:1013107/139‑1 (MQ=255)
                                                                                                   ccGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCg                              <  1:1054227/139‑1 (MQ=255)
                                                                                                   ccGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCg                              <  1:370610/139‑1 (MQ=255)
                                                                                                     ggTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGtt                            <  1:84864/139‑1 (MQ=255)
                                                                                                     ggTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGtt                            <  2:850372/139‑1 (MQ=255)
                                                                                                       tCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGc                          >  2:934751/1‑139 (MQ=255)
                                                                                                         gAACCACCCACGTATAAACCGTCGACGCCCTGCTGA‑‑ATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCAt                       <  1:1144326/138‑1 (MQ=255)
                                                                                                             cacccacGTATAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATaaa                    <  1:1004154/139‑1 (MQ=255)
                                                                                                             cacccacGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATaaa                    <  1:109097/139‑1 (MQ=255)
                                                                                                                    gTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAg                                                           >  1:600205/1‑93 (MQ=255)
                                                                                                                    gTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAg                                                           <  2:600205/93‑1 (MQ=255)
                                                                                                                         aaCCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGt        <  2:1106572/139‑1 (MQ=255)
                                                                                                                             gTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATg    <  2:525872/139‑1 (MQ=255)
                                                                                                                              tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGa                                                                                                    <  1:854679/42‑1 (MQ=255)
                                                                                                                              tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGa                                                                                                    >  2:854679/1‑42 (MQ=255)
                                                                                                                              tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCttgttgt                                                                            >  1:58997/1‑66 (MQ=255)
                                                                                                                              tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCttgttgt                                                                            <  2:58997/66‑1 (MQ=255)
                                                                                                                              tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGc   >  2:657998/1‑139 (MQ=255)
                                                                                                                              tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGc   >  2:373970/1‑139 (MQ=255)
                                                                                                                               cGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCg                              <  2:616301/111‑1 (MQ=255)
                                                                                                                               cGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCg                              >  1:616301/1‑111 (MQ=255)
                                                                                                                               cGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCt  >  2:417147/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                       
CAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCT  >  NZ_CP009273/3366691‑3366956

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑