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_CP009273	3,366,821	0	T	C	100.0%	140.3 / NA	39	I39V (ATC→GTC)	nanA	N‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/21);  total (18/21)

CAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGC  >  NZ_CP009273/3366691‑3366955                                                                                                                                   |                                                                                                                                       cAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCctgct                                                                                                                                <  2:709281/139‑1 (MQ=255) cAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCctgct                                                                                                                                <  1:445502/139‑1 (MQ=255)            ggCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATtgaa                                                                                                                      <  1:606158/138‑1 (MQ=255)                 tGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATtgaa                                                                                                                      >  2:546898/1‑133 (MQ=255)                 tGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATtgaa                                                                                                                      <  1:546898/133‑1 (MQ=255)                             tACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACg                                                                                                       >  2:639182/1‑136 (MQ=255)                               ccTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGa                                                                                                   <  2:400447/138‑1 (MQ=255)                               ccTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGa                                                                                                   <  1:639182/138‑1 (MQ=255)                               ccTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGAc                                                                                                  >  2:143030/1‑139 (MQ=255)                                     gCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCtt                                                                                            <  1:761079/139‑1 (MQ=255)                                     gCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCt                                                                                             <  1:409939/138‑1 (MQ=255)                                            cGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATc                                                                                        <  2:475466/136‑1 (MQ=255)                                            cGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCa                                                                                      <  1:842122/138‑1 (MQ=255)                                                         cAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCttgttgttg                                                                         <  1:443898/138‑1 (MQ=255)                                                                  ttCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCaa                                                                    >  1:411413/1‑134 (MQ=255)                                                                  ttCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGATTGTTGTTGGTCaa                                                                    <  2:411413/134‑1 (MQ=255)                                                                        cTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAg                                                          >  2:16835/1‑138 (MQ=255)                                                                            gAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTg                                                     >  2:691075/1‑139 (MQ=255)                                                                            gAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTg                                                     <  1:143030/139‑1 (MQ=255)                                                                             aaaGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGc                                                    <  1:754631/139‑1 (MQ=255)                                                                                     tGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTa                                            >  2:535927/1‑139 (MQ=255)                                                                                           aaGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCAc                                      >  1:604776/1‑139 (MQ=255)                                                                                                 tCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGAc                                                                                                  >  2:736145/1‑73 (MQ=255)                                                                                                 tCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGAc                                                                                                  <  1:736145/73‑1 (MQ=255)                                                                                                 tCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCaaaa                                                                  <  1:58844/105‑1 (MQ=255)                                                                                                 tCGCCGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCaaaa                                                                  >  2:58844/1‑105 (MQ=255)                                                                                                    ccGGTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCg                             <  2:604776/139‑1 (MQ=255)                                                                                                      ggTCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGtt                           <  1:273412/139‑1 (MQ=255)                                                                                                        tCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTg                          <  2:45818/138‑1 (MQ=255)                                                                                                        tCGAACCACCCACGTATAAACCGTCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTg                          >  1:45818/1‑138 (MQ=255)                                                                                                                               tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATaaa                   <  1:718508/122‑1 (MQ=255)                                                                                                                               tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATaaa                   >  1:488803/1‑122 (MQ=255)                                                                                                                               tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATaaa                   <  2:488803/122‑1 (MQ=255)                                                                                                                               tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATaaa                   >  2:718508/1‑122 (MQ=255)                                                                                                                               tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGc  >  2:414033/1‑139 (MQ=255)                                                                                                                               tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGc  >  1:468268/1‑139 (MQ=255)                                                                                                                               tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGc  >  1:423192/1‑139 (MQ=255)                                                                                                                               tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGc  >  1:19093/1‑139 (MQ=255)                                                                                                                               tCGACGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGc  >  1:147888/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                       CAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGC  >  NZ_CP009273/3366691‑3366955

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