Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,937,210	C→T	L543F (CTC→TTC)	fucI →	L‑fucose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,937,210	0	C	T	98.7%	254.2 / NA	74	L543F (CTC→TTC)	fucI	L‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/43);  minor base G (0/1);  total (30/44)
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.91e-01

GACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTCTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTAT  >  NC_000913/2936963‑2937444                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           gACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTCCTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCg                                                                                                                                                                                                                                         >  1:29903/1‑251 (MQ=255)  aCTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGACTTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGc                                                                                                                                                                                                                                        >  1:47315/1‑251 (MQ=255)       tgtCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTcc                                                                                                                                                                                                                                    >  1:289448/1‑250 (MQ=255)          cAACATCATCACAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCAtgc                                                                                                                                                                                                                                >  2:64829/1‑251 (MQ=255)                   caAAGGCATGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATcc                                                                                                                                                                                                                       >  2:247894/1‑251 (MQ=255)                     aaGGCCTGGGACCGGGACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTCTCCCg                                                                                                                                                                                                                     >  2:71727/1‑251 (MQ=255)                        gCCTGGTACCGGTACTGCAAATCGCGGAAGGCTGGAGCGGGGAAGTGCCGAAGGATGTTCATGACAGCCTTAACAAACGCACCAACTCAACCTGGCCAGCGTCCTGGTTTGCCCCGCGCCTCCCCGGTAAAGGGCCGTTTACGGATGTGTACTCGGGAATGGCGAACTGGGGTTCTAACCATGGGGTTCTGACCAGCGGCTACGTTGGCGCATACTGTATCACCTTCGCCGCCATGCTGCTTAGCCCGgta                                                                                                                                                                                                                  <  2:353468/251‑1 (MQ=255)                            gggACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGgtatgt                                                                                                                                                                                                               >  2:495079/1‑250 (MQ=255)                             ggACCGGTACTGCAAATCGCGGACTCCTGGAGCGTGGAATTGCCGAAGGATGGGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGAGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTAGCCGCCACGCTGCGTAACACGgcatgtat                                                                                                                                                                                                             >  2:210455/1‑251 (MQ=255)                             ggACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGgtatgtat                                                                                                                                                                                                             <  1:572659/251‑1 (MQ=255)                               aacGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGc                                                                                                                                                                                                           >  1:231150/3‑251 (MQ=255)                                     aCTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACg                                                                                                                                                                                                     <  1:78122/251‑1 (MQ=255)                                        ggaaaTCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAAGAAACGCCCCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGTCGCATGCTTTATCTCTGTCGCCTCCATGCTGCGTATCCCGGTAGGTATGCACAACGTTg                                                                                                                                                                                                  <  2:538000/249‑1 (MQ=255)                                           aaTCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAag                                                                                                                                                                                               <  1:132577/251‑1 (MQ=255)                                              cgcgGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAag                                                                                                                                                                                               >  1:335916/1‑248 (MQ=255)                                                 ggAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGAcc                                                                                                                                                                                          >  1:94568/1‑250 (MQ=255)                                                   aaGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCaa                                                                                                                                                                                        >  2:129298/1‑250 (MQ=255)                                                      gCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAgtg                                                                                                                                                                                    <  1:342436/251‑1 (MQ=255)                                                       cTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAgtgt                                                                                                                                                                                   >  1:506030/1‑251 (MQ=255)                                                             cGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTTGCCAACCACCTTGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCTCAGCGTTGAAGAGACCAAAGTGTATCGTc                                                                                                                                                                             <  1:180080/251‑1 (MQ=255)                                                                ggAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGTCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCtt                                                                                                                                                                          <  1:540895/251‑1 (MQ=255)                                                                   aTTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTg                                                                                                                                                                       <  1:463/251‑1 (MQ=255)                                                                   aTTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTg                                                                                                                                                                       >  1:10995/1‑251 (MQ=255)                                                                     ggCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCACCCTGGCCAACCGCCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCTCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGTCCAAAGTGTATCGTCCTTCTGcc                                                                                                                                                                     <  2:367011/250‑1 (MQ=255)                                                                       ccGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTg                                                                                                                                                                   >  1:178768/1‑251 (MQ=255)                                                                       ccGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTg                                                                                                                                                                   <  2:90521/251‑1 (MQ=255)                                                                           aGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCt                                                                                                                                                               <  1:441972/251‑1 (MQ=255)                                                                                  gCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCAc                                                                                                                                                         <  2:94568/250‑1 (MQ=255)                                                                                    aTGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGc                                                                                                                                                      <  1:108638/251‑1 (MQ=255)                                                                                         aTCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATgg                                                                                                                                                  >  1:189745/1‑250 (MQ=255)                                                                                                acaaACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGaa                                                                                                                                          <  1:171933/251‑1 (MQ=255)                                                                                                 caaACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGa                                                                                                                                           >  1:312771/1‑249 (MQ=255)                                                                                                   aaCGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGc                                                                                                                                       <  1:368927/251‑1 (MQ=255)                                                                                                    aCGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGcc                                                                                                                                      <  2:46565/251‑1 (MQ=255)                                                                                                      gCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAg                                                                                                                                    <  2:350235/251‑1 (MQ=255)                                                                                                        aCCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGa                                                                                                                                  <  1:188825/251‑1 (MQ=255)                                                                                                         ccAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGAt                                                                                                                                 >  2:103671/1‑251 (MQ=255)                                                                                                           aaCTCAACCTGGCCAGCCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTa                                                                                                                               <  1:387844/251‑1 (MQ=255)                                                                                                               cAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTAcc                                                                                                                             <  1:309895/249‑1 (MQ=255)                                                                                                               cAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTAcc                                                                                                                             <  1:543249/249‑1 (MQ=255)                                                                                                                  ccTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTAACGCGCt                                                                                                                        >  1:424849/1‑251 (MQ=255)                                                                                                                   cTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCtt                                                                                                                       >  2:572260/1‑251 (MQ=255)                                                                                                                    tGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTg                                                                                                                      >  1:397239/1‑251 (MQ=255)                                                                                                                            caccTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACt                                                                                                              <  1:495079/251‑1 (MQ=255)                                                                                                                             accTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACt                                                                                                              >  2:485468/1‑250 (MQ=255)                                                                                                                              ccTGTTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTAGCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCGGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCATAACTAc                                                                                                            <  2:244220/251‑1 (MQ=255)                                                                                                                              ccTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTAc                                                                                                            >  1:174022/1‑251 (MQ=255)                                                                                                                                  gTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTc                                                                                                        >  2:142175/1‑251 (MQ=255)                                                                                                                                     tGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGt                                                                                                     >  2:236415/1‑251 (MQ=255)                                                                                                                                         ccGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTa                                                                                                 <  2:150064/251‑1 (MQ=255)                                                                                                                                            cgcCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACa                                                                                               >  2:307545/1‑250 (MQ=255)                                                                                                                                                      gTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATAc                                                                                    <  1:268252/251‑1 (MQ=255)                                                                                                                                                      gTAAAGGGCAGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCAGCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGATGCGCACGGCAGGGATATTGAAGGACAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATAc                                                                                    <  2:133707/251‑1 (MQ=255)                                                                                                                                                          aGGGCCGTTTACGGCTGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACAtt                                                                                 >  2:442706/1‑250 (MQ=255)                                                                                                                                                          aGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTa                                                                                <  1:485468/251‑1 (MQ=255)                                                                                                                                                              ccGTTTCCGAATGGGTCCTCGGTATTGGCGACCTGGGCGGTTAACCATGGGGTCCTGCCCATCGGCCCCGTTGGGGCAGCCTTTATCCCTTCCGCCCCCTGGCTGCTTATCCCGGTATTTTTGCCAAACGTTAAAAAGACCAAAGTTTATCGTCCTTCTCCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTcc                                                                            <  1:132085/251‑1 (MQ=255)                                                                                                                                                                gTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCt                                                                          <  1:544801/251‑1 (MQ=255)                                                                                                                                                                        aTGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGTCATGTATATTGAAGGCCAGGATTACCGCGCTTGGCAGAACTACGGGCCGTTGTACAAGCGTTAATACATTACTCCCtgatgtgat                                                                  <  1:569214/251‑1 (MQ=255)                                                                                                                                                                            gtACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGccc                                                              <  1:512579/251‑1 (MQ=255)                                                                                                                                                                                        ggCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGTCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCgt                                                   <  2:369467/250‑2 (MQ=255)                                                                                                                                                                                            aaCTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTAcc                                              <  1:247894/251‑1 (MQ=255)                                                                                                                                                                                             aCTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTTGCTCCGGCTACCg                                             >  2:362061/1‑251 (MQ=255)                                                                                                                                                                                                        aaCCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAAcaa                                  <  2:508955/251‑1 (MQ=255)                                                                                                                                                                                                        aaCCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTAGGCGCAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACTGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCGTATGTGATGCCCGGTCGCTCCGGCTTCCGGGCCTGAAcaa                                  <  2:508980/251‑1 (MQ=255)                                                                                                                                                                                                         aCCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAAcaag                                 <  1:537213/251‑1 (MQ=255)                                                                                                                                                                                                                ggTTCTGACCATCGGCCACGTTGGCGCAGACGTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGt                          <  2:459717/251‑1 (MQ=255)                                                                                                                                                                                                                  ttCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTgg                        >  1:505193/1‑251 (MQ=255)                                                                                                                                                                                                                      gACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGtt                      <  1:220926/249‑1 (MQ=255)                                                                                                                                                                                                                          aTCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCgg                <  1:131389/251‑1 (MQ=255)                                                                                                                                                                                                                            cGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGAt              <  1:406896/251‑1 (MQ=255)                                                                                                                                                                                                                              gCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATaa            <  2:302844/251‑1 (MQ=255)                                                                                                                                                                                                                               ccACGTTGGCGAAGACTTTATCAATTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATa             >  2:29441/1‑249 (MQ=255)                                                                                                                                                                                                                                 aCGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAGCTACGGTCCGTTGTACAAGCGGTAATACATGACTGCCTGATGTGATGGCCGGTCGCTCCGGCTGCCGGGTCGGAACAAGCAAGAGGGGTGAGCCGGATAAGCa         <  2:29322/251‑1 (MQ=255)                                                                                                                                                                                                                                        cgcAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGAGAAGCAATGTTAt  <  1:444784/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           GACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTCTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTAT  >  NC_000913/2936963‑2937444

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