Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_LN832404	4,023,962	T→C	37.5%	intergenic (+261/+120)	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,962	0	T	C	37.5%	28.0 / 22.0	24	intergenic (+261/+120)	EK12AG100_RS19960/fadA	NAD(P)H‑flavin reductase/acetyl‑CoA C‑acyltransferase FadA
Reads supporting (aligned to +/- strand):  ref base T (11/4);  new base C (4/5);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

CCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAAT‑AGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGAT  >  NZ_LN832404/4023671‑4024252                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                    ccTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATT‑GGTGCCGGATCGGTAGACCGGATAAGGCGTTCAcg                                                                                                                                                                                                                                                                                                  >  1:359295/1‑296 (MQ=18) ccTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATT‑GGTGCCGGATCGGTAGACCGGATAAGGCGTTCAcg                                                                                                                                                                                                                                                                                                  >  1:725445/1‑296 (MQ=18)   tGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATT‑GGTGCCGGATCGGTAGACCGGATAAGGCGTTCAcgcc                                                                                                                                                                                                                                                                                                >  1:360762/1‑296 (MQ=18)     tttGGCGATGCTATCGCATTTAGCTGGGATAGAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGGTATTCAGTGGCATTTAGATCTATGACGTATCTGGCGAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATT‑GGTGCCGGATCGGTAGACCGGATAAGGCGTTCAcgccgc                                                                                                                                                                                                                                                                                              <  2:449977/296‑1 (MQ=17)         gCGATGCGTTTGCATTTAGCTTAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACGAAACTGGTATTCAGTGGCAGTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCGGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑AGT‑GGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                                                                                                                                                                          <  1:479416/296‑1 (MQ=17)                gttgtcATTTATCTGGGATAACAAAAACCCCGCCCCTGACGTGCGGGAAGAACGGCAACTAACCTGTGATTCAGTGGCATTTAGATCTATGACGTACCTGGCAAACCATGCCCGATGCGACGCGGTCGCGTCTGATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGAGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATT‑GGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCaagt                                                                                                                                                                                                                                                                                   <  2:180999/291‑3 (MQ=14)                  ttGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATT‑GGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTgg                                                                                                                                                                                                                                                                                 >  1:796293/1‑296 (MQ=17)                  gtGCATTTAACTGGGATATAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATT‑GGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTgg                                                                                                                                                                                                                                                                                 <  1:111429/295‑1 (MQ=17)                         tATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGGGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATT‑GGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGttaaacc                                                                                                                                                                                                                                                                          <  1:208737/296‑7 (MQ=14)                           tCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCAATTGGTGCATGATGCCTGATGCGACGCCGGCGCGTCTTCTCAGGCCTAC‑ATT‑GGTACCGGCTCGGTAGACCGGCAAAGGCGTTCCCGCCCCATCCGGCaagcggttaaacccg                                                                                                                                                                                                                                                                        >  2:293229/1‑283 (MQ=17)                            cTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTTGTGC‑CG‑GCTCGGTAGACCGGATCAGGCGGTTCCGCCGCATCCGGCCAGTGGTTCaccccgc                                                                                                                                                                                                                                                                       >  2:35177/1‑290 (MQ=17)                                                                                                                                                                                                               gtatcagtcaaaccatgCCCGATGCGACGCTGTCGCGTCTGTTCGGGCCTACAAAT‑AGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCAGCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGa                                                                                   <  2:417838/282‑1 (MQ=255)                                                                                                                                                                                                                   ctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAAT‑AGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCCCAACACCCTCACCCCGACCGATACACATCCTCCCCCGACCCACCTGACCGGCTTTGCCTTCCCCCACATTcc                                                                               >  2:351348/11‑295 (MQ=18)                                                                                                                                                                                                                        aaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAAT‑AGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCagcagc                                                                          >  1:433954/6‑296 (MQ=255)                                                                                                                                                                                                                              tgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAAT‑AGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGCACGTCTTTGCGTTCCATCAGATTCAGCAGCgtggtg                                                                    >  1:754063/1‑296 (MQ=255)                                                                                                                                                                                                                              tgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAAT‑AGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCgtggtg                                                                    >  1:99895/1‑296 (MQ=255)                                                                                                                                                                                                                              tgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAAT‑AGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCCCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCCAACCGAACGTCTTTGCCTTCAATCAGCTTCAGCAGCgcggtg                                                                    >  2:526412/1‑296 (MQ=255)                                                                                                                                                                                                                                  cgATGCGACGCTGTCGCGTCTTATCGTGCCTACAAAT‑AGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGCCCGACACCCATCGTCCCCCGACCAAACCGACCGTCCTTGCCTTCCATCAGCTTCAGCACAGTGGTGCtcc                                                                >  2:820045/2‑294 (MQ=25)                                                                                                                                                                                                                                  cgATGCGACGCTGTCGCGTCTTATCGTGCCTACAAAT‑AGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGa                                                                <  1:654197/295‑1 (MQ=255)                                                                                                                                                                                                                                         aCGCTGTCGCGTCTTATCGTGCCTACAAAT‑AGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACg‑cgc                                                        >  1:54881/1‑296 (MQ=255)                                                                                                                                                                                                                                                  cgtcgttgcgggCCTACAAAC‑AGC‑CC‑GAACCGTAGGCAGGGTAAGGCGTTTACGCCTCATCCGGCAATTTGTGCAGGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGaaca                                               <  2:449904/285‑1 (MQ=255)                                                                                                                                                                                                                                                    tCTTATCGTGCCTACAAAT‑AGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTCAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGCTACACATCGTCGCCAGCCCAAACTGAACGCCTTTGCGTTACACCAGCTCCAGCACCGTGGCGCGGATCCG‑CGCACCgtccccac                                             >  2:771717/1‑289 (MQ=25)                                                                                                                                                                                                                                                              ggccgacgatgaggcc‑GAACCGTGGGCCGGTTAAGGCGTTTACGCCGCAGCCGGCAATTGGGGCAGGATGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGAACAACCCAGCGGa                                     <  2:725445/282‑1 (MQ=255)                                                                                                                                                                                                                                                                 aCAAAG‑AGT‑CC‑GGACCGTGGGCCGGATACGGGGTTTACGCCGCATCCGGAAATTGGGGCATGAGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGAACAACCCAGCGGATGAcc                                <  2:360762/296‑1 (MQ=255)                                                                                                                                                                                                                                                                          t‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGAACAACCCAGCGGATGACCCAGCGCGa                        <  1:841998/296‑1 (MQ=255)                                                                                                                                                                                                                                                                               gAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGAACAACCCAGCGGATGACCCAg                             <  2:390889/288‑1 (MQ=255)                                                                                                                                                                                                                                                                               gAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGAACAACCCAGCGGATGACCCAg                             >  1:390889/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                   cGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGG‑GCCGAATCGCCCCCCCGGACCAACCAAGCGGATGACCCacccggaccgccca                 >  1:254561/1‑283 (MQ=37)                                                                                                                                                                                                                                                                                    ggaGGGCGGAGAATGGGGTTACGCCGCACCCGTGAATTGGTGCATGTTGCCTGATGCGACGTGGGGGCGTCTTAGCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACGTCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCa                <  2:774318/294‑1 (MQ=255)                                                                                                                                                                                                                                                                                        gtcGGAGAAGGCGTTGACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGGCGGTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTATACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGTTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGt            <  1:406391/294‑1 (MQ=255)                                                                                                                                                                                                                                                                                          cGGATAGGGCATTTGCGCCGTATCGCGCAATTGGTGCATGGTGGCGTGTGAGAAGCTGGCGCGGCATATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACGCCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTg          <  2:118166/296‑1 (MQ=255)                                                                                                                                                                                                                                                                                            gATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAg        >  1:654205/1‑296 (MQ=255)                                                                                                                                                                                                                                                                                                  gCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGAt  <  1:646809/296‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                    CCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAAT‑AGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACG‑CGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGAT  >  NZ_LN832404/4023671‑4024252

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