Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,749,123	C→T	G146D (GGT→GAT)	ydhT ←	uncharacterized protein YdhT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,749,123	0	C	T	97.8%	153.5 / NA	45	G146D (GGT→GAT)	ydhT	uncharacterized protein YdhT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/22);  minor base G (0/1);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTACCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGACTGCAAAATGGTCTCTCTCCGGGGTCTGTGGCCAGTTTTTGCATAGAGGAAACTTCCTGATGCGC  >  NC_000913/1748883‑1749367                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     ccGAGAATAATTGAGTCGACCACTCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACACACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAGGAATTATCGGTTACAGc                                                                                                                                                                                                                                            >  1:101733/1‑251 (MQ=255)              aGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGt                                                                                                                                                                                                                               >  2:52814/1‑251 (MQ=255)               gTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGtt                                                                                                                                                                                                                              >  1:149994/1‑251 (MQ=255)                              gCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCAGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCa                                                                                                                                                                                                                >  1:249516/1‑250 (MQ=255)                                          tGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTa                                                                                                                                                                                                    <  2:174508/250‑1 (MQ=255)                                               tttCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCCCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGCCACCGTTCACCCAATGTACAAACGCGCACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCg                                                                                                                                                                                              >  1:258969/1‑251 (MQ=255)                                                      gCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTc                                                                                                                                                                                       <  1:48834/251‑1 (MQ=255)                                                       cTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTcc                                                                                                                                                                                      >  2:320705/1‑251 (MQ=255)                                                         ggCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTTGCGATTTGCCCCAGGTCACCCTTACTGGCGACATGACACCGTTCACCCAAGGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCaa                                                                                                                                                                                    <  2:1118/251‑1 (MQ=255)                                                          gCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAAt                                                                                                                                                                                   <  1:88634/251‑1 (MQ=255)                                                                       cATTGTCGCCGCTGTTGGCACTTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCtata                                                                                                                                                                       >  1:148335/1‑250 (MQ=255)                                                                        aTTGTCGCCGCCATCGGCAATTCTGGCGTTAGCGCAAAAGCTTGGGATTTGCACCAGGTCATCATTACTGGCGCCATGACACCGTTCACACAATGTACTAACGCGTACACGCCTTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCTAAGTATTATCGGGTACAGCAACCCGTGCGCGTTCTCCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCACTccttattctatacc                                                                                                                                                                     <  2:272780/251‑13 (MQ=255)                                                                          tGTCGCCGATGTTGGCAATTCTGGAGTTTGCGCCAAAGCTGGAGATTTGAACAAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAAccc                                                                                                                                                                   >  2:203657/1‑251 (MQ=255)                                                                              gccgcTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATAc                                                                                                                                                               >  2:290198/1‑251 (MQ=255)                                                                                           aTTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTgcgc                                                                                                                                                  >  1:306806/1‑251 (MQ=255)                                                                                              cTGGCGTTTGCGCCAAAACTGGCGATTTGCACCAGTTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGCATTATCGGTTACAGCAACCAGTGCGCGTTCACCGGAGTTGCCACTGTGACTGTTAAAGCCGACCCTTCCAATCAGCACGCTATAACCCATACTGGCAATTTGCGCGGc                                                                                                                                               >  2:301383/1‑251 (MQ=255)                                                                                               tGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGCTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGc                                                                                                                                               >  2:23441/1‑250 (MQ=255)                                                                                                 gCGTTTGCGCCAAAGCTTGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCGCGCTATAACCCATACTGGCAATTTGCGCGGCATa                                                                                                                                            <  1:216021/251‑1 (MQ=255)                                                                                                     ttGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATAccc                                                                                                                                         <  1:22421/250‑1 (MQ=255)                                                                                                                   gCGNTTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATCCCCAGCGGTTGCGATTCt                                                                                                                          <  2:377115/251‑1 (MQ=255)                                                                                                                    cGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTg                                                                                                                         <  2:244582/251‑1 (MQ=255)                                                                                                                               gggtatcctgtccTGGGGACATGACACCCTTCACCCAATGTACAAACGCGTGCTCGCATTCCGGTATTCGCAATCCGACTGCTATCACCAGCGCTGCTTATCCGCGAAGAACTAGCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACGCTTCCAAGCCGCACGCTATAACCCATACGGGCAAGTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTc                                                                                                              <  2:101733/239‑1 (MQ=255)                                                                                                                               aGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTc                                                                                                              <  1:28106/251‑1 (MQ=255)                                                                                                                               aGGNCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTTCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGTCGTTATCCTc                                                                                                              <  2:381646/251‑1 (MQ=255)                                                                                                                                      cATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAg                                                                                                       <  2:306806/251‑1 (MQ=255)                                                                                                                                          aCTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCg                                                                                                    <  1:52814/250‑1 (MQ=255)                                                                                                                                                       cacCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGAtc                                                                                      <  1:238944/251‑1 (MQ=255)                                                                                                                                                         ccGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGAtcgt                                                                                    >  1:16105/1‑251 (MQ=255)                                                                                                                                                              cACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCTCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCt                                                                                <  2:172548/250‑1 (MQ=255)                                                                                                                                                                 ccAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCgg                                                                            >  1:154914/1‑251 (MQ=255)                                                                                                                                                                    aTGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGc                                                                          >  2:301712/1‑250 (MQ=255)                                                                                                                                                                                        aTTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGACTGCAAAATGGtct                                                     >  1:18935/1‑251 (MQ=255)                                                                                                                                                                                        aTTCCGGTATTCCAAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCCATACGCACGCTCTAACCCATACTAGCAATTTTCGCGGCATACCCAGCGGTTGAGATTCTGGCGGTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGAGCAAACTGACTGCAAACTGGtct                                                     >  2:222581/1‑251 (MQ=255)                                                                                                                                                                                                         ccGACTGCGATCTCCAGCACTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTGCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGACTGCAAAATGGTCTCTCTCCGGGGTCtgtg                                     <  2:258969/250‑1 (MQ=255)                                                                                                                                                                                                              tGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCTGTTACAGCAACCCGACCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGCCACTCCCAATCCGATCGATATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGAAAGTCGCTGCGGGCAAACTGACTGAAAAATGGTCTCTCTACGGCATCTGTGGCCAGt                               >  2:265044/1‑251 (MQ=255)                                                                                                                                                                                                               gCTATCTCCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGACTGCAAAATGGTCTCTCTCCGGGGTCTGTGGCCAGtt                              >  2:344497/1‑251 (MQ=255)                                                                                                                                                                                                                     tcCAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGAGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGAATGCAAAATGGGCTCTCTCCGGGGTCTGTGGCCCGTTTTTGc                         >  2:253712/1‑250 (MQ=255)                                                                                                                                                                                                                      ccAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGACTGCAAAATGGTCTCTCTCCGGGGTCTGTGGCCAGTTTTTGCAt                       <  1:320705/251‑1 (MQ=255)                                                                                                                                                                                                                       cAGCGCTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGCGCTGCCAATGTTACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGACTGCAAAATGGTCTCTCTCCGTGGTCTGTGGCCAGTTTTTGCATa                      <  2:249516/251‑1 (MQ=255)                                                                                                                                                                                                                          cgcTGCTTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGACTGCAAAATGGTCTCTCTCCGGGGTCTGTGGCCAGTTTTTGCATag                     >  2:197237/1‑249 (MQ=255)                                                                                                                                                                                                                              gctTATCCGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGACTGCAAAATGGTCTCTCTCCGGGGTCTGTGGCCAGTTTTTGCATAGAGGa                 <  1:64998/249‑1 (MQ=255)                                                                                                                                                                                                                                   tCCGCGCAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAAACCATACTGGCAATTTGCGCGGCATACCCAGCGGTCGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGCAAAACTGACTGCAAAATGGTCTCTCTCCGGGGTCTGTGGCCAGTTTTTGCATAGAGGAAACTTcc          >  2:266266/1‑251 (MQ=255)                                                                                                                                                                                                                                    ccGCGAAGAATTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGACTGCAAAATGGTCTCTCTCCGGGGTCTGTGGCCAGTTTTTGCATAGAGGAAACTTcc          <  1:137037/250‑1 (MQ=255)                                                                                                                                                                                                                                         aagaaTTATCGGTTCAATCAACCCGTTCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAAGCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGACTGCAAAATGGTCTCTCTCCGGGGTCTGTGGCCAGTTTTTGCATAGAGGAAACTTCCTGATgc    <  2:282662/251‑1 (MQ=255)                                                                                                                                                                                                                                            aaTTATCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCTCAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGACTGCAAAATGGTCTCTCTCCGGGGTCTGTGGCCAGTTTTTGCATAGAGGAAACTTCCTGATgcgc  <  1:23441/250‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     CCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGACTGCTATCTCCAGCGCTGCTTATCCGCGAAGAATTACCGGTTACAGCAACCCGTGCGCGTTCACCGGAGCTGCCAATGTGACTGTTAAAGCCGACACTTCCAATCCGCACGCTATAACCCATACTGGCAATTTGCGCGGCATACCCAGCGGTTGCGATTCTGGCGTTATCCTCACAACAGCCGATGCGAGCATGATCGTCGCTGCGGGCAAACTGACTGCAAAATGGTCTCTCTCCGGGGTCTGTGGCCAGTTTTTGCATAGAGGAAACTTCCTGATGCGC  >  NC_000913/1748883‑1749367

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