Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_LN832404	4,023,921	G→A	58.0%	intergenic (+220/+161)	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,921	0	G	A	58.0%	‑6.2 / 34.9	31	intergenic (+220/+161)	EK12AG100_RS19960/fadA	NAD(P)H‑flavin reductase/acetyl‑CoA C‑acyltransferase FadA
Reads supporting (aligned to +/- strand):  ref base G (7/6);  new base A (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 = 2.21e-01

TCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCA  >  NZ_LN832404/4023635‑4024212                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                      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)                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                      TCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCA  >  NZ_LN832404/4023635‑4024212

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: ‑