BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,752	T→A	G99G (GGA→GGT)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,752	0	T	A	96.5%	247.9 / NA	87	G99G (GGA→GGT)	araB	ribulokinase
Reads supporting (aligned to +/- strand): ref base T (1/0); major base A (47/36); minor base C (0/3); total (48/39)
Fisher's exact test for biased strand distribution p-value = 8.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TCCGCTTCTT‑CAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCCAC  >  NC_000913/69507‑69994

cacaGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TA‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACACT‑ACGAACATCGCGTTCGGGTTTTCGGCACACTCCGGGCGCAGCGCCAGCACGTTACCGTc                                                                                                                                                                                                                                                 >  1:147477/1‑251 (MQ=255)

 acaGCCCAATCCACGAGGCGGCAGATGGCCCCCCGGCGCTGTCCGGGCGAGTCACATGCAGGATTTTTGCCAAAACCCATCCGCGGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCg                                                                                                                                                                                                                                                <  1:293988/251‑1 (MQ=255)

               gATGCGGCAGAATGAGAAACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATGCGCTGGAATAAATACAGCCGATATAGCGGGAGTAGTCAACATTGCCCGGGGCGAGACACAAACGGTTCATCTAT‑TC‑GCATCTTTCAACAGAAGTGT‑GCACTTTTCCCAAA‑ACGAAAATCGTGTTAGGTGTTTCGTGAAACTATGGACGCAGCGCCAGAACGTTTACGTCGGGCTCACTCGGTGc                                                                                                                                                                                                                                  >  2:745035/1‑251 (MQ=255)

                     gCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACAGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCATCACGTTACCGTCGGCATCAATCGGTGCGGGCgg                                                                                                                                                                                                                            >  2:530673/1‑250 (MQ=255)

                      cAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCTTCAATCGGTGCGGGCG‑Tc                                                                                                                                                                                                                          >  2:93996/1‑251 (MQ=255)

                             ggccgcGGCGTTGTCGTGGCTGGTCACAGGCAGGTTTTTTGCCCAGACCCATCCGCGGAAATAAATCCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCGCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCgg                                                                                                                                                                                                                   <  1:789772/246‑1 (MQ=255)

                                     cgcTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGAGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATAAATCGGTGCGGGCG‑TCGAGCCGGTACTGTc                                                                                                                                                                                                            >  1:503330/1‑250 (MQ=255)

                                      gcTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGAGTGGCACAAACGGGTAATCTCT‑TC‑GCGTCTTTCAACAGCAGGGT‑GGGATTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCAAGTTACCGTGGGAATCAATCGGGCCGGGCG‑TCTAGCCGGTCCTggca                                                                                                                                                                                                           >  2:209287/1‑247 (MQ=255)

                                         gTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACAGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACg                                                                                                                                                                                                       <  1:530673/250‑1 (MQ=255)

                                             tggcggggcacaTGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCGATATAGCGGGAGTAGACAACCTTGCCCGGCGCGTGGAACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGGGTTCGGGTTTTCGGCAAACTCAGGGCGCAGCGCAAGCAGGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGAACTGTCA‑ACGcgca                                                                                                                                                                                                   >  2:41344/9‑247 (MQ=255)

                                                cGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATccc                                                                                                                                                                                               >  1:887059/1‑251 (MQ=255)

                                                       aaTGCAGGATTTTTGCCCAGAACCATTCGCTGGAATCAATACCGCCAATCTAGCGGGAATAGTCAACATTGCCCGGCGCGAGGCACAAACGGGTCATCTCT‑TC‑GCATCTTTAAACCGCAGTGT‑GGTCTTTCCCCAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTGCGTGCG‑TCGAGCCGGAACAGTCA‑ACTCCAATCCCGCCcac                                                                                                                                                                                         >  2:849121/2‑250 (MQ=255)

                                                         tGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCAGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGc                                                                                                                                                                                      >  1:957688/1‑251 (MQ=255)

                                                            aGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTgc                                                                                                                                                                                   >  1:787385/1‑251 (MQ=255)

                                                                 ttttGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCt                                                                                                                                                                              >  1:647870/1‑251 (MQ=255)

                                                                     tcccaaaacctgtcgcggaaaaaaataCCGCCAATATAGCGGGAGTAGGCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGTCG‑TCGAGCCGGTACTGGCA‑AAGCGAATCCCGACCACAGCTGCGCGCTGtt                                                                                                                                                                           <  1:803287/228‑1 (MQ=255)

                                                                           aaCCATTCCCTGGACTAAAAACCGCCAATATAGCGGGAGTAGTCAACATTCACCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTACACAAT‑CCGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGCCGCAGCGCCAGCCCGGTACCGTCGGCATCCACAGGCGAGGGCG‑TAGCGCCGGTACTGTAC‑ACACCAATCCCGCCcaacgatgcgcgcagtacgaagcc                                                                                                                                                                    >  2:458370/1‑229 (MQ=255)

                                                                             ccttctgttgaaaaaaataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTaa                                                                                                                                                                  <  1:978456/237‑1 (MQ=255)

                                                                             ccATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTa                                                                                                                                                                   >  2:756455/1‑250 (MQ=255)

                                                                               aTTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCAAGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGCTCGACGCTAAGc                                                                                                                                                                >  1:745034/1‑251 (MQ=255)

                                                                                            aatcccgccaaaatatcGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGTTAAT‑TCC‑TTCGCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCTCAGCGCCAGCACGTTACCGTCGGCATCAATCGTTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGTCGCTAAGCTCTGCAAGCACgg                                                                                                                                                   <  2:754382/235‑1 (MQ=255)

                                                                                            aatcccgccaATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACgg                                                                                                                                                   <  1:347679/244‑1 (MQ=255)

                                                                                            aataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACgg                                                                                                                                                   <  1:668023/251‑1 (MQ=255)

                                                                                               accgccaATATAGCGGGCGGAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATAGCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCCCAAG‑ACGAACATCGGTTTCGGATTTTCGGAAAACTCGGTGCGCAGCCCAAGCACGGTACAGTCGCCTTCAATCGACGCGG‑CGATCGAGCCGGAACTG‑CACACGCCAGTCCCTACCACATCGGCGCCCGGTTCGCGGCTAAGCTCAGCAAACACGGttt                                                                                                                                                >  2:587592/5‑251 (MQ=255)

                                                                                                ccgccaATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGtttt                                                                                                                                               >  1:691152/4‑251 (MQ=255)

                                                                                                 cgccaATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTc                                                                                                                                              >  1:270450/3‑251 (MQ=255)

                                                                                                     aagagagcgggggtatccAACATTGCCCGGCGCGGTGCCCAAACGGGTAATCTCT‑TC‑GCGTCTTTCAACCGCAGGGT‑GGTCTTTCCACAAT‑CCGAACACCGCGTGCTGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCCCGTTACCGTCGGCATAAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGTCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCGCTGCAAGCACGGTTTTCAGTg                                                                                                                                          <  2:975851/234‑1 (MQ=255)

                                                                                                      atatAGCGGGAGTAGGCAACATTGCCCGGCGCGTGGCAAAAACGGGTAATCTCT‑TC‑GTTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGc                                                                                                                                         <  1:552/248‑1 (MQ=255)

                                                                                                       tatAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGcc                                                                                                                                        >  1:30865/3‑251 (MQ=255)

                                                                                                        atAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCg                                                                                                                                       <  1:728413/250‑1 (MQ=255)

                                                                                                        atAGCGGGAGTAGACAACATTGCCCGTCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACGCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGcc                                                                                                                                        <  2:667874/249‑1 (MQ=255)

                                                                                                        atAGCGGGAGTAGACAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑CCACCAATCCCGACCACCGCTGCGCGCTGTTCGACGCTAAGCTCTGCTAACACGGTTTTCAGTGccc                                                                                                                                       >  2:931220/2‑250 (MQ=255)

                                                                                                        agagggggAGTAGTCAAAATTGGCCGGGGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCGTTCCACAAT‑ACGAACCTCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCGGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCg                                                                                                                                       <  2:760786/246‑1 (MQ=255)

                                                                                                           gCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCtt                                                                                                                                    >  1:845617/1‑251 (MQ=255)

                                                                                                             gggAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGAAATCTCT‑TC‑GCTTCTTTGAACCGCAGTGT‑GGTCTATCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGATACGGGCG‑TCGAGCCGGAACTGTCA‑ACGACAATCACGACCACAGCTGCGCGCTATTCGACGCTAAGCTCTGCCAGCACGGTTTTCGGTGCCGCttac                                                                                                                                  >  1:919566/1‑249 (MQ=255)

                                                                                                             gggAGTAGTCAACATTGCCAGGAGCGGGGACCAAACGGGTGATCTCT‑TA‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGGTACCGTCGGCATCAATCGGTACGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAGTCCCGCCCACAGCAGCGAGCTGTTCGACGCTAAGCTCTGCCAACACGGTTTTCAGAGCCGCTTcc                                                                                                                                  >  2:181629/1‑251 (MQ=255)

                                                                                                              ggAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGCACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGAATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCCAGGTGTTCAGTGCCGCTTccc                                                                                                                                 >  2:814803/1‑250 (MQ=255)

                                                                                                               gAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCGGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCAt                                                                                                                                >  2:265623/1‑251 (MQ=255)

                                                                                                                 gtagtCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTg                                                                                                                              <  1:771548/251‑1 (MQ=255)

                                                                                                                    gtCAACATTGACCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGATTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑GCGAGCCGGTACTGTCA‑CCGCCAATCCCGACCACAGCTGCGCGCGGTTCGCCGCTAAGCTCTGCAAGCAAGGCTTTCAGTGCCGCTTCAATTGAc                                                                                                                            >  2:780488/1‑250 (MQ=255)

                                                                                                                     tCAACATTGACCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGACCATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑CCGCCAATCCCGACCACAGATGCGCGCAGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCCATGCCGCTTCCCTTGAccc                                                                                                                          >  2:645120/1‑249 (MQ=255)

                                                                                                                     tCAACATTGACCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGCCCCCAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCACATCCATTGtctc                                                                                                                          >  2:891180/1‑251 (MQ=255)

                                                                                                                         cattGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGCGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATg                                                                                                                      >  2:285799/3‑251 (MQ=255)

                                                                                                                          attGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGt                                                                                                                     <  1:397100/250‑1 (MQ=255)

                                                                                                                                cGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCAc                                                                                                               <  1:422019/251‑1 (MQ=255)

                                                                                                                                cGGCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCAc                                                                                                               <  1:648851/251‑1 (MQ=255)

                                                                                                                                 ggCGCGTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCAcg                                                                                                              >  1:788428/1‑251 (MQ=255)

                                                                                                                                  cagcgTGGCACAAACGGGTAAGCTCTGTGCG‑TCCTT‑CAACCGCAGGGT‑GGTCTTTCCAAAATCAC‑ACCATCGCGTGCGTGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAAGCGGTGCGGGCG‑CCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGATTCCATTGCCTCAATGTAGTCAcgc                                                                                                             <  2:908698/249‑1 (MQ=255)

                                                                                                                                   cgcgTGGTACAAACGGGTAATCTAT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTAACGTCGGCATCAATCGTTGCGGGCG‑TCGCGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCCCGCCAGGTTTTAAGTGCAACTTCCATTGACTCAATGTAGTCAcgcg                                                                                                            >  2:984034/1‑251 (MQ=255)

                                                                                                                                   caggtgtcataaaCGGGTAACCTCT‑TC‑ACTTCTTTCAACCGCAGTGTTGGT‑TGTCCACAAT‑ACGAACTCCGCGTGCGGGTTTTCGGCAAGCTCCGGGCGCAGCGCCAGCACGTTCCCGTCGGCATCAATAGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACACAATGTAGTCAcgcg                                                                                                            <  2:8073/241‑1 (MQ=255)

                                                                                                                                     cgTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCCTCCATTGACTCAATGTAGTCACGCgg                                                                                                           >  2:667881/1‑250 (MQ=255)

                                                                                                                                      gTGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGgat                                                                                                         >  2:876233/1‑251 (MQ=255)

                                                                                                                                      gTGGCACAAACGGGTAATCTCT‑TC‑ACTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCTCGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGgat                                                                                                         <  1:235364/251‑1 (MQ=255)

                                                                                                                                       tGGCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGgatg                                                                                                        >  1:490262/1‑251 (MQ=255)

                                                                                                                                       aggCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCCATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCAGCAAGCACGGTTTTCATTGCCGCTTCCATTGACTCAATGTAGTAACGCggttg                                                                                                        >  2:220480/2‑248 (MQ=255)

                                                                                                                                        ggCACAAACGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGgatga                                                                                                       <  1:399096/251‑1 (MQ=255)

                                                                                                                                               aCGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATGGACTCAATGTAGTCACGCGGATGATGACGGa                                                                                                <  1:908508/251‑1 (MQ=255)

                                                                                                                                                cGGGTAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGaa                                                                                               <  1:459766/251‑1 (MQ=255)

                                                                                                                                                    tAATCTCT‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGAAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTgg                                                                                           >  2:409896/1‑251 (MQ=255)

                                                                                                                                                        ctct‑TC‑GCTTCTTTCAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTAt                                                                                       >  1:379193/1‑251 (MQ=255)

                                                                                                                                                              tcGCTTCTTGCAACCGCAGTGT‑GGTCTTTCCCCAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTc                                                                                     <  2:295173/248‑1 (MQ=255)

                                                                                                                                                                  tcctt‑cAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTAAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGc                                                                                <  1:598789/247‑1 (MQ=255)

                                                                                                                                                                        cAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATcaca                                                                          <  1:756455/250‑1 (MQ=255)

                                                                                                                                                                          aCCGAAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGAGGATGATGACGGAACTGGTTATTAGGGGCATCACaaa                                                                        >  1:784505/1‑250 (MQ=255)

                                                                                                                                                                                gttt‑gGTCTTTCCCCACT‑ACGAACAGCGCGTTCTGGTTTCCGACAAACCCCGGGCGCAGCGCCAGCACGTTACCTCCGGCATCAACCGGGGCTGGCG‑TCGAGCCGTTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTTTCCGACGCTAAGCTCTGCAAGCACGGTTTGCAGTGCGGCTTCCATGGACTCAATGTAGTCACGCGGATGATGACCGAACTGGTTATTCGGTGCATCACaaaactgtcc                                                                 <  2:641292/248‑7 (MQ=255)

                                                                                                                                                                                  gt‑gGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCt                                                                <  1:93996/250‑1 (MQ=255)

                                                                                                                                                                                     ggtcttctcccaAT‑ACGATCATCGCGTTCGGGTTGTCGGCAAACTCCGGGCGCCGCGCCATTCTGTTACCGTCGGCAGCAATCGGTGTGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGGCCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACGGTTTATTCGGGGCATCACAAAATTGCCCtttt                                                             <  2:95661/241‑2 (MQ=255)

                                                                                                                                                                                     ggTCCTTCCGCAAT‑ACGACCATCGCGTTCGGGTTTTCGGCAAACTCCGGCCGCAGCGCCAGCACGTTACCGTCGTCATCAATCTTTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTGTCAGTGCCGCTTCCATTGACTCAATGTAGGCACTCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGcccgttt                                                             <  1:849121/251‑5 (MQ=255)

                                                                                                                                                                                      cccTTTCCCCAAT‑ACGAACATCGCGTTCGGGTTTCAGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGGCGGCATCAAGCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCTAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTACAGGGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCttttt                                                            <  2:342718/249‑3 (MQ=255)

                                                                                                                                                                                             cacaAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCGGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACg                                                     <  1:953606/251‑1 (MQ=255)

                                                                                                                                                                                                 aT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGa                                                 <  1:99883/251‑1 (MQ=255)

                                                                                                                                                                                                     cGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGACCTGGTTATTCGGGGCATCACACAATTGCCCTTTTTGCCAACGGGGATAc                                              >  1:738659/1‑251 (MQ=255)

                                                                                                                                                                                                           ccgcgTGCGGGTTTCCGGCAAACTCCGGGCGCAGCGCCAGCACGTACCCGGCGGCATCAACCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCActct                                        <  2:213160/250‑1 (MQ=255)

                                                                                                                                                                                                             gcgTTCGAGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTAGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTAc                                      >  1:924787/1‑251 (MQ=255)

                                                                                                                                                                                                              cgTTTGGGTTTGCGGCAAACACCGGGCGCCGCGCCAGCACGTTACCGTCGCCATCAAGCGGTGCGGGAG‑TCGAGCCGGTTCTGTCA‑ACGCCAATCCCGACCACCGCTGCGCGCTGTTCGGCGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTGTTGCCAACTGTGATACCACGCTACg                                     <  2:616301/251‑1 (MQ=255)

                                                                                                                                                                                                                ttCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCt                                   >  2:48662/1‑251 (MQ=255)

                                                                                                                                                                                                                ttCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCt                                   >  1:28553/1‑251 (MQ=255)

                                                                                                                                                                                                                  cGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCACGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGAAACCACTCTACGCtgg                                 >  2:312944/1‑251 (MQ=255)

                                                                                                                                                                                                                      ttttCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGc                             >  1:831218/1‑251 (MQ=255)

                                                                                                                                                                                                                       tttCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCg                            >  1:666888/1‑251 (MQ=255)

                                                                                                                                                                                                                             cAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATctct                      >  1:404450/1‑251 (MQ=255)

                                                                                                                                                                                                                               aaCTTCGGGCGCAGCGCCAGCACGTTACCTCCGGCATAAATGGGTGTGGGCG‑GCGAGCCGGTACTGTCA‑ACGCCAACCCCGACCACAGCGGCGCGCTGTTCGACGCTAAGCGCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACGGGTTATTCGGGGCATCACAAAATTGCCCTGTTTGCCAACGGGGATTCCACTCTTCGCTGGGTGCGATCTCTTc                    <  2:56239/251‑1 (MQ=255)

                                                                                                                                                                                                                                accccGGGCGCGGGGCCAGCACGTTACCGTCGGAATCAATGGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCTCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATGGCCCTTTTTGCCAACGGGGATTCCACTCTACGCTGGTGGCGATCTCTTCa                   <  2:887059/248‑1 (MQ=255)

                                                                                                                                                                                                                                             cgcCAGCACGTTCCCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGgat                                                                                                         <  2:363651/152‑1 (MQ=255)

                                                                                                                                                                                                                                             cgcCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGgat                                                                                                         >  1:363651/1‑152 (MQ=255)

                                                                                                                                                                                                                                              gcCAGCACGTTACCGGCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTc     <  1:763796/251‑1 (MQ=255)

                                                                                                                                                                                                                                                 aGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCAc  >  1:879763/1‑251 (MQ=255)

CACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑TCCGCTTCTT‑CAACCGCAGTGT‑GGTCTTTCCACAAT‑ACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCCAC  >  NC_000913/69507‑69994

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