BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,135,444	G→T	11.2%	A309E (GCA→GAA)	pitB ←	phosphate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,135,444	0	G	T	11.2%	344.5 / 2.8	132	A309E (GCA→GAA)	pitB	phosphate transporter
Reads supporting (aligned to +/- strand): ref base G (34/80); new base T (15/3); total (49/83)
Fisher's exact test for biased strand distribution p-value = 2.28e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 134 total aligned reads displayed.

CATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCC  >  NC_000913/3135206‑3135682

cATGGTGATCCCCCCCGTGGGGTTCGCAATGGTGCTTAACGTTGCGCTGCGAAGTTTTTTCAGCAGGTTCTGGGCTGCTTGAATGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGGATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑Tg                                                                                                                                                                                                                                                      <  1:1253581/239‑1 (MQ=255)

                ggtgcgttttcatttgtccTTTACATTGCGCTGCGAAGTTGTTGCAGCAAGTTGTGGTCTTTTTTACTGACGCCTGTGAGTTTGGCTAGCTTCGCGGAGGGATCAGTGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCAGATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGTATGAC‑A‑GT‑Gaaa                                                                                                                                                                                                                                   <  1:900571/221‑1 (MQ=255)

                 tggggttcTCGATGGTGTTTAACATAGTGCTGGGGAGTTTTTTCAGCAGGTTCTGGTCTTCGTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCACCAAAGGTATT‑TGCCGGATGAC‑C‑GG‑GAAAc                                                                                                                                                                                                                                  <  1:313202/232‑1 (MQ=255)

                  aGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTGACTTTACATATTGCTTGGCAGC‑ACCGTTTTTACGCGCCCAATAAACTCCAAGGGATT‑TTACGGAtttacctgaaa                                                                                                                                                                                                                                      >  2:268262/1‑227 (MQ=255)

                  aGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAAc                                                                                                                                                                                                                                  <  2:866367/238‑1 (MQ=255)

                   gCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCAATTGCCCGGCAGA‑ATCGCTTTAACTCGCGCAATATCATCAAAGGGATT‑TTCCCGGTGAC‑A‑GT‑Gaaaa                                                                                                                                                                                                                                  >  2:222813/1‑236 (MQ=255)

                       aCTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCTATTGCCTGGCAGC‑ACCGTTTTAACCCGCGCCAAACCATCAAAAGGATT‑TTCCGGGTGAccggggaaacccttt                                                                                                                                                                                                                               >  2:1212493/1‑225 (MQ=255)

                       aCTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAAATTTCCAAATTTCCTGGCAACAACCGTTTTAA‑GCACGCAAATGCATCAAAGGGTTT‑TTCCGGATGACAA‑GT‑Gaacccgatt                                                                                                                                                                                                                             >  2:827941/1‑232 (MQ=255)

                            aTGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGGAACTTTCCATATTTCCTGGCAGC‑ATCGTTTTAACGCCCGCAATAGCATCAAAGGGATT‑TTCCCGATGAC‑A‑GT‑GAAACTCTCTTAC‑tt                                                                                                                                                                                                                       >  2:1058622/1‑238 (MQ=255)

                             tGGTGCTTAACGTTTCCCGGGGAAGTTTTTTCAGCAGGTTCTGGGCTTCGTTTCTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGc                                                                                                                                                                                                                     <  1:904456/239‑1 (MQ=255)

                             tGGTGCTTAACATTTCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGc                                                                                                                                                                                                                     <  2:238818/239‑1 (MQ=255)

                                     aaCATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGGTTTAACGCGCGCAATAGCATCAAAGGGATT‑TGCCGGATTAC‑A‑GT‑GAAACTCTTTTAC‑TTGGCTGccct                                                                                                                                                                                                              >  2:904456/1‑236 (MQ=255)

                                        ttttcgcgtggAAGTTTTTTCAGCAGGTTCTGGGCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGt                                                                                                                                                                                                          <  1:222813/229‑1 (MQ=255)

                                        ctttcgtttgcaaGTTTTTTTAGCGGGTTGTGGTCTTGTTTCCTGACGCCTGGCATTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGt                                                                                                                                                                                                          <  1:1316717/228‑1 (MQ=255)

                                          atCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGCCGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCACATTACTTTCCCTATAGCTTGGCAGC‑ATCGGTTTAACGCGCGACATAACTTCAAAGGAATTGCTCCGTAT‑AC‑A‑GT‑TAAACTCTTTTAC‑TTgtatttcttacatac                                                                                                                                                                                                        >  1:1040001/1‑225 (MQ=255)

                                          atCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCTATTGCCTGGCAGC‑ATCGGTTTAAAGCGCGCAAAAGCATCAAAGGGATT‑TTCCGGGTGGA‑A‑GG‑GAAAACCTGTTTA‑TTGGGTGCCACCAg                                                                                                                                                                                                           >  2:863237/1‑236 (MQ=255)

                                            ccctgggcttttttttCGGCGGGTTCTGGGTTTTTTTAATGACGCCTGGCAGTTTCGGTGGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGAt                                                                                                                                                                                                      <  1:1212493/229‑1 (MQ=255)

                                               tGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGAtgct                                                                                                                                                                                                   <  2:1109500/239‑1 (MQ=255)

                                                 cGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCCGC‑ATCGTTTTAACGCGGGCAATAGCATCAAAGGTATT‑TGCCCGATGAC‑A‑GG‑GGAACTCTGTTAC‑TTTCTGGCCATAAGTCGAagctgc                                                                                                                                                                                                 >  2:1169957/1‑239 (MQ=255)

                                                 cGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCCTAAATTTTCCTATTGACCGGCAGC‑ATCGCTTTAAACCGCGCAATAACATCAAAAGGATT‑TTCCGTGAGAC‑A‑GT‑GAAACTCTGTGAC‑tttctggcgctaagcgttggcgtc                                                                                                                                                                                                 >  2:1070764/1‑217 (MQ=255)

                                                     gTTTTTTCAGCCGGTTTTGGTCTTCTTTACTGACGCCTGGGCGTTTCGGCAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGc                                                                                                                                                                                            <  1:1193045/239‑1 (MQ=255)

                                                                             cttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTTCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGGATT‑TGCCGGATGAC‑A‑AT‑GAAAAACTGTTAC‑TTGCGTGGCATCAGTCGATGCTGC‑AGGCAATTGGGGGTCCCTCCCaaacc                                                                                                                                                                      >  2:1174015/1‑234 (MQ=255)

                                                                                tACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑CTCGTTTTTAAACGCGCAAAAGCATCCAAGGACTT‑TGCCCGATG‑C‑A‑CTGGAAAACCTGTTACACTGTGT‑CCAACAGCCATTGCCGC‑CGGCAATGGCGGTGCAACCGCAaacaac                                                                                                                                                                    >  1:125648/1‑236 (MQ=255)

                                                                                    gccGCCTGGCAGTTTCGGTAGCGTCGGGGAGGTAATAGAGATGTCCAGCAAAATGCGGCGCAGCTGGCTGCGCTTAATCACGCTTAACGGCTCGTAACTTTCCAAATTGACTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCAGCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACGCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCg                                                                                                                                                             <  1:1284486/237‑1 (MQ=255)

                                                                                     cggcctgtgattttGGCTACGTTCGCGGGGTTGTCAGAGATGCCCAGCATAATGCGGCGCAGCTGGCTGCGCTGATTCACGCTTAACGGCCCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCgg                                                                                                                                                            <  1:1026440/228‑1 (MQ=255)

                                                                                     aCGCCTGGCAGTTTCGCTGCGTTCGCGGGGGTACCAGAAATGCACAGCATAATGCGGGGCAGCTGGCTGCGCTGACTGACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCgg                                                                                                                                                            <  1:863237/239‑1 (MQ=255)

                                                                                     aCGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCCTCAGTCGATGCTGC‑AGGCAATGGAggctgtcgcttttacaccactccga                                                                                                                                                             >  2:1097391/1‑215 (MQ=255)

                                                                                      cGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCgg                                                                                                                                                            <  2:188071/238‑1 (MQ=255)

                                                                                                 ttCGCTAGGTGCGCGGAGGTATCAGAGATTCACAGGATAATGCGGCGCAGCTGGCTGCGCTGACTCCCGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATg                                                                                                                                               <  2:1195296/239‑1 (MQ=255)

                                                                                                 ttCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGa                                                                                                                                                 <  2:452717/237‑1 (MQ=255)

                                                                                                 ttCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCAAATTTCCTGGCACC‑CTCGGTTTTACGCCCGCCATAGCCACAAAGGGATC‑TGTCCGTAGAC‑A‑AT‑GAACACCCGTGAC‑CTTAGTGCCAACAGTCGA‑GTTCCGTGGCCACGGAGGTTCTATCGCAATCacaacttcgtcccgatccggcgg                                                                                                                                                >  1:848044/1‑217 (MQ=255)

                                                                                                   cGCTAGCTTCGCGGAGGTATCAGGGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCt                                                                                                                                             <  2:649040/239‑1 (MQ=255)

                                                                                                            cggggAGGTTTCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGTTTGCTCCCGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGt                                                                                                                                    <  1:65907/236‑1 (MQ=255)

                                                                                                            cccgGGGGTTTCAGAGATGCCCATCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGt                                                                                                                                    <  1:586837/237‑1 (MQ=255)

                                                                                                             gaggAGGTATCAGAGATGCAAAGCATAATGCGGCGCAGCTGGCTGCGCTTACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTa                                                                                                                                   <  2:1280345/237‑1 (MQ=255)

                                                                                                                   gtttCAGAGATGCACGGCATAATGCGGCTCAGTTGGCTGGGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGtt                                                                                                                              <  1:866472/235‑1 (MQ=255)

                                                                                                                   gTATCGGAGTTGCACCGCATAATGCGGCGCGGCTGGCGTCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCCGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGtt                                                                                                                              <  1:107312/238‑1 (MQ=255)

                                                                                                                   gTATCAGAGATGCACAGCATAAGGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGtt                                                                                                                              <  2:420489/238‑1 (MQ=255)

                                                                                                                     aTCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGGGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCg                                                                                                                            <  2:66192/238‑1 (MQ=255)

                                                                                                                             ttcacaGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTCGGTa                                                                                                                   <  2:1145052/237‑1 (MQ=255)

                                                                                                                              gCACAGCATAATGCGGCGCAGGTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCg                                                                                                                                                                                                        <  1:197762/155‑1 (MQ=255)

                                                                                                                              gCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TTCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCg                                                                                                                                                                                                        >  2:197762/1‑155 (MQ=255)

                                                                                                                                      tatgggggggcgggtgggctggcTGACTTACGTTTAACGGTTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATcg                                                                                                          <  1:582318/218‑1 (MQ=255)

                                                                                                                                         tGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGt                                                                                                                                                                                                                   <  2:513677/133‑1 (MQ=255)

                                                                                                                                             gcgcAGCTGGCTGCGGTGACTCACGCTTAACGGGTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggta                                                                                                    <  2:1040001/238‑1 (MQ=255)

                                                                                                                                             gcccAGCTGGCTGGGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggta                                                                                                    <  2:125648/235‑1 (MQ=255)

                                                                                                                                              cgcAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑CTCGTTTTTAACCGCGCAATAACCTCAAAAGTATT‑TTCCGGGTGAC‑C‑GT‑GAAACTCTGTTAC‑TTTCGGGCCCTCAACTGATGCTGC‑CGGCAATGGGGGGTCCATCTCCATCACCTTCTGCGCCAGTTCC‑GGATGCCTTTGCGGGTAGTGTTTTCAGTTGGTAACGGCACCGcgggggag                                                                                                  >  1:240223/1‑235 (MQ=255)

                                                                                                                                                 aactggctgCCCTGTCTCACGCTTAACGGTTCGTAACTTTCCATAGTGCCAGGCAGC‑ATCGTTTTAACGCACGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑TAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggtacgggt                                                                                               <  2:783067/237‑1 (MQ=255)

                                                                                                                                                               aCTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCAGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCgg                                                                                 <  2:1165225/239‑1 (MQ=255)

                                                                                                                                                                cccccgtttACCGGGTTGTAACTTTCCATATTTCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGa                                                                                <  1:1174015/232‑1 (MQ=255)

                                                                                                                                                                cTCCCGTTTAACGGCTGGTAACTTTCCATATTGCCGGGCAGC‑ATGGTTGTAACGCGCGCAATAGCATCAAAGGGATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGGTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AAGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGa                                                                                <  2:1256001/239‑1 (MQ=255)

                                                                                                                                                                cTCAAGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGa                                                                                <  2:279666/239‑1 (MQ=255)

                                                                                                                                                                  cccGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACg                                                                              <  2:293187/237‑1 (MQ=255)

                                                                                                                                                                    ctgttAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATTGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCa                                                                            <  2:927810/236‑1 (MQ=255)

                                                                                                                                                                    cGCTTACCGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCa                                                                            <  2:1280278/239‑1 (MQ=255)

                                                                                                                                                                    cGCTTAACGGCGCGTAAATTTCCCTATTGCCTCGCAGC‑CTCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TTCCCGATGAC‑A‑GT‑GAAAATACGTTAC‑TTGCGTTCCATCAGTCGATGCTTC‑AGGCAATGGAGGGTCCATCGCGAGCAACTACTGCGCCAGTTCACGGAT‑AGGTTGCACAGACTGGTCGAAGTTGGGAACCGCATCCcggggacgggtgaatataagaacggagcac                                                                            >  1:66880/1‑220 (MQ=255)

                                                                                                                                                                         aaCGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTTCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGGTCCATCGCAATCAACTTCTGCGGCCGTTCC‑GGAGGCTGTTTCAGGTGTGGTTCGAAGTTGGTAAAGGCTTCCCGGGGACGGGGAAATTctatgcccgggccgtta                                                                         >  2:125855/1‑221 (MQ=255)

                                                                                                                                                                         aaCGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATTGCAATCAACTTTTTCGGCAGTTCA‑GGATGCTGTTTCAGGTGGGGTTCGAGGTTGGGAACAGCATCCCGGGGAAGGGGAATTTTTAAGCCggggcgagtcac                                                                       >  2:537691/1‑229 (MQ=255)

                                                                                                                                                                            ggCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATAtt                                                                    <  2:553307/239‑1 (MQ=255)

                                                                                                                                                                                     cTTACCATATTGCCTGGCAGC‑AACGTTTTAACGTGCGCAATGGCAGCAAAGGTAGT‑TGCCGGATCCC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTCGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGaa                                                           <  2:478325/239‑1 (MQ=255)

                                                                                                                                                                                      tttCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGGAGTGTTCGCAGTTGGTAACGGCATCGCCGGTACGGGGAATTTCATAGCCCGACCCCTTCATTTTGACGAGGaa                                                           >  2:1198102/1‑238 (MQ=255)

                                                                                                                                                                                        tCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGCTGCTGC‑AGGCAATGGAGGTTCCACCGCAATCCACTTCTGCGGCGGTTAC‑GGATGCTGTTGCAGGTGTGGTTTCCACTTGGTAACGTCCTCCCGGGGCCGGGCACTTTCAGAGCCGGAGGCCTTACTAttttggcacaacg                                                         >  1:30293/1‑227 (MQ=255)

                                                                                                                                                                                            tatTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAACCAACTTCTGCGGCAGTTCC‑GGATGCTTTTGCAGGTAGTGTTCGAAGTTGGCAACGGCCTCGGGGGTACGGGGAATTTCAAAACCGGGCGCTTTCACAtttaggcagaaaacagcg                                                    >  1:710521/1‑223 (MQ=255)

                                                                                                                                                                                             atTGCCTGGCAGC‑ATCGTGTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGAGCAATCAGTCGATGCTGC‑AGGCATTGGAGGTTCCATCGCAATCAACTTCTGCGTCAGTCCA‑GGCTGCTGTTGCAGGTAGAGTTCGCACTTTGGAACCGCCCCCCCGGTAAGGGCAATTTCATTGCCGGATTCATTTAAAttaaaagagcagccgaaag                                                   >  1:1016002/1‑222 (MQ=255)

                                                                                                                                                                                               tGCCTGGCAGC‑ATCGTTTTAATGCGCGCAATAGCATCAAAGGTATT‑TTCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATAAGTCTGTGCTGC‑ACACACTGGTGGCTCCCTCGCAATCAACTTCTGCGTCAGTTCA‑GGCTGCCGCTGCAGGTAGTGATGTCAGTAGGTACGTGCCTCCCGGGCACGGGCAATTTTTTAGTCGGACGTCTTCACAATTACGACAAAGcccgacgg                                                  >  1:441791/1‑232 (MQ=255)

                                                                                                                                                                                               tGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCa                                                                            <  2:1092142/212‑1 (MQ=255)

                                                                                                                                                                                                gCCTGGCATG‑ATCGTTTTAACGCGCGCAATAGCCACAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAgggg                                                <  2:1294757/239‑1 (MQ=255)

                                                                                                                                                                                                gCCTGGCAGC‑ATCGTTTTAACGCGCGGAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAgggg                                                <  2:419652/239‑1 (MQ=255)

                                                                                                                                                                                                gCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCACCTTCTTCGGCAGTTCA‑GGAGGCTGTTGCTGGTAGGGTTTGAAATTGGTAACGGCCTCTCGGGGACGGGGAATTTCATAGCCGCGGGCATTTATATTGTCGACAAAAACAACAgggg                                                >  1:65027/1‑239 (MQ=255)

                                                                                                                                                                                                 ccTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGCCAGTTCA‑GGATGCGGTTGCCGGTAGTGTTCGAAGTTGGTAACGGCGTCCGGGGTACGGGGAATTTCAATGCCGGGGGCATTTAAAtttccgcgcaaaccaccaggggg                                               >  1:1263145/1‑218 (MQ=255)

                                                                                                                                                                                                    ggCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggta                                                                                                    >  2:903238/1‑183 (MQ=255)

                                                                                                                                                                                                    ggCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCaa                                             >  2:1072245/1‑238 (MQ=255)

                                                                                                                                                                                                      cagc‑aTCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGCGGCGTTCAAATTGTACACGAAGCCAGCAGGGGCaaac                                           >  2:507660/1‑236 (MQ=255)

                                                                                                                                                                                                         c‑aTCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGCTGCGGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGGACGGGTAATTTCATAACCGGAGGCATTCATATTTACGCCGAAGCCAGCAGGGGCAATcccc                                        >  2:1041172/1‑238 (MQ=255)

                                                                                                                                                                                                           tgttttttAACGCGGTCAATAGGATCAAAGGTATT‑TGCCGGATGAA‑C‑GT‑GAAACTCTGTTAC‑TTGAGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCTACGCAATCAACGTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCccac                                      <  1:1080283/235‑1 (MQ=255)

                                                                                                                                                                                                           aTCTTTTTAACCCGGGGAATAGCTACAAAGGTTGT‑TGCCGGCTTAC‑A‑GG‑TAAACTCTGTTAC‑GTGCGTGCAATTGGTCGGTGCTTC‑AGGCAAGGGAGGTCCTATCCCGATCAACTTTGGCAGCAGTTCA‑GGAGGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTTCGGGTAATTTCATAGCCGGACGGATTCATATTGACGACGAAGCCAGCAGGGGCAATCCcca                                       <  1:869651/238‑1 (MQ=255)

                                                                                                                                                                                                           aTCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GAATGCTGTTGCAGGTAGAGTTCGAAGTTGGTAACGGCATCGCGGGAACGGGACATTTCATAGCCGGACGAATTCATATTGAGACCGAAGCCAGCAGGGGCAATcccacc                                      >  2:644455/1‑236 (MQ=255)

                                                                                                                                                                                                            gcGTTTTAACGCGCGCAATAGAGTCAAAGGTATT‑TGCCGGATAAC‑A‑GT‑GAAACTCTGTTGC‑TTGCATGCCAGCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGGTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCccac                                      <  2:505003/237‑1 (MQ=255)

                                                                                                                                                                                                               ttttAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAgggg                                                <  2:1072646/225‑1 (MQ=255)

                                                                                                                                                                                                               ttttAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGAAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGCTGCTGTTGCCGGACAGGTTCGAAGTTGGGCACGGCCTCGCGGGGCCGGGGACATTCCTAGCCCGGCCCATTCTAATTGACGGCTAAGCCAACAGGGGGaaaaccccacaaa                                  >  1:71948/1‑228 (MQ=255)

                                                                                                                                                                                                                tttAACGCGCGCAATAGCATCAAAGGTATT‑GGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑GTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTa                                 <  2:71584/239‑1 (MQ=255)

                                                                                                                                                                                                                tttAAAGCTCGCCATTGCATCAAAGGGATT‑TGCCGGTTGAC‑A‑GG‑TAAACTTTGTTAC‑GTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTa                                 <  2:748520/239‑1 (MQ=255)

                                                                                                                                                                                                                   accgcgcg‑AATAGCAACAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTTTTTTAC‑TTGCGTGCCATCAGTCGTGGCTGG‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAg                             <  1:907969/237‑1 (MQ=255)

                                                                                                                                                                                                                    acgggcggcaTAGCATCAAAGGTTTT‑TGCCGGAATAC‑A‑GT‑GAAACTTTGTTAC‑TTGCGTGCCATCAGTCCATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGTTAGTGTTCGAAGTTGGTAACGGCATCGCTGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAg                             <  2:828450/230‑1 (MQ=255)

                                                                                                                                                                                                                     ccccccgaaTAGCTTCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑AAAACTCTGTTAC‑TTGCGTGCCATCATTCGAGGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGc                            <  1:1175216/232‑1 (MQ=255)

                                                                                                                                                                                                                         cgcAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGCCGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCAtt                         >  2:227610/1‑238 (MQ=255)

                                                                                                                                                                                                                         cgcAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCCtt                         >  2:918236/1‑238 (MQ=255)

                                                                                                                                                                                                                         cgcAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGGACGGGGAATTTCCTAGCCGGGCCCATTCATATTTACGACGAAGCCGGCAGGGGAAAAcccccc                                      >  1:1072646/1‑223 (MQ=255)

                                                                                                                                                                                                                            aaTAGAATCAAAGGTATT‑TTCCGGATGAC‑A‑GT‑GAAACTTTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCGGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTAcc                      <  2:702877/238‑1 (MQ=255)

                                                                                                                                                                                                                                  tccAAGGGTTTT‑TG‑CGGATGAA‑AGGT‑GAAACTCT‑TTAC‑TTTGGGTCCATTAGTCGAGGCGTC‑AGGCAATGGGGGTTCCATCGCAATCAACTTGTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTTCCAGCATTACCAGGCCGAt              <  1:933911/237‑1 (MQ=255)

                                                                                                                                                                                                                                  aTCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCTCCACCAGTACCAGCATTACCAGGCCg                <  2:535461/238‑1 (MQ=255)

                                                                                                                                                                                                                                   tCAAATGTATT‑TGCCGGATGAC‑A‑GT‑TAAACTCGGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGAt              <  2:1144324/239‑1 (MQ=255)

                                                                                                                                                                                                                                   tCAAAGGTATT‑TGCCGGATGAC‑A‑GG‑GAAACTCTGTTAC‑GTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGGTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGAt              <  2:848044/239‑1 (MQ=255)

                                                                                                                                                                                                                                    cAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTACTTTCCTAGCAGGACGCATCCATTTTTACACCGAAATCGCCAGGGGCAATCCCCACCAGACCCCCCATTCCCAGGCCgca              >  1:802490/1‑236 (MQ=255)

                                                                                                                                                                                                                                      aaGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGCTGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGGAACGGAATCTCGGGTCCGGGTAATTTCATAGCCGCACGCCTTCATATTGACGACGAACCCAGCAGGGGCAAACCCCCCCAGTAACACCATTACCACGCCGATc             >  1:282293/1‑237 (MQ=255)

                                                                                                                                                                                                                                      aaGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCACTACCAGCATTACCAGGCCGATc             >  2:867545/1‑237 (MQ=255)

                                                                                                                                                                                                                                           aTT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGCAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATACCtcttt      >  2:200938/1‑239 (MQ=255)

                                                                                                                                                                                                                                            tt‑tGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGTAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggta                                                                                                    <  1:488440/144‑1 (MQ=255)

                                                                                                                                                                                                                                            tt‑tGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggta                                                                                                    >  2:488440/1‑144 (MQ=255)

                                                                                                                                                                                                                                                gCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑GTGCGTGCCATCAGGCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTcc  <  2:428035/239‑1 (MQ=255)

CATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGC‑ATCGTTTTAACGCGCGCAATAGCATCAAAGGTATT‑TGCCGGATGAC‑A‑GT‑GAAACTCTGTTAC‑TTGCGTGCCATCAGTCGATGCTGC‑AGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCA‑GGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCC  >  NC_000913/3135206‑3135682

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