Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,577,937	T→C	E147E (GAA→GAG)	maeB ←	malate dehydrogenase (oxaloacetate‑decarboxylating) (NADP(+))

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,577,937	0	T	C	98.3%	206.7 / NA	59	E147E (GAA→GAG)	maeB	malate dehydrogenase (oxaloacetate‑decarboxylating) (NADP(+))
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (30/28);  minor base A (1/0);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CGGCCCTGATAGATAACGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACC  >  NC_000913/2577692‑2578174                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              cGGCCCTGATAGATAACTCCTTTTTAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCTGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCTTTGAGGATGGCGTCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGTTCGCGCATTTTCTGCTCAat                                                                                                                                                                                                                                          <  1:206617/251‑1 (MQ=255) cGGCCCTGATAGATAACGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAat                                                                                                                                                                                                                                          <  2:492589/251‑1 (MQ=255)       tGATAGATAACGCCTTTTGAATCGCAAAACACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATCATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGTAATATTCATCCGCTCGCGCAGTTTCTGCTCAAtataaa                                                                                                                                                                                                                                     >  2:420522/1‑248 (MQ=255)        gatagataACGCCTTTTGAAACGCAAACCACGATGTTATGTCTTTGCAGACCCAGCGCTACCAGCATGTTCAGACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGaaa                                                                                                                                                                                                                                   <  1:11842/251‑1 (MQ=255)               aaCGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCACGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCACTCGCGCAGTTTCTGCTCAATATAGAACCCTTCTg                                                                                                                                                                                                                            >  2:434248/1‑251 (MQ=255)                      tttGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCt                                                                                                                                                                                                                       <  2:186744/249‑1 (MQ=255)                         gtaccgataaccccgatatgatgttTTTGCAGACCCAGCGCTACCACCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGTAGATGTTTTTCGCCACCACGCGCCAGCCTTTGAGGTTGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATTCCGGAATATTCATTCGCTCGCTCAGTTTCTGCTCAATATAGAACCATTCTGGCTCTTTAAt                                                                                                                                                                                                                  <  2:273381/231‑1 (MQ=255)                             cGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCGCCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTGCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTc                                                                                                                                                                                                               <  1:375071/250‑1 (MQ=255)                               cAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTc                                                                                                                                                                                                            >  1:47057/1‑251 (MQ=255)                                    cACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACCGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGACGGAGATGTTTTTCTCCACCACGAGCAAGACGTTGAGAATGGCGGCAATGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGCATATTCATCCGCTCGCGCAGTTTATGCTCCATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGt                                                                                                                                                                                                       >  1:275389/1‑251 (MQ=255)                                        atgttatgttTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGAt                                                                                                                                                                                                   <  1:487670/251‑1 (MQ=255)                                            tatgttTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGCTGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCCATATAGAAACATTCTGGCGCTTTAATGTCTTCGCGGTTGTTgccg                                                                                                                                                                                               >  1:280698/1‑251 (MQ=255)                                             attgttttGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTTGGAGATGTTTTGCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTTCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATgccg                                                                                                                                                                                               <  1:383121/246‑1 (MQ=255)                                              tgttTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATgccgcc                                                                                                                                                                                             >  1:473118/1‑251 (MQ=255)                                                 ttttGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCCGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATCTTTTTCTCCACCACGCGCAACCCGTTGCGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTACTGTCTTCGAGGTTGATGCCGCCGa                                                                                                                                                                                           >  1:77489/1‑250 (MQ=255)                                                    tGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAggt                                                                                                                                                                                       >  2:398589/1‑251 (MQ=255)                                                         aCCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAggttggtt                                                                                                                                                                                  >  2:429584/1‑251 (MQ=255)                                                                gcTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAgcgc                                                                                                                                                                           >  1:195890/1‑251 (MQ=255)                                                                   aCCAGCAGGTTCAGACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGAGGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCGCGCGCAGTTTCTGCTCAATATAGAAACATGCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGcgcg                                                                                                                                                                          <  1:317651/249‑1 (MQ=255)                                                                     cagcagGTTCATACAGGCGATTGCTGCGGCACCAGGGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGAATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCCGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCgcggcg                                                                                                                                                                       >  2:262903/1‑250 (MQ=255)                                                                      agcagGTTCATACAGGCGATGGCTGCGGCTCCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGAc                                                                                                                                                                     <  2:47057/251‑1 (MQ=255)                                                                        cagGTTCATACAGGCGATTGCTGCGGCACACGCGCCGGAAACCACGATCCGCCCGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGATGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGCTGCCGCCGAAGGATGGTTCGCGCGCGGCGACac                                                                                                                                                                   >  2:386188/1‑250 (MQ=255)                                                                                       gATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGt                                                                                                                                                    >  1:43298/1‑251 (MQ=255)                                                                                              gcGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCggg                                                                                                                                               <  1:103079/249‑1 (MQ=255)                                                                                               cGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGGTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAAGAAATTTGTCCGGGTCg                                                                                                                                            <  1:162490/251‑1 (MQ=255)                                                                                                    cccGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTc                                                                                                                                       <  2:246332/251‑1 (MQ=255)                                                                                                    cccGCGCCGGAAACCACCATCCGCACGGCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAAGATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTc                                                                                                                                       <  1:57686/251‑1 (MQ=255)                                                                                                     ccGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCg                                                                                                                                      <  1:300600/251‑1 (MQ=255)                                                                                                      ggcgcCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGt                                                                                                                                     >  2:197079/2‑251 (MQ=255)                                                                                                        cgcCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTc                                                                                                                                    >  2:168047/1‑250 (MQ=255)                                                                                                           cGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACt                                                                                                                                >  2:397924/1‑251 (MQ=255)                                                                                                                 ccacccTCCGCACGTCGGAGATGTTTTTCGCCACCTCGCGCAATCCGTTGAGGATGGCGGCAGTGCCGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATGCAGCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGTCGAAGGTTGGTTCGACCGCGGCGACAACTTCAATAAATTTGTCCGTGTCGAGTTCGTCAACTTCAATg                                                                                                                          <  1:314333/251‑1 (MQ=255)                                                                                                                    ccaTCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCa                                                                                                                       <  1:31847/251‑1 (MQ=255)                                                                                                                          gctcGTCGGAGATGTTTTTCTCCGCCGCGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATGGCCGTGCCGTGCTGTTCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCTACAACTTCAATAAATTTGTCCGGGTCGTGTTCGTCAACTTCAATGTCAAAttca                                                                                                                 <  2:374549/248‑4 (MQ=255)                                                                                                                                    gaTGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCgg                                                                                                       <  2:288276/251‑1 (MQ=255)                                                                                                                                                         cAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGCGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGACCGc                                                                                  >  2:334073/1‑251 (MQ=255)                                                                                                                                                          aaGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGcc                                                                                 <  1:220747/251‑1 (MQ=255)                                                                                                                                                          aaGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCGCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACCGAACGcc                                                                                 <  2:195890/251‑1 (MQ=255)                                                                                                                                                           aGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAACCATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGCACGccc                                                                                >  2:166181/1‑251 (MQ=255)                                                                                                                                                              cGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTg                                                                             <  1:397924/251‑1 (MQ=255)                                                                                                                                                                 tGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCt                                                                          >  2:28172/1‑251 (MQ=255)                                                                                                                                                                 ggAGGATGGCGGCAGTGCTGATAATTGCCGGGCCGTGCTGATCGTCGTGGAAGACCGGAATATTCATCCGCTCGCGCAGTGTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTGCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATTCCGGCGAAGTTCTTAAACAGAACTCCCTTGCCt                                                                          <  1:386188/250‑1 (MQ=255)                                                                                                                                                                  gAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCtt                                                                         <  1:429584/251‑1 (MQ=255)                                                                                                                                                                      aTGGCGGCAGGGCTGATAATTGGCGTGCCGTGCTGATCGTCGTGGAAAACCGGAATATTCATCCGCACGCGCAGTTTCCGCTCCATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGCTGCCGCCGCAGGTTTGTTCGAGCGCGGCGACAACTTCAATAACTTTGTCCGGGTCGAGTTCTTCAACTTCAATGTCAAATACATCAATCCCGGCGACTTTCTTAAACAGAACGCCCTTGCCTTCCAt                                                                     >  2:60887/1‑251 (MQ=255)                                                                                                                                                                        ggcggcAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCCACTTCAATGTCAAATACATCAATCCCGGCGACTTTCTTAAACAGAACGCCCTTGCCTTCCATc                                                                    >  2:467550/1‑250 (MQ=255)                                                                                                                                                                           ggcAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCg                                                                <  1:251691/251‑1 (MQ=255)                                                                                                                                                                                tGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTtg                                                           >  1:293637/1‑251 (MQ=255)                                                                                                                                                                                  cTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTtgcc                                                         >  1:240651/1‑251 (MQ=255)                                                                                                                                                                                     aTAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTACACAGAACGCCCTTACCTTCCATCACCGGTTtgcctgc                                                      >  2:481104/1‑251 (MQ=255)                                                                                                                                                                                          tgccgtgccgTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAgcg                                                 >  2:113953/1‑251 (MQ=255)                                                                                                                                                                                            ccgtgccgtgCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAgcgcg                                               >  1:395313/1‑251 (MQ=255)                                                                                                                                                                                                     gcTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGtt                                       <  1:28172/250‑1 (MQ=255)                                                                                                                                                                                                             tcgtGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAAccc                              >  1:54291/1‑251 (MQ=255)                                                                                                                                                                                                             tcgtGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAAccc                              >  2:155457/1‑251 (MQ=255)                                                                                                                                                                                                                       ccGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCAccgcc                    <  2:514261/251‑1 (MQ=255)                                                                                                                                                                                                                          gAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGt                  >  1:269163/1‑250 (MQ=255)                                                                                                                                                                                                                           aaTATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTAc                <  1:268765/251‑1 (MQ=255)                                                                                                                                                                                                                                 catccgaacgcgCACTTACTGATCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGACGGTTGGTTCGAGCGCGGCGACACCTTCAATCAATTTGTCCGGGTCGCGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGtt            >  2:84326/9‑249 (MQ=255)                                                                                                                                                                                                                                                 ttCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCAcc  >  1:378174/1‑243 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              CGGCCCTGATAGATAACGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACC  >  NC_000913/2577692‑2578174

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: ‑