Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	214,743	A→C	Y174D (TAC→GAC)	proS ←	proline‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	214,743	0	A	C	100.0%	218.8 / NA	55	Y174D (TAC→GAC)	proS	proline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (29/26);  total (29/26)

ATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCG  >  NZ_CP009273/214658‑214829                                                                                      |                                                                                       aTAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGtt                                                                                    <  2:771035/90‑1 (MQ=255)    gAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTcc                                                                                 >  2:545634/1‑90 (MQ=255)         ggTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAg                                                                            >  2:375389/1‑90 (MQ=255)           tgtCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAggg                                                                          <  2:120030/90‑1 (MQ=255)           tgtCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAggg                                                                          <  2:1427336/90‑1 (MQ=255)           tgtCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAggg                                                                          <  2:1900511/90‑1 (MQ=255)           tgtCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAggg                                                                          <  1:568055/90‑1 (MQ=255)             tCGGCTTGTACGGCGCGGAAATCCAGCCCCATTCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGAt                                                                        <  2:913915/90‑1 (MQ=255)             tCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCTTGCAGGGAt                                                                        <  1:1697236/90‑1 (MQ=255)             tCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGAt                                                                        <  1:1089674/90‑1 (MQ=255)             tCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGAt                                                                        <  1:1245560/90‑1 (MQ=255)             tCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGAt                                                                        <  1:756639/90‑1 (MQ=255)                    gTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgag                                                                 <  2:597466/90‑1 (MQ=255)                    gTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgag                                                                 <  1:314512/90‑1 (MQ=255)                    gTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgag                                                                 >  2:196342/1‑90 (MQ=255)                    gTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgag                                                                 >  1:677454/1‑90 (MQ=255)                     tACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgaga                                                                <  2:1095774/90‑1 (MQ=255)                      aCGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                               <  2:675362/90‑1 (MQ=255)                      aCGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                               <  2:526999/90‑1 (MQ=255)                      aCGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                               <  2:617842/90‑1 (MQ=255)                         gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgaga                                                                <  1:266984/86‑1 (MQ=255)                         gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgaga                                                                >  2:266984/1‑86 (MQ=255)                         gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                            >  1:474334/1‑90 (MQ=255)                         gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                            >  2:91549/1‑90 (MQ=255)                         gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                            >  1:1172635/1‑90 (MQ=255)                         gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                            >  1:396242/1‑90 (MQ=255)                           gcgGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATg                                                          >  2:1266760/1‑90 (MQ=255)                              gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                               >  1:25344/1‑82 (MQ=255)                              gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                               <  2:25344/82‑1 (MQ=255)                                    cAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTa                                                 <  1:471580/90‑1 (MQ=255)                                          cATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATc                                           >  1:1112201/1‑90 (MQ=255)                                           aTGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                               >  1:1617621/1‑69 (MQ=255)                                           aTGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                               <  2:1617621/69‑1 (MQ=255)                                           aTGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCt                                          >  2:458497/1‑90 (MQ=255)                                           aTGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCt                                          >  1:1509056/1‑90 (MQ=255)                                            tGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCtt                                         <  2:1164317/90‑1 (MQ=255)                                            tGCGGCCGAAGATTTTGCTGTAGGCCGCATACAATGCATCGTCGGTTTCCTGCAGGGAATCCCCAAAAGTAAGGAACGCGTTAGCAttct                                         <  1:1648785/90‑4 (MQ=255)                                                gCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTtcat                                     >  2:1301612/1‑90 (MQ=255)                                                          ttGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTcg                           >  1:515669/1‑90 (MQ=255)                                                          ttGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTcg                           >  2:1301969/1‑90 (MQ=255)                                                          ttGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTcg                           >  2:344313/1‑90 (MQ=255)                                                          ttGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTcg                           >  2:51164/1‑90 (MQ=255)                                                               gTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGa                      >  2:411922/1‑90 (MQ=255)                                                                    ccGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCa                 >  2:1554036/1‑90 (MQ=255)                                                                      gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAAc                    >  1:703965/1‑85 (MQ=255)                                                                      gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAAc                    <  2:703965/85‑1 (MQ=255)                                                                      gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATg               >  1:1482908/1‑90 (MQ=255)                                                                       cATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGa              <  2:1509056/90‑1 (MQ=255)                                                                           cATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGcc          <  2:1638309/90‑1 (MQ=255)                                                                            aTTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCg         >  1:299170/1‑90 (MQ=255)                                                                             ttGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa        >  1:684650/1‑90 (MQ=255)                                                                              tGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGaa       >  1:1640633/1‑90 (MQ=255)                                                                               gCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGaaa      >  2:1506404/1‑90 (MQ=255)                                                                                 aTCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAcg    >  1:1859874/1‑90 (MQ=255)                                                                                  tcgtcgGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAcgc   <  2:1742818/90‑1 (MQ=255)                                                                                   cGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAcgcg  <  2:1112201/90‑1 (MQ=255)                                                                                      |                                                                                       ATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCG  >  NZ_CP009273/214658‑214829

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