Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,560,329	G→T	49.9%	pseudogene (125/606 nt)	glpR ←	pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,560,329	0	G→T	49.9%	46.1 ‑0.6 ‑0.4	48	pseudogene (125/606 nt)	glpR	pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  new base (11/12):  ref base (9/16):  total (20/28)
Fisher's exact test strand distribution p-value = 5.59e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.03e-01

TTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAGCCGTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAT  >  NC_000913/3560031‑3560613                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                              tttGGTGCGAACTTCGTGGGAATCGAAGCCCAGAAGGGAG‑‑GTAGCGATGGATGCCTCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGAGAAAATCGAGCGGCGCTTCGCCAATGAGCCCGCCATCGCGGCGTCGGAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGc                                                                                                                                                                                                                                                                                            <  1:414343/300‑1 (MQ=255)    ggTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAAGTCGCAAATTGCTGTGGTTGAGCAGTGCGTGCGCTACCGCTTACGGCGGGCTGCCGCTATCGAGAAACAGCGTAGCACCAtt                                                                                                                                                                                                                                                                                        >  2:217357/1‑301 (MQ=255)                 tgttaaccacactacggcatcGAGCCGTCTCTATCGATGCAGCTTATCCGCAGAATGCCGAAATCACGGCGGAACTGGGAGATAAAAGCGAGCGTCGCTGCGCAAATGCGCCCGTCATCGCGGCTGCGTAATTCGCCACCGGCGAGAAGGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTTTTGGTGACAATGCGCAAAGTGCTGTTATTGAGCAGTGCGTGCGATACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGc                                                                                                                                                                                                                                                                              <  2:327467/284‑1 (MQ=255)                 tGGCAAACGAACACCAGCAGCGAGCCGTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGc                                                                                                                                                                                                                                                                              <  1:379027/300‑1 (MQ=255)                          aaCTCCAGCAGCGAGCCGTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTTCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACtt                                                                                                                                                                                                                                                                    >  2:261851/1‑301 (MQ=255)                            atcaagctgcGAGCCGTCGCTATCGATGCCGCCAAGCCCCAGAAGGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGATGGTCGCTTCGCAAATGATCCCGCCATCGCGGATGCGTAATTCGCCACCGGCGAGAATGATGCGACAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTg                                                                                                                                                                                                                                                                  <  2:408373/297‑1 (MQ=255)                                                  acGAAGCCGCTGAGCGCCAGAAGGCCGAAAGAAAGGCGGAACTGGGGGATAAAATCGGGCGTTGCATCGCTGATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGGATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTttctt                                                                                                                                                                                                                                            <  2:130522/300‑1 (MQ=255)                                                               tCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTTACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGCTAAACAGCGTCTTGCCATTGCGGATTTGCTCCGCCACTTTGCGTGCGATGCGCTCTTTTTCTTCGGTCTGGGTgg                                                                                                                                                                                                                                >  2:73068/1‑300 (MQ=255)                                                                          ccGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTAACGCTTCCGGCGTGGTGCCGATATCGCTCAACAGCGTCTCGCCCTTGGGGATTTGCTCCGCCACTTTGCGGTCGATGCGATCTTTTTCTTCGGTCTGGGTGGCCTTGTGATCGt                                                                                                                                                                                                                    >  2:66576/1‑301 (MQ=255)                                                                                        gAACTGGGAGATAAAATCGAGCGTTGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTTCGTGCGCTACCGCTTCCGGAGTAGTTCCGCTATCGATAAACCGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCTGTCTTATTCTTCGGTCTTGGTGGCATTGCGATCGTGCCCCGGTGTGttt                                                                                                                                                                                                      >  2:171449/1‑300 (MQ=255)                                                                                        gAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCACCGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGCTTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCTTTGGGGCTTTGCTCCGCCACTTTCCGTGCGATGCGATCTTTTTCTTCGGTATGGGTGGCCTTGCGAGCGTGCCACGCCGTGTTa                                                                                                                                                                                                      >  2:109178/1‑301 (MQ=255)                                                                                           cTGGGGGATAAAATCGAGCGGCGCTTCGCCAAGGAGCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAcc                                                                                                                                                                                                   <  1:344193/301‑1 (MQ=255)                                                                                            tGGGAGATAAAATCGGGGGGCGCTGCGCCAATGATCCCGCCATCGGGGGTGCGTAATTCGCCACCGGCGAGAAGGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAcc                                                                                                                                                                                                   <  1:433326/300‑1 (MQ=255)                                                                                                                         aaggATCCCGCCAGCGCGGCGGCGTAAGTCGCCACGGGCGAGAATGATGCGAAAATCTGCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATGGCTGTGATTGGGCAGTGCGTGCGCTACCGCGTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCatga                                                                                                                                                                     <  1:20542/298‑1 (MQ=255)                                                                                                                         aaTGAACCCACCATAGTGGCTGCTTAATGCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTGAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCatga                                                                                                                                                                     <  2:223979/301‑1 (MQ=255)                                                                                                                                    cATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAATGCAGCGCCGCACCGCCATGCTGGCGCAGGAt                                                                                                                                                          >  2:358050/1‑301 (MQ=255)                                                                                                                                          ggtggcggaagtggacACCGGAGAGAAGGATGCGAAAATCTGCGTTGACAACAAACGGGTTGGAAACGGTGAGATGGTTGGTGACAATGCGCAAATTGCGGGGATGGAGCAGTGCGGGCGCTACCGCTTCCGGCGGGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGGTCCGCCACTTTGCGGGCGATGCGCTCTGTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGTTCAGGt                                                                                                                                                     <  1:399702/285‑1 (MQ=255)                                                                                                                                                 tAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGGTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAACGGGAGGGCATCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTGGcgc                                                                                                                                              >  2:40507/1‑298 (MQ=255)                                                                                                                                                           ccGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCCCCGCACCGCCCTGCTGGCGCAGGATCCGGTTTTGCTCCGCCAGCTc                                                                                                                                     >  1:348981/1‑299 (MQ=255)                                                                                                                                                             ggCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCCTTGGGGATTTGCTCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTCACCGCACTGGACGGTAGCGCCGCCCCGCCATGTTGGCGCCGGTACAGGTTTTGCTCCGCCCGCTcttcg                                                                                                                                 >  2:376692/1‑297 (MQ=255)                                                                                                                                                              gctaaaaTGATGCGAAAAGCGTCTTTGCACATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATCCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                 <  2:348990/295‑1 (MQ=255)                                                                                                                                                                 agaATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGGCGCAGGCTCAGGTTTTGCTCCGCCAGGTCCttacggt                                                                                                                             >  2:186204/1‑296 (MQ=255)                                                                                                                                                                                   tcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGCTGGCGCAGGCTCAGGTTTTGCTCCGCCGGCTCATTGAGGTCCCGGCGACTGGTCTggg                                                                                                           >  1:423470/1‑299 (MQ=255)                                                                                                                                                                                            atgagcGTGTTCGATACGTTGAGATGGTCGTTGACAATGCTCAAGTGGATGTGATTGCGCAGTGCCTGTGCCACCGCTGCTGGTTGTGTGCCGATGTCGATAACCAGCGTCGAGCCATGGGGGAGTTGTTCTGCCCAGTTGCGGGCGATGCGCTAGTTTTCGTCGGTCTGGGTGGCCTGGCGATCGGGCCACGGCGGGTTACCCGAACTGGAAGGCAGCGTCGCACCGCCATGATGTCGCAGGAGCAGGTTTGGCTCCGCAAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                  <  2:390340/296‑1 (MQ=255)                                                                                                                                                                                                 cGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGATGGCGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGGGGGCTGACGGAGaaa                                                                                             >  1:421519/1‑301 (MQ=255)                                                                                                                                                                                                 cGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGGTGGCGCAGGGTCCGGTTTTGCCCCGCCAGCTCCTTGAGGTCCCGGCGAATCGCCGGCGGGCTGCCGGGGaaa                                                                                             >  2:74461/1‑301 (MQ=255)                                                                                                                                                                                                               tggggtgggtgGTGACGATGGGCAAATGGCGGTGATTGGGCAGGGCGTGCGGTAGCGCTTCCGGAGTGGTGCCGATATCGGTAAATAGCGTCTGGCCATTGGGGATGTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctc                                                                               <  1:73068/292‑1 (MQ=255)                                                                                                                                                                                                               tgggaTGGTTGGGGCAAAGGCGCAAAGTGGTGGGATTGAGCAGTGCGTGGTCTACCGGTTCCGGCGTGGGGCCGAGAGCGAAAAACAGCTGCGAGCCATTGGGGATTGGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTGCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctc                                                                               <  1:217352/298‑1 (MQ=255)                                                                                                                                                                                                               ggtggttgagggttcaaaTGCGGAAATTGTTGTGGTGGAGCAGGGGGGTGGCTACCGTTGATGGGGTGGTGCCTATATCGATAAGAAGGGTCGTGCCAGTGGGGATTTGCTCCGAGACGTGGCGGTGTACGCGCTCTTGTCCTTAGGTCTGGGTGGCCTCGCGATCGCGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGct                                                                                <  1:433727/284‑1 (MQ=255)                                                                                                                                                                                                                            gACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCTTTGCGGTCCCGGCGAATAGTCTGCGGCCTGACGGAGGAATGATCTACCCGCCCTTCGGTCATGACa                                                                  >  2:238145/1‑301 (MQ=255)                                                                                                                                                                                                                                 tgcgtaattggatgtgATTGGGGAGTGCGTGAGCTACCGGTTCTGGCGGGGTGCCGAGAGCGATAAACAGGGTCTAGCCGGTGGGGGTTGGCTCCGCCACTTGGCGGGCGATGCGCTCTGTTTCTGCGGTCTGGGTTGCCTTGCGATCGTGCCACGGCGTGTTAAGCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                                             <  1:241798/289‑1 (MQ=255)                                                                                                                                                                                                                                 tgagcaatgtgttgtgATTGGCCAGTGCGGGCGCGACCGCTCGCGGCGGGGGTCCGATTGTGATAAACAGCGTCGAGCCATGGGTGAGTCGCTCCTCCACTGTGCGTGCGATGCGCTCGTGTTCTTCGGGCTGGGTGTCCGTGCGATGGTGCCACTGCGTGTTAATCGAACTGGAAGGCAGCGCCGCACCGCCATCAGGGCGCAGGCTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                                             <  1:187200/289‑1 (MQ=255)                                                                                                                                                                                                                                   ggcAAATTGCTGTGATTGGGAAGGGCGGGCGGGCGGGGTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTGTGCGGGCGATGCGGTCTTTTTCTTCGGTCTGGGTGGCCTTGCGTTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                            <  1:186203/299‑1 (MQ=255)                                                                                                                                                                                                                                   cgtcatgtttggtgATGGAGCAGTGCGTGTGCTACCGCTGCGGGGGTTGTGCCGATCGCGATAAACAGCGTCGAGCCATGGGGGCTTTCCTCCGCAACTTTGCTGGCGATGCGCTCTTTTTCTGCGGTCTGGGGGGCCTTCCGATCCTGCCACGGCGTGTTAACCGAACTTGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                            <  1:171448/289‑1 (MQ=255)                                                                                                                                                                                                                                   cgtaaaTTGGCGTGCTTGAGAAGTGCGTGGGCAACGCGTTCCGGGTTGGTGACGATACCGGTAACCAGCGGCTTGCCATTGGGGATTTGCTCCGCGACTGGGCGGGCTATGCGCTCTTTTTCTTGGGTCTGGGTGGCCTTGCGAGCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCACGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                            <  1:66576/297‑1 (MQ=255)                                                                                                                                                                                                                                   cggtaaTTGCTGTGATGGGGCAGTGTGTGCCGTACCGCTTCCGGCGTGGGGCCGATATCGATAACCAGCGTCGAGCCATTGGGGATTTGCTACGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                            <  1:304583/296‑1 (MQ=255)                                                                                                                                                                                                                                    gcAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGACATGCTCTAACAGCTCTTCGGTAc                                                                       >  2:404175/1‑288 (MQ=255)                                                                                                                                                                                                                                    gcAAATTGCTGTGATTGAGCAGGGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                                       <  1:404166/288‑1 (MQ=255)                                                                                                                                                                                                                                              tgggttGAGCAGGGGGGGGGGTACGGGTTCCGGCGGGGGGCCGATATCGATAAACAGCGTCGGGCAATTGGGGAGTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGAACTGACATAACCCTGCTGTTTAAcc                                                <  1:358041/297‑1 (MQ=255)                                                                                                                                                                                                                                              tggcttggtcAGTGCGTTCGGGTCCTCTGCCGGCTTGGTGCCTCTATCGATAACCACGCTCTGGCCATTGGGTATTTGCTCCGCCACTTTGCGGGCTATGCGGTCGTTTTCTTCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGAACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                                <  1:162958/292‑1 (MQ=255)                                                                                                                                                                                                                                                 gtttggagggcggggcggtcccGGTTCCGGCGTGGTGCCGAGTTCGATAAACGGCGTCTGGCAATTGGGGATTTGCGCCGCCACTTTGCGGGGGATGCGCTCTGTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGGGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTGGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                                              <  1:376683/280‑1 (MQ=255)                                                                                                                                                                                                                                                         gtgcgtgcGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGCTTTGCTCCGCCACTTTGCGGGCGCTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACCGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATCATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGCATAGTCCGGGGGCTGACGGCGAAATTCTCTACCAGCTCTTCGGTACTTCCATAACGCTGCTGTTCAACCCGTTCGaaaaa                                     >  1:288276/1‑300 (MQ=255)                                                                                                                                                                                                                                                            cgtgcgAGACCGCCTGCGGGGTGGTGCTGAGATGGCTAAAGACCGTTTAGCCATGGGGGATTGGCGCCGCACGTTTGCGGGCGAGCCGCTCTTTTTCTTCGGTGTGGGTGGCCTTGCGACCGTGCCACGGCGTGTTAACAGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCATGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGGCTGCTGGCTGACGGAGAAATGCTCTACCAGCTCGTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                   <  2:10348/300‑1 (MQ=255)                                                                                                                                                                                                                                                             gtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGATGGCGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATCGTCTGCGTGCTGACGTAGAAATGCTCTCCCAGCTCTTCGGTGCTGACATAACCCTGCTGGTTAACAATTTCGAAAATGCCg                                 >  1:134751/1‑301 (MQ=255)                                                                                                                                                                                                                                                                          ccccgCGTGGTGCCGATATCGAGAAAGAGCGTCGAGCGATGGGGGATTTGCTCCGCCATGTGGCGGGCGATGCGCTCTGTTTCTGCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                    <  1:142919/297‑1 (MQ=255)                                                                                                                                                                                                                                                                            cggggtggttccgcttggaaTAATCAGCGCGTAGTCATTGGGGTTTTGGTCGGCCACTGTGCGGTCGAGGGGCTGTTGTTGTTAGGTCTTGGTGTCCGTGCGAGCGGGCCACGGTGTGGTAACCGAACTGGAAGGCAGCGCCGCAACTCCATGATTGCTCAGGATCAGGTTTTGCCCCGCCAGCTCATCGAGGTCGCTGCGAATAGTCGGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtt                  <  2:423479/282‑1 (MQ=255)                                                                                                                                                                                                                                                                                            tCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGCATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGTCCTGACATCACCCTGCGGTTTAACCAGTTCGATCATACCGTTGTGCCTTTGTGTTTTTTTCATTTATAAAt  >  2:308052/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                            tCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGGCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCGTCGGTACTGACATAACCATGCTGGTTCACCCTGTGGATAATACCGTTGTGCCGTTGAGTTTGTTTCTTTtatacaa  >  2:28934/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                              TTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAGCCGTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAT  >  NC_000913/3560031‑3560613

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.