Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,728,611	C→T	G254E (GGG→GAG)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,611	0	C	T	90.8%	187.2 / ‑4.3	66	G254E (GGG→GAG)	xylB	xylulokinase
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base T (42/17);  total (43/23)
Fisher's exact test for biased strand distribution p-value = 1.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAAGATAAGGCAGAAACCAAACTGGCTCGGCACTTTCATCAGCCTGTTGAGCTGCAGCGATTAAAGCTGGGACATTGCTCAGGCCGGTTAATTTCGCGGCCCAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTG  >  NC_000913/3728370‑3728860                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                           ggAAAGATAAGGCAGAAACCAAACTGGCTCGGCACTTTCATCAGCCTGTTGAGCTGCAGCGATTAAAGCTGGGACATTGCTCAGGCCGGTTAATTTCGCGGCCCAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGAt                                                                                                                                                                                                                                                    >  1:262178/1‑250 (MQ=255)      gATAAGGCAGAAACCAAACTGGCTCGGCACTTTCATCAGCCTGTTGAGCTGCAGCGATTAAAGCTGGGACATTGCTCAGGCCGGTTAATTTCGCGGCCCAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACAtt                                                                                                                                                                                                                                              >  1:164346/1‑251 (MQ=255)                         tGGCTAGGCACTTTCATCAGCCTGTTGAGCTGCAGCGATTAAAGCTGGGACATTGCTCAGGCCGGTTAATTTCGCGGCCCAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAAccc                                                                                                                                                                                                                          >  2:452339/1‑250 (MQ=255)                                  aCTTTCATCAGCCTGTTGAGCTGCAGCGATTAAAGCTGGGACATTGCTCAGGCCGGTTAATTTCGCGGCCCAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCacac                                                                                                                                                                                                                 >  1:180003/1‑251 (MQ=255)                                     ttCATCAGCCTGTTGAGCTGCAGCGATTAAAGCTGGGACATTGCTCAGGCCGGTTAATTTCGCGGCCCAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCa                                                                                                                                                                                                               >  1:225722/1‑250 (MQ=255)                                           aGCCTGTTGAGCTGCAGCGATTAAAGCTGGGACATTGCTCAGGCCGGTTAATTTCGCGGCCCAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCAc                                                                                                                                                                                                        <  1:98397/251‑1 (MQ=255)                                                gTTGAGCTGCAGCGATTAAAGCTGGGACATTGCTCCGGCCGATTAATTTAACGACCCCATCCAGACCCGGCGCTGCACTCAGCCTCAAAGACATTAAATGCCACCGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCt                                                                                                                                                                                                   >  1:341171/1‑251 (MQ=255)                                                  tGAGCTGCAGCGATTAAAGCGGTGACATGGCTCAGGCCGGTGAATTTCGCGGCCCAATCCAGACACGACGCTGCCCTTAGCATCACAGACATTATATGCCAACGTTGCGGTAGCGCATGGCCAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTGCGCTGACAGCAAAATAGACCCCCGACTTCTCCAGCGATAACATTG‑CCTGATTAGCATCACCCATTCCCACACCATCTGCTCCAGCTg                                                                                                                                                                                                  <  2:37335/250‑1 (MQ=255)                                                   gAGCTGCAGCGATTAAAGCTGGGACATTGCTCAGGCCGGTTAATTTCGCGGCCCAATCCAGACACGACGCTGCACTCAGCATCACAGACATTCAATGCCAACGTTGCGGTAGCGCATGGCAACAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGc                                                                                                                                                                                                 >  2:159555/1‑250 (MQ=255)                                                             gATTAAAGCTGGGACATTGCTCAGGCCGGTTAATTTCGCGGCCCAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGATGCGGTAGCGCATGGCACAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCCGCGATAACATTG‑CCTGCTTAGCATCAACCATTCCCACACCAACTGCCCCAGATGCATTATCGCCAc                                                                                                                                                                                      >  2:394566/1‑251 (MQ=255)                                                                      tGGGACATTGCTCAGGCCGGTTAATTTCGCGGCCCAATCCCGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGCACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCACCTGCACCAGCTGCATTGTCGCCACCGCCt                                                                                                                                                                                 >  2:38660/1‑247 (MQ=255)                                                                      tGGGACATTGCTCAGGCCGGTTAATTTCGCGGCCCAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCt                                                                                                                                                                                 <  1:38660/247‑1 (MQ=255)                                                                         gACATTGCTCAGGCCGGTTAATTTCGCGGCCCACTCCAGCCACGACGCTGAACTCAGAATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAACTAGACCCCCGACGCCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCCTTCCCACACCAACTGCCCCAGCTGCATTGACGCCACCGCCGGCGACa                                                                                                                                                                           >  2:226207/1‑250 (MQ=255)                                                                                cTCAGGACGGTTAAATTCGCTGAACCATCCAGACACCACCAGGCACTCAGCATCCCAGAAATTAAATGCCACCGTTACGGTAGCCCATGGCCAAAGATATGTAAGGCGCTTTCTGGCTTAATTAAGAACCCTTCGCTGACAGCAAAATAGAACCCCGACGCCTACAACGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTCCACAAGCTGCATTGTCCCCACCGCCTGCGCCAACTGCCAc                                                                                                                                                                   >  2:469583/1‑251 (MQ=255)                                                                                            aaTTTCGCGGCCCAATCCAGACACGACGCTGCACTCCGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATAcc                                                                                                                                                       >  2:15014/1‑251 (MQ=255)                                                                                                ccgcgGCCCAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCGCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCAc                                                                                                                                                   <  1:458037/250‑1 (MQ=255)                                                                                                     ccccAATCCAGACACGACGCGGCCCTCAGAATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAAAAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCACCAACCATCCCCCCACCAACTGCCCCAGCTGCTTTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCttt                                                                                                                                              <  2:164346/250‑1 (MQ=255)                                                                                                       ccAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCg                                                                                                                                            >  1:32600/1‑251 (MQ=255)                                                                                                        cAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTACGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCGGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGc                                                                                                                                           >  2:243058/1‑251 (MQ=255)                                                                                                        cAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGc                                                                                                                                           >  2:243068/1‑251 (MQ=255)                                                                                                          aTCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCa                                                                                                                                          >  1:495819/1‑250 (MQ=255)                                                                                                             cAGACACGCCGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCAGGGCAAAAGCTATTTACGGCGCTTTCTGGTTTGCTTAAGAACCCTCCGCTGACAGCAAAATAGCCCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACtt                                                                                                                                      <  1:333244/251‑1 (MQ=255)                                                                                                                      cgCTGCAATCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGACCCCTTCGCTGAAAGCAAAATAGACCCCCGACGTCTCCAGCGATAACCTTG‑CCTGATTAGCATCAACCATTCCCACACCAACAGCACCAGCTGCATTGACGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCAAGCTTTCGCTCCTTCAGGTAAc                                                                                                                              >  2:303751/1‑250 (MQ=255)                                                                                                                       gCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACa                                                                                                                             >  1:422334/1‑250 (MQ=255)                                                                                                                          gCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGACCCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACGGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACaaa                                                                                                                           <  1:412790/249‑1 (MQ=255)                                                                                                                                 ccATCACAGAAATTAAAGGCCACCGTTGCGGGAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTGGCTTAAGAACCCTTCGCTGACAGCAAAACAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTATCATCAACCATTCCCACCCCAACTGCCCCAGCGGCATGGTCGCGACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTa                                                                                                                  <  2:512597/250‑1 (MQ=255)                                                                                                                                     cacaGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTg                                                                                                                                                                                >  1:2427/1‑185 (MQ=255)                                                                                                                                     cacaGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTg                                                                                                                                                                                <  2:2427/185‑1 (MQ=255)                                                                                                                                      acaGACATTAAATGCCAACGTTGCGGTAGCGCATAGCAAAAGCTATATACAGCGCTTTCTGGCTTGCTTAAGAACCCTACGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACCTTG‑CCTGACTAGCATCAACCATTCACACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGACACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGAAACAAAGCAACAGTAATTTc                                                                                                             >  1:252607/1‑251 (MQ=255)                                                                                                                                        aGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGc                                                                                                           >  2:115549/1‑251 (MQ=255)                                                                                                                                        aGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCACCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGc                                                                                                           >  1:60672/1‑251 (MQ=255)                                                                                                                                                   gCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCAc                                                                                                                                                                                                        <  1:479352/147‑1 (MQ=255)                                                                                                                                                   gCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCAc                                                                                                                                                                                                        >  2:479352/1‑147 (MQ=255)                                                                                                                                                        cGTTGCGGGAGCGCAGGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCGGACAGCAAAATAGCCCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTTCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATTCCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGc                                                                                           <  2:492949/251‑1 (MQ=255)                                                                                                                                                            gCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGc                                                                                       >  2:99792/1‑251 (MQ=255)                                                                                                                                                                 aGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTg                                                                                  >  2:81366/1‑251 (MQ=255)                                                                                                                                                                    gcATGGCAAAAGCTATGTACGGCGCTTTATGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGCCCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTCAGACAACTGGCACCGTCGCCATACCCCACGCt                                                                                                                                                >  2:283217/1‑186 (MQ=255)                                                                                                                                                                      aTGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCCGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCAc                                                                             >  2:276559/1‑251 (MQ=255)                                                                                                                                                                           aaaaGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACgaggt                                                                        >  2:276081/1‑249 (MQ=255)                                                                                                                                                                             aaGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATaa                                                                      <  1:99792/251‑1 (MQ=255)                                                                                                                                                                                 tATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCg                                                                  >  1:286289/1‑251 (MQ=255)                                                                                                                                                                                  aTGTACGCCGCTTTCTGGCTTGCTTAAGAACCCTCCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCg                                                                  <  1:501317/250‑1 (MQ=255)                                                                                                                                                                                        ggCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCt                                                           <  1:3465/251‑1 (MQ=255)                                                                                                                                                                                           gctttagggcctgcttAAGAACCCTACGCTGTCAGCAAAATAGCCCCCCGACGTCCCCAGCGATAACATTGACCT‑ATTAGCATCAACAATTCCCACACAAACTGCCCCAGCTGCATTGTCGCCACCGCCTGCGACAACTCGCACCGACGCCATCACCCACGCTTTCGCAACTTCCGGTACCAAAGAACCAGAAATTTCGCTGAATTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTgc                                                         >  2:65886/12‑250 (MQ=255)                                                                                                                                                                                                cgggctggctTAAGAACCCTTCGCGGACAGCAAAATAGCCCCCCGACGTCTCCAGCGATAACATTG‑CCTGATGAGCATCAACCATTCCCACACCAACTGCCCCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCAt                                                    <  2:262178/248‑1 (MQ=255)                                                                                                                                                                                                  ggCTTGCTTAACAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACCTTG‑CCTGATTCGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGAAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGAc                                                 >  2:457979/1‑251 (MQ=255)                                                                                                                                                                                                   gcttgcttAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACGGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACg                                                <  1:81366/251‑1 (MQ=255)                                                                                                                                                                                                    cttgcttAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCCTTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGt                                               >  2:51842/1‑251 (MQ=255)                                                                                                                                                                                                      ggattAAGAACCCTTCCCTGAAAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTG‑CCTGATGAGCACCACCCCTTCCCACACCATCTGCCCCAGCTGCATGGCCGCCCCCGCCTTCGCCAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTGCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTc                                              <  2:254062/247‑1 (MQ=255)                                                                                                                                                                                                        cttAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACt                                           >  2:164644/1‑251 (MQ=255)                                                                                                                                                                                                            aGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTcc                                         >  2:91339/1‑249 (MQ=255)                                                                                                                                                                                                            aGAACCATTAGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCAACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCAACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTcc                                         >  2:171809/1‑249 (MQ=255)                                                                                                                                                                                                            aGAACACTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCCCACCCACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTACCAAAGCACCAGTAATTTCGCTGCCTTCGTATACTGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTcc                                         >  2:443172/1‑249 (MQ=255)                                                                                                                                                                                                              accccTTCGCTGACAGCAAAAAAGACCCCCGCCGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTc                                     <  1:452339/249‑1 (MQ=255)                                                                                                                                                                                                                 ccTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACg                                  >  2:26896/1‑251 (MQ=255)                                                                                                                                                                                                                  cTTCGCTGACAGCAAAATAGACCCCCGACGTCTCCCGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACCCCACCTGCACCAGCTGCATTGTCGCCACCGCATGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCCACTTCAGGTAACAAAGCACCAGCAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATACGTCGCAAGCCTGCAGCATGACGTCCCTCCAGTCACGc                                 >  2:28158/1‑251 (MQ=255)                                                                                                                                                                                                                            aGCAAAATAGACCCCCTACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACAt                       <  1:510349/251‑1 (MQ=255)                                                                                                                                                                                                                              cAAAATAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATcc                     >  1:85424/1‑251 (MQ=255)                                                                                                                                                                                                                                aaaTAGAACCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATcc                     >  1:402191/1‑249 (MQ=255)                                                                                                                                                                                                                                  aTAGACCCCCGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGcc                 >  2:205518/1‑251 (MQ=255)                                                                                                                                                                                                                                      ccccccGACGTCCCAAGCGATAACATGG‑CCGTATGAGCATCACCCATTCCCACACCAGCTGCCCCAGCTGCATTGTCCCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGTTGCCTCCGTATAAGGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCATCCACAt             <  1:394566/250‑1 (MQ=255)                                                                                                                                                                                                                                       cccccGTCGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATg            >  1:150936/1‑251 (MQ=255)                                                                                                                                                                                                                                          ccGACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCCCCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTg         <  1:243068/251‑1 (MQ=255)                                                                                                                                                                                                                                          ccGACGTCCCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACCCCAACTGCCCCAGCTGCATTGTCGCCCCCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATggcg         <  1:243058/251‑3 (MQ=255)                                                                                                                                                                                                                                            gACGTCTCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGcc       >  2:25153/1‑251 (MQ=255)                                                                                                                                                                                                                                              cGTCCCCAGCGATAACATTG‑CCTGATTAGCATCACCCATTCCCACACCAACTGCCCCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAg     <  1:427477/251‑1 (MQ=255)                                                                                                                                                                                                                                                 ccccAGCGATAACATTG‑CCTGATTAGCATCACCCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTg  <  1:36083/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                           GGAAAGATAAGGCAGAAACCAAACTGGCTCGGCACTTTCATCAGCCTGTTGAGCTGCAGCGATTAAAGCTGGGACATTGCTCAGGCCGGTTAATTTCGCGGCCCAATCCAGACACGACGCTGCACTCAGCATCACAGACATTAAATGCCAACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTG‑CCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTG  >  NC_000913/3728370‑3728860

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