Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	480,295	2 bp→AA	coding (13‑14/219 nt)	hha ←	hemolysin expression‑modulating protein Hha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	480,295	0	G	A	100.0%	134.3 / NA	40	P5F (CCT→TTT)	hha	hemolysin expression‑modulating protein Hha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/16);  total (24/16)
*	NC_000913	480,296	0	G	A	97.5%	125.7 / NA	40	P5F (CCT→TTT)	hha	hemolysin expression‑modulating protein Hha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/16);  minor base C (1/0);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTTGTAATTAAGTGCG  >  NC_000913/480065‑480535                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                                tataGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCt                                                                                                                                                                                                                               <  1:283176/250‑1 (MQ=255)  ataGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTAc                                                                                                                                                                                                                             >  1:234388/1‑251 (MQ=255)    agtcccGAAAGCGAATTATTTTTTAACGAAAAAATTGCCATACTGAGGAAGGGATCTGTTCGTCTAGTTTATT‑CATGGCCAATTCGGCGGGGCGGGGCTCTGCGGCTGATAAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTAccg                                                                                                                                                                                                                           <  2:442279/247‑2 (MQ=255)                gAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGtaat                                                                                                                                                                                                                >  1:625400/1‑250 (MQ=255)                    tatttattAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGTACTGCGGCTGAGTTAAATACCGCCAGTTCATTATCTGGTAATTTGGATTTATTTTTCTCG‑ATAACACGCTCCCGCGTGTCAATTGTCTTGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTTGGTCATACTTCTACCTATGGTTGTAATAAttt                                                                                                                                                                                                           >  2:193518/1‑250 (MQ=255)                    tatttattAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAAttt                                                                                                                                                                                                           >  1:692981/1‑250 (MQ=255)                       ttattaGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTa                                                                                                                                                                                                         >  1:205959/1‑250 (MQ=255)                         attaGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCCTTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAAc                                                                                                                                                                                                       >  2:139233/1‑250 (MQ=255)                         attaGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGGTAATTCGTATTTATTTTTCTCG‑ATAACACGCCCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTTTAATAATTTTAAc                                                                                                                                                                                                       >  2:63316/1‑250 (MQ=255)                                        tCCATAGTGGGGAGGGGATCTTATCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTTAAATACCGCCCGTTCCTTATCTTGTCATTCGTGTTTATTTTTCTCGCATC‑CCCGCTCCAGCGTGGCCATTTTCTGGCAAAGGCGTAAACGCATTAAATAATCGGGTTTCGTTACAAATTTTTCGGACATACTTTTTCCTATGGTTGTAATTATTCTAACAAAGTAAACTCGCa                                                                                                                                                                                         >  2:534411/1‑249 (MQ=255)                                         ccATACTGAGGAAGGGATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCCGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGCTAAACTCGCAgg                                                                                                                                                                                       >  1:514345/1‑250 (MQ=255)                                                 gggAAGGTATCTTGTCGTAAAGTTTATT‑CATGTTCAATTCGGCGAGGCGGTGATCTGCGGCTGATTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCttt                                                                                                                                                                              <  2:682887/250‑1 (MQ=255)                                                  ggAAGGGATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCCTTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCtttt                                                                                                                                                                             >  1:394644/1‑250 (MQ=255)                                                  ggAAGGGATCTTGTCGTA‑AGTTTTTTTCAGGGTAAATTCGGCAGGGCTGTGTTCTGCGGTTGATTAAAATTCCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTtc                                                                                                                                                                             <  2:289659/251‑1 (MQ=255)                                                        gATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCCGTTCATTATCTGGTAATTCGTATTTTTTTTTCTCG‑ATAAAACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATTGTTGTAATTATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGcc                                                                                                                                                                       >  2:211399/1‑250 (MQ=255)                                                                                             gAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCCGTTCATTATCTGATAATTCGCATTTATTTTTCTCG‑CTAACACGCTCCAGCGTGTCAATTGTCTGGCAAAGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTAGGAAATACCTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAAAAACGGAATAGTCGATTACCTGACTTc                                                                                                                                   >  1:744569/1‑251 (MQ=255)                                                                                                 gggTGATCTGCGGTTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGcc                                                                                                                               <  2:212550/250‑1 (MQ=255)                                                                                                                          ccGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAAt                                                                                                      >  1:118003/1‑251 (MQ=255)                                                                                                                                        tCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTg                                                                                         >  2:709616/1‑250 (MQ=255)                                                                                                                                            aTAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGaaca                                                                                     >  1:25755/1‑250 (MQ=255)                                                                                                                                            aTAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTAAATTGTCTGGCAACGTCGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGaaca                                                                                     >  2:375342/1‑250 (MQ=255)                                                                                                                                                     atttatttTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAgg                                                                            <  1:709616/250‑1 (MQ=255)                                                                                                                                                      tttatttTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATATTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGgtg                                                                          >  2:59636/1‑251 (MQ=255)                                                                                                                                                             ttCTCG‑CTAACACGCTCCAGCGTGTCGATTGTCTGGCAACGACGTCAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATcc                                                                    >  2:332381/1‑250 (MQ=255)                                                                                                                                                               ctcG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGCCGTAAACGCATTAAATAATCGGTTTTCGTTAAAACTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATcc                                                                    >  2:807182/1‑248 (MQ=255)                                                                                                                                                                     tAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATtc                                                            >  1:322218/1‑251 (MQ=255)                                                                                                                                                                          cGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATtcttcg                                                        >  1:592040/1‑250 (MQ=255)                                                                                                                                                                                aGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCt                                                 >  2:227529/1‑251 (MQ=255)                                                                                                                                                                                  cGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCa                                               <  1:63316/251‑1 (MQ=255)                                                                                                                                                                                       cAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGa                                           <  2:205959/250‑1 (MQ=255)                                                                                                                                                                                           tGTCTGGCAACGACGTAAACGCATTACATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGCTAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTCCTGAACAACATAAAGGTGTCATCCAGATATTCTTCAATCTGCTCAATGAGCtt                                       >  2:745076/1‑250 (MQ=255)                                                                                                                                                                                                ggCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACGCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCAGTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTg                                   <  1:117737/249‑1 (MQ=255)                                                                                                                                                                                                          gaaaCGCAATAAATAATCGGTTTTCGTAAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCCCGTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCAATTCGGAAGATCCTGCAGATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGCCGAGCTGCTCAATGAGCTTATTGGCTTCATTATAc                       <  2:744375/250‑1 (MQ=255)                                                                                                                                                                                                             aCGCATTAAATAAGCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTa                    <  1:293097/251‑1 (MQ=255)                                                                                                                                                                                                              cGCATTAAATAATCGGTTTGTGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTaa                   <  2:234388/251‑1 (MQ=255)                                                                                                                                                                                                               gCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAAt                  <  2:394644/251‑1 (MQ=255)                                                                                                                                                                                                                cATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATTCCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAAtt                 <  2:25755/251‑1 (MQ=255)                                                                                                                                                                                                                    aaaTAATCGGTTTTCGTTAAAAAATTTTCGGACATTCTTCTACCTATGGTTGTAATAATTCTAACAAGAAAAACCCGTAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTTg              <  2:426288/250‑1 (MQ=255)                                                                                                                                                                                                                            ggTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTTGTAATTAAGt     <  1:250573/251‑1 (MQ=255)                                                                                                                                                                                                                               tttCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTTGTAATTAAGTGCg  >  1:318002/1‑251 (MQ=255)                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                                TATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATT‑CATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCG‑ATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTTGTAATTAAGTGCG  >  NC_000913/480065‑480535

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