Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,088,390	T→C	I36T (ATT→ACT)	galP →	galactose:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,088,390	0	T	C	100.0%	187.2 / NA	52	I36T (ATT→ACT)	galP	galactose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/22);  total (30/22)

ACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATTTTGTTTGTTGCCGGTTCGCTGTTCTCTGCGGCTGCGCCAAACGTTGAAGTACTGATTCTTTCCCGCGTTCTACTGG  >  NC_000913/3088170‑3088629                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                aCTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATg                                                                                                                                                                                                                    >  1:597783/1‑250 (MQ=255)      gATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCCGTGTAACCGCTGGCGCACTGCCGTTTTTTGCAGATGAATTcc                                                                                                                                                                                                              >  2:479423/1‑251 (MQ=255)       aTAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTTGCGCACTGCCGTTTATTGCCGATGAATTcc                                                                                                                                                                                                              >  2:634099/1‑250 (MQ=255)               tAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGCTATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACtt                                                                                                                                                                                                     >  1:600447/1‑251 (MQ=255)                     aaGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGTCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGTCCTGGATATCGGTGTAACTGCTTGCGCTCTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGcac                                                                                                                                                                                                <  1:496663/250‑1 (MQ=255)                     aaGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGcaca                                                                                                                                                                                               <  2:30569/251‑1 (MQ=255)                                   aTGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATACAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCCAACAAGGAAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTTGCCTGGATATCGGTGTAACTTCTGGCGCACTGCCGTTTATTTCAGATGAATTCCAGATTACTTCGCACACGCAATAATTGGtt                                                                                                                                                                                 >  1:95150/1‑250 (MQ=255)                                       gATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTa                                                                                                                                                                              >  1:510575/1‑250 (MQ=255)                                           tGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGGCGTAAGCTc                                                                                                                                                                         >  1:334300/1‑251 (MQ=255)                                               tCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCat                                                                                                                                                                      <  1:302560/250‑1 (MQ=255)                                                  aaTCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCatgatg                                                                                                                                                                  <  2:145407/250‑1 (MQ=255)                                                       ttttACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTC‑AGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTc                                                                                                                                                               <  1:244972/248‑1 (MQ=255)                                                          tACGTCGTACTCACCTATCTTAATTTACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTg                                                                                                                                                           <  2:176527/250‑1 (MQ=255)                                                          tACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTg                                                                                                                                                           <  2:255579/250‑1 (MQ=255)                                                             gtcgtACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGCATGGGTCGTAAGCTCCATGATGTTCGGTGCGGc                                                                                                                                                       >  1:205375/1‑251 (MQ=255)                                                                     cACCTATCTTAATTCACAATAAAAAATAACCAGATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTcggt                                                                                                                                                <  2:653486/250‑1 (MQ=255)                                                                                aTTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGGGTAACTTCCGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGGAAGCTCCATGATGTTCGGGGCGGCAGTCGGTGCGGGGGGCAgg                                                                                                                                    >  2:50370/1‑250 (MQ=255)                                                                                  tCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGTTGCGGTGGGCAGCg                                                                                                                                   <  1:155066/250‑1 (MQ=255)                                                                                     caATAAAAAATAACCATATTGGAGGGCATCATGCCTGACTCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCggctg                                                                                                                               <  2:374821/251‑1 (MQ=255)                                                                                           aaaaTAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGCTTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTTGGTGCGGCAGTCGGTGCGGTTGGCAGCGGCTGGCtctc                                                                                                                         >  1:322246/1‑251 (MQ=255)                                                                                           aaaaTAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGGGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGTTGCGGCAGGCGGTGCGGTGGGCAGCGGCTGGCtctc                                                                                                                         >  1:196810/1‑251 (MQ=255)                                                                                            aaaTAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTcc                                                                                                                        <  1:149264/251‑1 (MQ=255)                                                                                             aaTAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCt                                                                                                                       >  1:187148/1‑251 (MQ=255)                                                                                                   cATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAAc                                                                                                                 <  2:80535/251‑1 (MQ=255)                                                                                                      atTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCCTTGGCCTGGATATCGGTTTAACTGCGTGCGCAATGCCGGTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTAGTAAGCTTCATGAAGTTCGGTGCGGCAGGCGGTGCGGTGGGCCGCGGCCGGCTCCCCTTTTAACTCg                                                                                                              >  2:837/1‑251 (MQ=255)                                                                                                        tGGAGGGCATCATGCCTGCCGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCCTGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCggg                                                                                                            >  1:573377/1‑251 (MQ=255)                                                                                                          gAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGgcg                                                                                                          >  1:190460/1‑251 (MQ=255)                                                                                                             ggCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCa                                                                                                        >  1:138820/1‑250 (MQ=255)                                                                                                               catcatGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCCGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGGGCGGTGGGCAGCGGATGGCTCTCCTTTAAACTCGGGCGCaaa                                                                                                      >  2:56550/1‑250 (MQ=255)                                                                                                                atcatGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGCTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCaaaa                                                                                                     >  1:272145/1‑250 (MQ=255)                                                                                                                   atGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAagag                                                                                                 <  1:484963/251‑1 (MQ=255)                                                                                                                     gCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGTATGACTCTTTGTCCTGGATATCGGTGTAACTGCTGGCTCACTGCCGTTTATTGCAGATGa                                                                                                                                                                                                                   <  2:301907/135‑1 (MQ=255)                                                                                                                     gCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGGTTACTCTTTAGCCTGGCTATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGa                                                                                                                                                                                                                   >  1:301907/1‑135 (MQ=255)                                                                                                                     gCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGGTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGGGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCCCTCCTTTAAACTCGGGCGCAAAAAGAGcc                                                                                               >  2:656428/1‑251 (MQ=255)                                                                                                                               aaaaaaCAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTTGCCTGGATATCGGTGTAACTGCTGGCGCCCTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAACTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCgg                                                                                     >  1:419611/1‑251 (MQ=255)                                                                                                                                     cAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAAt                                                                               >  1:130306/1‑251 (MQ=255)                                                                                                                                            ggTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATtttgttt                                                                        <  1:661604/251‑1 (MQ=255)                                                                                                                                            ggTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATtttgttt                                                                        <  1:414873/251‑1 (MQ=255)                                                                                                                                            ggTCAAACAAGGCAATGACGTTTTTAGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATtttgttt                                                                        >  2:99772/1‑251 (MQ=255)                                                                                                                                                        cAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATTTTGTTTGTTGCCGGTTCg                                                            >  1:56454/1‑251 (MQ=255)                                                                                                                                                        cAATGACGTTTTTCGTCTGATTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTCTTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCAGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTTATGATCGGCGCAATTTTTTTTTTCGCCGGTTCg                                                            >  2:259751/1‑251 (MQ=255)                                                                                                                                                              cGTTTTTCGTCTTCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCGGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGAACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTAGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAACAAGAGCCTGATGATAGGCGCAATTTTGTTTTGTGCCGGATCGCTggtt                                                      >  2:356265/1‑248 (MQ=255)                                                                                                                                                                 ttttCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATTTTGTTTGTTGCCGGTTCGCTGTtctc                                                    <  1:515204/250‑1 (MQ=255)                                                                                                                                                                              ccNTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATTTTGTTTGTTGCCGGTTCGCTGTTCTCTGCGGCTGCGCCa                                       <  2:592662/250‑1 (MQ=255)                                                                                                                                                                                 tGCCGCTCTGGCTGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATTTTGTTTGTTGCCGGTTCACTGTTCTCTGCGGCTACGCCAAAc                                    <  2:322518/250‑1 (MQ=255)                                                                                                                                                                                  gCCGCTCTGGCGGGATTACTCGTGGGCCTGGATATCGGTGTAACTGCTGGCGCACTTCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACGCGGGCGCAAAAAGAGCCTGATGATCGGCGCAAGTTTGTTGGGTGCCGGTTCGCTGTGCTCTGCGGCTGCGCCAAACGt                                  <  2:95150/251‑1 (MQ=255)                                                                                                                                                                                    cGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATTTTGTTTGTTGCCGGTTCGCTGTTCTCTGCGGCTGCGCCAAACGTTg                                >  1:321807/1‑251 (MQ=255)                                                                                                                                                                                        ctGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATTTTGTTTGTTGCCGGTTCGCTGTTCTCTGCGGCTGCGCCAAACGTTGAAGt                            >  1:200073/1‑251 (MQ=255)                                                                                                                                                                                                        ttGGCCTGGATATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATTTTGTTTGTTGCCGGTTCGCTGTTCTCTGCGGCTGCGCCAAACGTTGAAGTACTGATTCTTTcc               <  1:56550/248‑1 (MQ=255)                                                                                                                                                                                                           gCCTGGATATAGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGAGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATTTTGTTTGTTGCCGGTTCGCTGTTCTCTGCGGCTGCGCCAAACGTTGAAGTACTGGTTCTTTCCCGCGtt         >  1:107890/1‑251 (MQ=255)                                                                                                                                                                                                               ggANATCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGGCGTAAGCTCCATGAGGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATTTTGTTTGTTGCCGGTTCGCTGTGCTCTGCGGCTGCGCCAAACGTTGAAGTACTGATTCTTTCCCGCGTTCTAc     <  2:600447/251‑1 (MQ=255)                                                                                                                                                                                                                  tatCGGTGTAACTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTAGTAAGCTCCATGCTGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATTTTGTTTGTTGCCGGTTCGCTGTTCTCTGCGGCTGCGCCAAACGTTGAAGTACTGATTCTTTCCCGCGTTCTACTgg  >  1:64803/1‑251 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                ACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGGTGCGGTGGGCAGCGGCTGGCTCTCCTTTAAACTCGGGCGCAAAAAGAGCCTGATGATCGGCGCAATTTTGTTTGTTGCCGGTTCGCTGTTCTCTGCGGCTGCGCCAAACGTTGAAGTACTGATTCTTTCCCGCGTTCTACTGG  >  NC_000913/3088170‑3088629

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