Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_LN832404	4,023,922	T→G	58.1%	intergenic (+221/+160)	EK12AG100_RS19960 → / ← fadA	NAD(P)H‑flavin reductase/acetyl‑CoA C‑acyltransferase FadA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_LN832404	4,023,922	0	T	G	58.1%	‑4.2 / 36.6	31	intergenic (+221/+160)	EK12AG100_RS19960/fadA	NAD(P)H‑flavin reductase/acetyl‑CoA C‑acyltransferase FadA
Reads supporting (aligned to +/- strand):  ref base T (7/6);  new base G (14/4);  total (21/10)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.24e-01

TCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCG  >  NZ_LN832404/4023635‑4024215                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                        tCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAc                                                                                                                                                                                                                                                                                                   >  1:899503/1‑294 (MQ=255)  cTGTTTTGCAGTGAGCCTAGTGCGCGGGAAGATCGCCGGTTTGGCGATGCGTTTGCCTTTATCTTCGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCACCTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACa                                                                                                                                                                                                                                                                                                  <  1:279751/294‑1 (MQ=25)     ttttGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACattg                                                                                                                                                                                                                                                                                               >  1:324017/1‑291 (MQ=255)      tttGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACattgg                                                                                                                                                                                                                                                                                              >  1:796291/1‑290 (MQ=255)      ttggCAGGGAGCGTACGGCGCGGTAAGGGTGCTGGTTTGGTGGTGGGTTTGTAGTTATCTGAGATATGAAAAAACCCGCCCCTGACAGGCGTGAAGAACGGCAACTAAACTGTTAGTCAGTGGCATTTAGATCTATGACGTAGCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACattgg                                                                                                                                                                                                                                                                                              <  2:218361/291‑5 (MQ=18)                                   cGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGtt                                                                                                                                                                                                                                                                <  1:72316/294‑1 (MQ=18)                                       tGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTGTCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCAcg                                                                                                                                                                                                                                                            >  1:591008/1‑294 (MQ=18)                                           tGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCAcgccgc                                                                                                                                                                                                                                                        >  1:74365/1‑294 (MQ=18)                                            ggCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCa                                                                                                                                                                                                                                                       <  1:419790/294‑1 (MQ=18)                                             gCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACCAATAGTCCGAACCGTAGGCCGGATAAGCCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGCAGACCGGATACGGCGTTCACGCCGCAt                                                                                                                                                                                                                                                      >  1:173220/1‑294 (MQ=18)                                                  gCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                                                                                                                                 >  1:757098/1‑294 (MQ=17)                                                     tttGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGACTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGCCCGGCTAAGGCGTCCACGCTGCATCCGGCaag                                                                                                                                                                                                                                              >  1:50528/1‑293 (MQ=17)                                                              aTCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGttaaac                                                                                                                                                                                                                                     >  1:896115/1‑289 (MQ=17)                                                                 tGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGCTAAGGCGTTCACGCCGCATCCGGCAAGTGGttaaacccg                                                                                                                                                                                                                                  >  1:224694/1‑286 (MQ=14)                                                                 tGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTACTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGCTCAGGCGTTCACGCCGCATCCGGCAAGTGGttaaacccg                                                                                                                                                                                                                                  >  1:300035/1‑286 (MQ=14)                                                                  gagaTATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGCACGGTAGCCCGGACAAGGCGTTCCCGCCGCATCCGGCAAGTggctacacccgt                                                                                                                                                                                                                                 >  2:483878/1‑284 (MQ=17)                                                                      tataAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACCAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTCTCAGGCCTAC‑ATTGGTGCCGGATCGGTCGCCCGGATCAGGCGTTCACCCCCCATCCGGCCAGTGGttaaacccgctcac                                                                                                                                                                                                                             >  2:689046/1‑281 (MQ=14)                                                                      tataAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGttaaccccgctcaa                                                                                                                                                                                                                             >  1:130427/1‑281 (MQ=14)                                                                                                                                                                                                                                                     atctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACAACGTCGCAATACCCTGACCCAGAACGATACACATCGTCGCCAGCCCAAACTACACGTCTTTGCGTTCCATCAGa                                             >  1:666810/13‑294 (MQ=21)                                                                                                                                                                                                                                                           caaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCCCGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCCATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAACCTGAACGTCTTTGCTTCCCATCACATTcacc                                       >  2:417520/7‑292 (MQ=255)                                                                                                                                                                                                                                                            aaaccatgCCCGAGTCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTcagca                                      <  1:980872/289‑1 (MQ=255)                                                                                                                                                                                                                                                                 attcccgATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCAGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCgtg                                 <  1:986336/288‑1 (MQ=255)                                                                                                                                                                                                                                                                       gATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGCCCAAACTGCACGTCCTTTCGTTCCATCCGATTCAGCCGCGTGGTGCTg                           >  2:170813/1‑294 (MQ=255)                                                                                                                                                                                                                                                                        aTGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCTTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCCCAAACACCGTCGCCATACCCCGACCCAGACCGATACCAACCGCCGCCCGACCAAACTGAACGCCATTGCGTTCCTTCACATTCCGCCGCGTGGCGCTGa                          >  2:330600/1‑294 (MQ=21)                                                                                                                                                                                                                                                                         tGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGACCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGCCGCACTACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACAGAACGTCATTGCCTTCCATCAGATTCACCCGCGTGCTGCCGAt                         >  2:128163/1‑294 (MQ=25)                                                                                                                                                                                                                                                                          gCGACGCTGTCGCGTCTGATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATa                        <  1:915003/294‑1 (MQ=255)                                                                                                                                                                                                                                                                            gaggCTGTTGCGTTTTATCGGGCCTACAAATAGT‑CC‑GAACCGGAGGCCGGATAAGGCGTTTACGCCGCATCAGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATAcg                      <  2:510248/291‑1 (MQ=255)                                                                                                                                                                                                                                                                             aCGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCACATTCAACAGCgtgggtcgtctacgc                     >  2:854816/1‑283 (MQ=255)                                                                                                                                                                                                                                                                               gCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACgcgc                   <  1:307763/294‑1 (MQ=255)                                                                                                                                                                                                                                                                                    cgcgTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCgg              <  1:746721/294‑1 (MQ=255)                                                                                                                                                                                                                                                                                    agcgTCTTAGCGGGCGTACAAAGAGG‑CC‑GAACCGAGGGCGGGATAAGTGTGTGAGGGCGAATCCGGGAATGTGTGCATGATGCCTGATGGGACGCTGGCGCGTCGTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTGCACGCCGCAGCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCGTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCgg              <  2:1078348/293‑1 (MQ=255)                                                                                                                                                                                                                                                                                           tATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGCCCAAACTGAACGTCTTTGCGTTCCCTCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGACCAAcc       >  1:556579/1‑294 (MQ=255)                                                                                                                                                                                                                                                                                             tCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCa     <  1:1010096/294‑1 (MQ=255)                                                                                                                                                                                                                                                                                                tGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCg  >  1:188962/1‑294 (MQ=255)                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                        TCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCG  >  NZ_LN832404/4023635‑4024215

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