Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,103,240	G→A	L76L (TTG→TTA)	mutY →	adenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,103,240	0	G	A	98.5%	222.3 / NA	65	L76L (TTG→TTA)	mutY	adenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/32);  minor base C (0/1);  total (32/33)
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.96e-01

GTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGA  >  NC_000913/3103007‑3103481                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   gTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGCACGCTTTATGGCGCGCTTCCCGACGGTGACCGATATCGCCAATGCGCCGCTCGCCGCAGCTCTCCACTTATGGACCGGGCTTGGCtc                                                                                                                                                                                                                                   >  2:147821/1‑250 (MQ=255) gTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCt                                                                                                                                                                                                                                    >  1:128446/1‑250 (MQ=255)  ggACCATGCAAGCTTCGCAATTTTCAGCCCATGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACGTCTGGACCGGGCTTGGCTAt                                                                                                                                                                                                                                  <  2:410198/250‑1 (MQ=255)    aCCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTa                                                                                                                                                                                                                                <  2:262710/251‑1 (MQ=255)                   cAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCaa                                                                                                                                                                                                                 <  1:126759/251‑1 (MQ=255)                       tttCAGCCCAGGTTCTGGACTGGTACGATCAATACGGGCGAAACACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGCCCGGGCTTGGCTATTACGCCCGCGCGCGCCATCTg                                                                                                                                                                                                             >  2:59869/1‑251 (MQ=255)                       tttCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCt                                                                                                                                                                                                              <  1:195038/250‑1 (MQ=255)                          cAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCa                                                                                                                                                                                                           <  1:102372/250‑1 (MQ=255)                                               aCGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTGATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGTCCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGc                                                                                                                                                                                      <  1:484695/250‑1 (MQ=255)                                                cGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTTACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCg                                                                                                                                                                                     >  1:354053/1‑250 (MQ=255)                                                      aTACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTac                                                                                                                                                                              <  1:494614/251‑1 (MQ=255)                                                      aTACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTac                                                                                                                                                                              >  1:159452/1‑251 (MQ=255)                                                        aCGGGCGAAAAACTCTGCCCTGGCAAATTGANAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAATCTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTGTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTacac                                                                                                                                                                            <  2:238983/251‑1 (MQ=255)                                                              gAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGCGGCGACCTTACAc‑ggcggg                                                                                                                                                                     >  1:318825/1‑250 (MQ=255)                                                                      ctGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACGCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTc                                                                                                                                                               <  1:326135/249‑1 (MQ=255)                                                                        gCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTccc                                                                                                                                                             <  1:139203/249‑1 (MQ=255)                                                                           cTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGaa                                                                                                                                                         >  1:396057/1‑250 (MQ=255)                                                                                 aaTTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCCGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTT‑CACCGGCGGTAACTTCCCGGAAACCTTTg                                                                                                                                                  >  2:209010/1‑251 (MQ=255)                                                                                 aaTTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTg                                                                                                                                                  >  2:210001/1‑251 (MQ=255)                                                                                      aCAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCg                                                                                                                                                            >  1:182442/1‑236 (MQ=255)                                                                                      aCAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCg                                                                                                                                                            <  2:182442/236‑1 (MQ=255)                                                                                       cAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCttttcgga                                                                                                                                              >  2:165734/1‑244 (MQ=255)                                                                                          gACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTg                                                                                                                                         >  2:12228/1‑251 (MQ=255)                                                                                               ccTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGTCGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAATCCTTTGAGGAAGTTgcagca                                                                                                                                    <  1:447097/251‑1 (MQ=255)                                                                                                 tACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCCCAACAAGTGGAGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGGTGCAGCAcc                                                                                                                                  >  2:187694/1‑250 (MQ=255)                                                                                                  aCAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCCCTg                                                                                                                                 >  2:494200/1‑251 (MQ=255)                                                                                                  aCAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTg                                                                                                                                 >  2:356951/1‑251 (MQ=255)                                                                                                  aCAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTg                                                                                                                                 >  1:330309/1‑251 (MQ=255)                                                                                                                 cAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCcgggcgtcgggcgt                                                                                                                  >  1:313849/1‑251 (MQ=255)                                                                                                                 cAGAAGGGATGTTGAAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCcgggcggcgggcg                                                                                                                   <  2:41971/250‑1 (MQ=255)                                                                                                                  aGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCTATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGtt                                                                                                                 <  1:461599/251‑1 (MQ=255)                                                                                                                  aGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTa                                                                                                                                                                    <  1:225805/200‑1 (MQ=255)                                                                                                                  aGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTa                                                                                                                                                                    >  2:225805/1‑200 (MQ=255)                                                                                                                        gATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCCATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTGGCAGCACTGCCGGGCGTCGGGCGTCCCACCg                                                                                                           >  1:481232/1‑251 (MQ=255)                                                                                                                         aTGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGc                                                                                                          >  1:13353/1‑251 (MQ=255)                                                                                                                                  caaACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGAt                                                                                                 <  2:121593/251‑1 (MQ=255)                                                                                                                                    aaCTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATtc                                                                                               <  1:347214/251‑1 (MQ=255)                                                                                                                                             gtcGACCGTTATCCCCTATTTTGATCGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCttt                                                                                      <  2:354053/249‑1 (MQ=255)                                                                                                                                                gACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATTTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGCTTCTCTCGCCTtctc                                                                                   >  2:48197/1‑251 (MQ=255)                                                                                                                                                              tttttGACCGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGGTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACt                                                                       <  2:396057/248‑1 (MQ=255)                                                                                                                                                                 ttGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGGTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCg                                                                  >  1:309128/1‑251 (MQ=255)                                                                                                                                                                    aaCGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGAtt                                                               <  1:187694/251‑1 (MQ=255)                                                                                                                                                                    aaCGCTTTATGGCGCGCTTCCCGACGGCGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTCAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCCCTCGCTTTCTCggggcactgcatttccgttt                                                               >  2:4277/1‑235 (MQ=255)                                                                                                                                                                        cTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCg                                                           >  1:333972/1‑251 (MQ=255)                                                                                                                                                                             tGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTACGCACTTTCCGATTCTCGACggg                                                      >  2:91185/1‑250 (MQ=255)                                                                                                                                                                              ggCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCCCTTTCCGCTTCTCGACGGt                                                      >  2:326678/1‑250 (MQ=255)                                                                                                                                                                               gcgcgcTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTaa                                                    <  2:171693/251‑1 (MQ=255)                                                                                                                                                                                         cGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAAGCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAAcgc                                          <  1:59869/251‑1 (MQ=255)                                                                                                                                                                                             ggTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTtct                                                                                    >  2:367970/1‑205 (MQ=255)                                                                                                                                                                                             ggTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTtct                                                                                    <  1:367970/205‑1 (MQ=255)                                                                                                                                                                                                     aTCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGgcgcgc                              >  2:298915/1‑251 (MQ=255)                                                                                                                                                                                                         cGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCgctg                            <  1:263500/249‑1 (MQ=255)                                                                                                                                                                                                         cGCCAATGCGCCGCTCGACGAAGTTCCCCACTTATGGACCGGGCTTGTCTATTACGCCCGCGCGCGCATTCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCTGAAACCTTTGAGGAAGTTGCAGCACTGCCGTGGGTCGGGCGTTCCACCGCAGTCGCGATGCTCTCGCTTTCTCTGGGTACGCACTTTCTGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCgct                             <  1:296095/248‑1 (MQ=255)                                                                                                                                                                                                          gCCAATGCGCCGCACGACAAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCACGCAATCTGCATAAAGCGGCACAACAAGTGGAGACCTTAAAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGGCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCCTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGCGCTGGAGCgctgctc                         >  2:126761/1‑250 (MQ=255)                                                                                                                                                                                                             aaTGCGCCGCGCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAAACTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGc                      <  2:128446/251‑1 (MQ=255)                                                                                                                                                                                                              aTGCGCCGCTCGACGAAGTTCTCCACTTAGGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCt                     <  2:89376/251‑1 (MQ=255)                                                                                                                                                                                                                gcgcCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCCCCGCAGGCGCGATTCTCTCGCTTTCTCTGGGGAAGCACTTTCCGATTCTCGACGGTAACNTCAAACGCGTGCTGGCGCGCTGCTATGCTGt                   >  2:244304/1‑251 (MQ=255)                                                                                                                                                                                                                    cgcTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGc               <  1:392750/251‑1 (MQ=255)                                                                                                                                                                                                                        cgacgaAGTTCTCCACTTATGGACCTGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGTCGTCTGGCGTTCCTCCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCt           <  1:168994/251‑1 (MQ=255)                                                                                                                                                                                                                        cgacgaAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTCTGCTGTAAGCGGCt           >  2:465163/1‑251 (MQ=255)                                                                                                                                                                                                                          acgaAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCGACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCt                                     <  2:150907/223‑1 (MQ=255)                                                                                                                                                                                                                          acgaAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCt                                     >  1:150907/1‑223 (MQ=255)                                                                                                                                                                                                                            gaAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGcc       <  2:313849/251‑1 (MQ=255)                                                                                                                                                                                                                               gTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTgg    <  1:6806/251‑1 (MQ=255)                                                                                                                                                                                                                                    ccACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGa  <  1:12228/248‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   GTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACAC‑GGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGA  >  NC_000913/3103007‑3103481

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