Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,893,647	C→T	Q331* (CAG→TAG)	ygcN →	putative oxidoreductase with FAD/NAD(P)‑binding domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,893,647	0	C	T	98.4%	191.8 / ‑4.5	63	Q331* (CAG→TAG)	ygcN	putative oxidoreductase with FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (34/27);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

GCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCCAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATC  >  NC_000913/2893400‑2893890                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    gatgACTCGCGGTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCCTGCACAGTAGGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGCCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTTACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGc                                                                                                                                                                                                                                                  <  2:48553/249‑1 (MQ=255)    gACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCAc                                                                                                                                                                                                                                                <  2:217515/250‑1 (MQ=255)             tAAAGCGCATCCGGCAGTTCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCg                                                                                                                                                                                                                                       <  1:53639/250‑1 (MQ=255)                aGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCaa                                                                                                                                                                                                                                    >  1:326870/1‑250 (MQ=255)                           cAGTGCGCCCTCTCATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGACGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGCCTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCCCCGCGAGCCGCAAAATCTGTTTCCg                                                                                                                                                                                                                        >  2:206178/1‑251 (MQ=255)                           cAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCCCCAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTccc                                                                                                                                                                                                                        >  2:400038/1‑250 (MQ=255)                            aGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGc                                                                                                                                                                                                                       >  1:367897/1‑251 (MQ=255)                                   ccGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACCAACGCTGATAGACGCCTGCTAGCCCCGCGGGCCGCAAAATCTGGTTcccatttctc                                                                                                                                                                                                                 >  2:253756/1‑242 (MQ=255)                                       ttATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATcac                                                                                                                                                                                                            >  1:425741/1‑251 (MQ=255)                                        tATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATcac                                                                                                                                                                                                            >  1:110320/1‑250 (MQ=255)                                        tATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATcac                                                                                                                                                                                                            >  1:207061/1‑250 (MQ=255)                                        tATCAAAAACACGGAACCCCTGGGGTATGGGGCGCAGCTGGGGCCAGAAGGTGGCTTGCACAGTAGGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATAGGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATcac                                                                                                                                                                                                            <  2:158626/250‑1 (MQ=255)                                         aTCAAAAACACGGACTCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCCCAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCCCAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTACGCTTTATCATcacac                                                                                                                                                                                                          >  2:454569/1‑250 (MQ=255)                                           cAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACg                                                                                                                                                                                                        >  1:19891/1‑251 (MQ=255)                                            aaaaaCACGGAATCCGTGGTGTATGGTGCGCTTCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCATTTGCGTCAATCCCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGt                                                                                                                                                                                                       <  2:129752/251‑1 (MQ=255)                                               aaCACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGt                                                                                                                                                                                                       >  1:281119/1‑248 (MQ=255)                                                                    ggTGCGCATCTGGTTCCAGAAGGTGGCTGGCCCAGTATGCCGGGTCAATCCGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATg                                                                                                                                                                                <  2:110938/250‑1 (MQ=255)                                                                    ggTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCTGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGt                                                                                                                                                                               >  1:350700/1‑251 (MQ=255)                                                                    ggTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTatcatc                                                                                                                                                                                                              >  1:251777/1‑220 (MQ=255)                                                                           aTCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGtt                                                                                                                                                                                                                            >  1:324627/1‑199 (MQ=255)                                                                           aTCTGGTGCCAGAAGGGGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGtt                                                                                                                                                                                                                            <  2:324627/199‑1 (MQ=255)                                                                             cTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCCCAAAAGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGTTGTTCTACAGCg                                                                                                                                                                      >  2:255223/1‑251 (MQ=255)                                                                               ggTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGCGCCGCAAAATCTGTTTCCGCTTTTTCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGtt                                                                                                                                                                    >  2:22503/1‑251 (MQ=255)                                                                                 tGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACACCGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTAt                                                                                                                                                                  >  2:318800/1‑251 (MQ=255)                                                                                    cAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATTACAACGTAGAGCGCTGCCTGCTGTGGGATGTTCTACAGCGTTATc                                                                                                                                                                 >  2:269687/1‑249 (MQ=255)                                                                                       aaGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCAt                                                                                                                                                            >  2:154067/1‑251 (MQ=255)                                                                                       aaGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCAt                                                                                                                                                            <  1:154067/251‑1 (MQ=255)                                                                                        aGGTGGCTTGCACAGTATGCCGGTGCAATTCGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATg                                                                                                                                                           <  1:400038/251‑1 (MQ=255)                                                                                        aGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGGGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGGGCGCCGCCTGCCGGGGGGTGTTCTAACGCGTTATCAGCATg                                                                                                                                                           >  2:430298/1‑251 (MQ=255)                                                                                         gggggtttGCACGGTAGGCCGTTGCAATACGCCTGTAGCGGCTGCCTGCTGGTGGGCGATGCGTTGCGCAGGTGCGTCAATACCGGAATTTCCGTGTGCGGCATGGATAGGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCGGATAACCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGt                                                                                                                                                          <  2:207061/245‑1 (MQ=255)                                                                                           tGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGATGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAACGTTATCAGCATGTTc                                                                                                                                                        >  1:207537/1‑251 (MQ=255)                                                                                           gggCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGGCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTc                                                                                                                                                        <  1:143820/250‑1 (MQ=255)                                                                                            ggCTGGCACATTAGGCCGGTGCAATACGCCGGTAACGGCTGGTTGCGGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATAGGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTcc                                                                                                                                                       <  2:293482/251‑1 (MQ=255)                                                                                             gCTTGCACAGTATGCCGGTGCAATCCGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAAAACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCg                                                                                                                                                      <  2:476182/251‑1 (MQ=255)                                                                                                tGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGgcg                                                                                                                                                   >  1:311537/1‑251 (MQ=255)                                                                                                    caGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCttt                                                                                                                                               <  1:271868/251‑1 (MQ=255)                                                                                                    caGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCttt                                                                                                                                               >  1:350574/1‑251 (MQ=255)                                                                                                                      cGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGt                                                                                                                             <  1:11364/251‑1 (MQ=255)                                                                                                                                      ctgctgGTGGGCGATGCGTTGCGCCGTTGTGTAATTACCGGAATTTCCGTGCGCGGCATGGATATGGCGATGACTGGCGCGCCGGCGGCCACACAAACGCTGATAAGCGCCTGCTAGCACCGAGAGCCGCAAAATCCGTTTCCGCTTTATCATCACAACGTAGCGCGCAGCCTGCTGTGGGAGGTTCTACAGCGTTATAAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGGACGAGGCCTGc                                                                                                             >  2:453447/1‑251 (MQ=255)                                                                                                                                            gtggGGGATGCGTTGCGCAGTTGCGTCAATCCCGGAATTTCCGTGCGGGGCATGGATATGGCGCTGACGGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAAt                                                                                                       <  2:281119/251‑1 (MQ=255)                                                                                                                                               ggCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAAAGCCAGGGATGGTACCGTACGTGGACTGCGTTAATGc                                                                                                     >  1:275836/1‑250 (MQ=255)                                                                                                                                                   attggTTGCGCAGTTGCGTCAATCCCGGAATTGCCGTGCGCGGCATGGATATGGCGCTGACTGGCGGGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGat                                                                                                <  2:311537/247‑1 (MQ=255)                                                                                                                                                                  gCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGAttt                                                                                  >  2:273736/1‑250 (MQ=255)                                                                                                                                                                    gtcacgaCCGGAATGTCCGTGCGCGGCATGGACATGGCGCTGACGGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCCCCGCGAGCCGCAAAATCTGTTTCCGCTTTAGCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGAtttctg                                                                               <  2:102638/245‑4 (MQ=255)                                                                                                                                                                         tACCGGCATTTCCGTGCGCGGCTTGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGGACCGTACGTGGCCTGCGTTATTGCAGGATATTTCCCGCGATTTATGGGAt                                                                           >  2:378441/1‑250 (MQ=255)                                                                                                                                                                          aCCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACGGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCCCCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCa                                                                         <  2:19891/251‑1 (MQ=255)                                                                                                                                                                                               aTGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCAcc                                                    >  1:481373/1‑251 (MQ=255)                                                                                                                                                                                               aTGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTAGCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCAcc                                                    <  1:454569/251‑1 (MQ=255)                                                                                                                                                                                               aTGGATAGGGCGCTGACTGGCGCGCGGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCAcc                                                    <  2:110320/251‑1 (MQ=255)                                                                                                                                                                                                 ggATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTaa                                                                                                        <  2:239467/197‑1 (MQ=255)                                                                                                                                                                                                 ggATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTaa                                                                                                        >  1:239467/1‑197 (MQ=255)                                                                                                                                                                                                             gACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTa                                       <  1:422167/250‑1 (MQ=255)                                                                                                                                                                                                                tGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCACGATTTATGGGATCAGGGTGCTAAACCTGTTCCACCGCTGCGCCAGTTATTCt                                   >  2:296185/1‑251 (MQ=255)                                                                                                                                                                                                                         ggcggcggcACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCAtt                          >  1:218591/1‑251 (MQ=255)                                                                                                                                                                                                                          gcggcggcACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCAttt                         <  1:329026/251‑1 (MQ=255)                                                                                                                                                                                                                           cggcggcACAAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTa                        <  1:238637/251‑1 (MQ=255)                                                                                                                                                                                                                           cggcggcACAAAATATGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCCTGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCAttt                         >  2:317844/1‑250 (MQ=255)                                                                                                                                                                                                                               ggcACAAACTCTGATAAGGGCCGGCGAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTAcgtc                    <  2:252116/251‑1 (MQ=255)                                                                                                                                                                                                                                 cacaAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGt                                                                                                                        <  2:141670/149‑1 (MQ=255)                                                                                                                                                                                                                                 cacaAACGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGt                                                                                                                        >  1:141670/1‑149 (MQ=255)                                                                                                                                                                                                                                    aaaCGCTGATAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCAc               <  1:433856/251‑1 (MQ=255)                                                                                                                                                                                                                                            aTAAGCGCCTGCTAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCtgtg       >  1:79861/1‑251 (MQ=255)                                                                                                                                                                                                                                                 cgcCTGCCAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATc  >  1:90248/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    GCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGATGCGTTGCGCAGTTGCGTCAATACCGGAATTTCCGTGCGCGGCATGGATATGGCGCTGACTGGCGCGCAGGCGGCGGCACAAACGCTGATAAGCGCCTGCCAGCACCGCGAGCCGCAAAATCTGTTTCCGCTTTATCATCACAACGTAGAGCGCAGCCTGCTGTGGGATGTTCTACAGCGTTATCAGCATGTTCCGGCGCTTTTGCAACGCCCGGGATGGTACCGTACGTGGCCTGCGTTAATGCAGGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATC  >  NC_000913/2893400‑2893890

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