Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_LN832404	2,771,251	A→G	81.8%	S544G (AGT→GGT)	EK12AG100_RS13785 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_LN832404	2,771,251	0	A	G	81.8%	68.8 / 10.7	33	S544G (AGT→GGT)	EK12AG100_RS13785	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (17/9);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

ATGATGACATAAGCAATGGAAAATTATCTTCACTGAGAGCAACAGGTATCTCCCAGTTTAGTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATGGAA  >  NZ_LN832404/2770961‑2771537                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                               atgatgACATAAGCAATGGAAAATTATCTTCACTGAGAGCAACAGGTATCTCCCAGTTTAGTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCt                                                                                                                                                                                                                                                                                            <  1:673097/295‑1 (MQ=255)       aCATAAGCAATGGAAAATTATCTTCACTGAGAGCAACAGGTATCTCCCAGTTTAGTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAACGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGc                                                                                                                                                                                                                                                                                      >  2:510996/1‑295 (MQ=255)                          aTCTTCACTGAGAGCAACAGGTATCTCCCAGTTTAGTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAattcatt                                                                                                                                                                                                                                                                   >  1:216334/1‑295 (MQ=255)                                               tATCTCCCAGTTTAGTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAatgat                                                                                                                                                                                                                                              <  1:349419/295‑1 (MQ=255)                                                                       gAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTaa                                                                                                                                                                                                                      >  1:707105/1‑295 (MQ=255)                                                                                        ggCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGACTGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAAtt                                                                                                                                                                                                     >  1:539252/1‑295 (MQ=255)                                                                                                tCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCtt                                                                                                                                                                                             <  1:293642/295‑1 (MQ=255)                                                                                                                cTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACtatttatt                                                                                                                                                                             <  1:32748/295‑1 (MQ=255)                                                                                                                              tATGAAATAGATAATAAAAGAAAAAAAAACAAAATATAAGGTACAGCAGTTATTGTTTTTGTGTATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGt                                                                                                                                                               <  2:287457/295‑1 (MQ=255)                                                                                                                                      aGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAATTATCCTAAGTTAATCCTg                                                                                                                                                       >  2:164844/1‑295 (MQ=255)                                                                                                                                       gTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGCGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGt                                                                                                                                                      >  1:603623/1‑295 (MQ=255)                                                                                                                                              aaaGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTaaa                                                                                                                                               >  1:125332/1‑295 (MQ=255)                                                                                                                                                       aaaTCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCACTTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAACTATCCTAAGTTAATCCTGTTTTTAAAGCAGTCTTg                                                                                                                                      >  2:709667/1‑295 (MQ=255)                                                                                                                                                            aaaaTATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGt                                                                                                                                 <  1:155611/295‑1 (MQ=255)                                                                                                                                                                                tCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTg                                                                                                             <  1:510996/295‑1 (MQ=255)                                                                                                                                                                                          gttgATAAAAATATAATGAACCTGCGTAAAATAACTTATTTGCAATCTAAGCTAGTAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGGCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTg                                                                                                   <  2:539252/295‑1 (MQ=255)                                                                                                                                                                                                aaCAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGttt                                                                                             <  2:137660/295‑1 (MQ=255)                                                                                                                                                                                                 aCAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCTGACACTCATACCTGCTTTAAAGTGGCTtttt                                                                                            >  2:216293/1‑295 (MQ=255)                                                                                                                                                                                                     gATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAACTTACGCCTGTTTTTAAATCAGGCTTGAGCGTTAAGCTGACACTCATACCTTCTTTAAAGTGCCATTTTTTTc                                                                                        >  2:671939/1‑295 (MQ=255)                                                                                                                                                                                                       tAATGAACCTGCTTAAAAGAACTTATTTTCAAGTGAAGCTAGAAATTTTTGGTTCTTCTATTATTGATATATTAAAGTATATATTTGGTTTATGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACGTTAAGAATGATGTGAATTAAAATTTATTCGTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGt                                                                                      <  2:532207/295‑1 (MQ=255)                                                                                                                                                                                                              ccTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCAt                                                                               >  2:546312/1‑295 (MQ=255)                                                                                                                                                                                                                      aaaTAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCTTTACCAATACCTatat                                                                       >  2:281778/1‑295 (MQ=255)                                                                                                                                                                                                                        aTAACTTATTTGCAATCTAAGCTAGTAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGCCTATTTATTACAATATCCTACGTTAAGCCTGTTTTAAAAGCTGTCTTGATCGTGAAGCGGCCACTACTACCTGCTTTAAAGTGGCGTATTTGTCGTTCCTCATAACTATAttt                                                                     >  2:401183/1‑295 (MQ=255)                                                                                                                                                                                                                        aTAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGATCGTGAAGCGGACACTCATACCTTCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTTTAttt                                                                     >  2:286664/1‑295 (MQ=255)                                                                                                                                                                                                                                tttGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCtgttgt                                                             <  2:216334/295‑1 (MQ=255)                                                                                                                                                                                                                                 ttGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGGCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTgttgtt                                                            <  2:212986/295‑1 (MQ=255)                                                                                                                                                                                                                                     aaTCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAAc                                                        >  1:502391/1‑295 (MQ=255)                                                                                                                                                                                                                                      aTCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACt                                                       >  2:565405/1‑295 (MQ=255)                                                                                                                                                                                                                                                               tctATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAATATATCCTAAGTTAAGCCTGTTTTTAAAGAAGGCTTGAGAGTGAAGAGGACAACCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTAACAATCCCTATATTTTCTGTTTTTTCACTGCCAGAACTGAGCGGTTGAAAAgac                              >  2:526834/1‑293 (MQ=255)                                                                                                                                                                                                                                                                     atGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAat                        >  1:391672/1‑295 (MQ=255)                                                                                                                                                                                                                                                                               tAAAGTATATATTTGGTTTAGGTCTTCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGcc              <  2:125332/295‑1 (MQ=255)                                                                                                                                                                                                                                                                                 aaGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAg            >  1:70775/1‑295 (MQ=255)                                                                                                                                                                                                                                                                                           ttGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATGGaa  >  1:508707/1‑295 (MQ=255)                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                               ATGATGACATAAGCAATGGAAAATTATCTTCACTGAGAGCAACAGGTATCTCCCAGTTTAGTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATGGAA  >  NZ_LN832404/2770961‑2771537

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