Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,436,956	2 bp→TT	coding (202‑203/741 nt)	cysQ →	3'(2'),5'‑bisphosphate nucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,436,956	0	C	T	100.0%	199.0 / NA	55	P68F (CCC→TTC)	cysQ	3'(2'),5'‑bisphosphate nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/24);  total (31/24)
*	NC_000913	4,436,957	0	C	T	98.2%	190.4 / ‑6.3	55	P68F (CCC→TTC)	cysQ	3'(2'),5'‑bisphosphate nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (31/23);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATGC  >  NC_000913/4436722‑4437194                                                                                                                                                                                                                                            ||                                                                                                                                                                                                                                              aGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGCAGTCCg                                                                                                                                                                                                                                >  2:287224/1‑251 (MQ=255)          cGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTg                                                                                                                                                                                                                       >  1:207554/1‑251 (MQ=255)           gCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTgg                                                                                                                                                                                                                      >  1:127704/1‑251 (MQ=255)              cacaGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAg                                                                                                                                                                                                                   <  1:352733/251‑1 (MQ=255)               acaGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGc                                                                                                                                                                                                                  >  1:216080/1‑251 (MQ=255)                 agagACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGt                                                                                                                                                                                                                >  1:402832/1‑251 (MQ=255)                   agaCGCGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTATCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGCAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGtt                                                                                                                                                                                                               >  2:428608/1‑250 (MQ=255)                   agaCGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGtt                                                                                                                                                                                                               >  1:53931/1‑250 (MQ=255)                       gAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTActgg                                                                                                                                                                                                          <  2:413558/251‑1 (MQ=255)                       gAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGAGGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGGTActgg                                                                                                                                                                                                          <  2:166957/251‑1 (MQ=255)                         ggTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTActggct                                                                                                                                                                                                        <  2:233539/251‑1 (MQ=255)                          gTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGAGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACACTTCTCCGGTAACGGCAGCGGATATTGCCGCTCCCACCGTTATCCTGGACGGTTTACATACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑CGCCTTTCGGTTGGAAAGTAAGTCAGAACTGGCAGCGTTActggatg                                                                                                                                                                                                       >  1:448115/1‑248 (MQ=255)                                      tAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCt                                                                                                                                                                                            <  1:31340/250‑1 (MQ=255)                                            aaGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGt                                                                                                                                                                                     <  2:385326/251‑1 (MQ=255)                                              gTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTAc                                                                                                                                                                                   <  1:232841/251‑1 (MQ=255)                                                ctGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGGCGACCCGCCGGATGGTACt                                                                                                                                                                                  >  2:458098/2‑250 (MQ=255)                                                    cAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAaga                                                                                                                                                                             <  2:133916/251‑1 (MQ=255)                                                    cAGCTTGCACGGAATGCAGGCGATGCCATTAGGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACTCGCTGGATGGTACTAAaga                                                                                                                                                                             <  1:384598/251‑1 (MQ=255)                                                       cTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGtt                                                                                                                                                                          >  2:185748/1‑251 (MQ=255)                                                            aCGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGGTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTttatta                                                                                                                                                                     <  1:287224/251‑1 (MQ=255)                                                              ggAATGCAGGCGATGACATTATGCAAGTCTACCACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCCCCGGTACAGGCAGCGGATAATGCCGCTCACACCGTTATCCTGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCCGCGTTACTGGCTGGTAGACCCGCCGGGTGGTACTAAAGAGTTTATTaa                                                                                                                                                                    >  2:464279/1‑250 (MQ=255)                                                                         gATGCCATTATGCAGGCCTACGACGGGACGACACCGATGGACGTCGTCAGCACAGCGGACCATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCACGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATACTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCCCCGGGCGGTACAAAAGAAGTTATTAAACGTACTGGCg                                                                                                                                                         >  2:294539/1‑250 (MQ=255)                                                                            gCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAAtt                                                                                                                                                     <  1:310497/251‑1 (MQ=255)                                                                               aTTATGCAGGTCTACGACGGGACGAAACCGATGGCCGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGCTGGTACTAAAGAGTTTATTACACGTCATGGCGAATTCAc                                                                                                                                                  >  2:167652/1‑251 (MQ=255)                                                                                   tGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTACCGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGtt                                                                                                                                              >  1:67936/1‑251 (MQ=255)                                                                                      aGGTCGACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCTCTCACACCGTTATCAGGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTACTTTCTGAAGAAG‑AGGCTTTCGGTTGGGAAGTCCTTCAGCTCTGGCAGCGCGACTGGCTGGTAGACCGGCTGTATGGTACTAAAGAGTTTATCAAACGTAATGGCGAATTCCCCTTTAAc                                                                                                                                           <  2:221642/251‑1 (MQ=255)                                                                                        gtctccgccgGGCCGAAACCGATGGACTTCGTCAGCAAACCGGACATTTCTCCGTTAACGGCACCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTCCTTCCGCTGACACCGGTTGTCCCGTTCCTTCCTGA‑GAAGGTTCCTTCCGTTTGGGAGGTCCTTCACCCCTGGCGCCTTTACTGGTTGTTGGCCCCGCTGGATG‑TTCTTAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACAtt                                                                                                                                        <  2:2171/243‑1 (MQ=255)                                                                                            acgacgGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTgcg                                                                                                                                     <  1:350458/251‑1 (MQ=255)                                                                                                       aaaCCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGAcc                                                                                                                           >  1:322118/1‑250 (MQ=255)                                                                                                                 aCGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCg                                                                                                                <  1:92035/251‑1 (MQ=255)                                                                                                                   gtcgtcAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGAt                                                                                                              <  2:455337/251‑1 (MQ=255)                                                                                                                   gtcgtcAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGAt                                                                                                              <  1:214857/251‑1 (MQ=255)                                                                                                                        cAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAg                                                                                                         <  1:87908/251‑1 (MQ=255)                                                                                                                          gCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGc                                                                                                       <  1:376513/251‑1 (MQ=255)                                                                                                                                  gACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGAACATGGCAAACCGAttt                                                                                                            >  1:317199/1‑238 (MQ=255)                                                                                                                                  gACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGAGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGAttt                                                                                                            <  2:317199/238‑1 (MQ=255)                                                                                                                                   aCAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTAt                                                                                              >  1:454740/1‑251 (MQ=255)                                                                                                                                               tAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATg                                                                                  >  1:175461/1‑251 (MQ=255)                                                                                                                                                   ggCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGGAATGAACg                                                                              >  1:403813/1‑251 (MQ=255)                                                                                                                                                                   gCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGc                                                              >  2:192626/1‑251 (MQ=255)                                                                                                                                                                      cacaCCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCa                                                             <  1:228516/249‑1 (MQ=255)                                                                                                                                                                            gTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCa                                                      >  1:323114/1‑250 (MQ=255)                                                                                                                                                                            gTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCa                                                      >  1:112203/1‑250 (MQ=255)                                                                                                                                                                              tATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGTGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCaa                                                     >  2:445233/1‑249 (MQ=255)                                                                                                                                                                                         gTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAaaga                                         >  2:129568/1‑250 (MQ=255)                                                                                                                                                                                         gTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAaaga                                         >  1:379141/1‑250 (MQ=255)                                                                                                                                                                                                      gACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAAGGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGtgtg                            <  1:302740/250‑1 (MQ=255)                                                                                                                                                                                                         acCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATGCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGTCAAAGCCTGGAAAGAAGAGTGCGGTGTGCGCa                        <  2:291986/251‑1 (MQ=255)                                                                                                                                                                                                           cGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGGGGAAAGAAGAGTGCGGTGTGCGCa                        >  2:257736/1‑249 (MQ=255)                                                                                                                                                                                                            ggATGTTCCGGTCCTTTCTGAAGACG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCCTTGACCATGGCAAACCGATTTTAGGCGTGGTTTATGCGCCGGTAATGCACGTCATGTACCGCGCGGCAGAAGGCAAAGCGTGGACAGAAAAGTGCGGTGTGCGCAAGc                     >  2:412629/1‑251 (MQ=255)                                                                                                                                                                                                            ggATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGc                     >  1:184961/1‑251 (MQ=255)                                                                                                                                                                                                                    cGGTCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTc               >  1:192083/1‑249 (MQ=255)                                                                                                                                                                                                                       tCCTTTCTGAAGAAG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAAAGTAATGGCGAATTCACCGTTAACATTGCGCTCCTTGACCCTGGCAAACCGATTTTAGGCGTGGGGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAACGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTc          >  2:311085/1‑251 (MQ=255)                                                                                                                                                                                                                            tCTGAAGAAG‑ATCCTTTCGGTTGGGAAGGCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCTCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAATAGTGCGGTGTGCGCAAGCAGATTCAGGTCcgcg      <  1:428608/250‑1 (MQ=255)                                                                                                                                                                                                                               gaagaaG‑ATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCCTTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATTAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATgg  >  2:318556/1‑250 (MQ=255)                                                                                                                                                                                                                                            ||                                                                                                                                                                                                                                              AGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAG‑ATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATGC  >  NC_000913/4436722‑4437194

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