BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	553,887	G→C	6.4%	I18M (ATC→ATG) ‡	lpxH ←	UDP‑2,3‑diacylglucosamine pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	553,887	0	G	C	6.4%	859.0 / 40.5	251	I18M (ATC→ATG) ‡	lpxH	UDP‑2,3‑diacylglucosamine pyrophosphohydrolase
Reads supporting (aligned to +/- strand): ref base G (104/130); new base C (6/10); total (111/140)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-03

Only 100 of 254 total aligned reads displayed.

ACCCGGCGACCATAAAGTTCGAGCACCTTTTCTTCCGGCAATAACGTCATGCCACTTTCACGGGCAAAGCGTTTGCCGAGCAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTG  >  NC_000913/553612‑554149

acCCGGCGACCATAAAGTTCGAGCACCTTTTCTTCCGGCAATAACGTCATGCCACTTTCACGGGCAAAGCGTTTGCCGAGCAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGAt                                                                                                                                                                                                                                                                      >  1:226388/1‑278 (MQ=255)

     gCGACCATAAAGTTCGAGCACCTTTTCTTCCGGCAATAACGTCATGCCACTTTCACGGGCAAAGCGTTTGCCGAGCAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCCAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAAACGGCGGTGCTCCCCg                                                                                                                                                                                                                                                                 >  1:182384/1‑278 (MQ=255)

                     agtaCCTTTTCGTCCGGCAATAACGTCATGCGACTTTCACGGGCAAAGCGTTTGCCGAGCAGAAAATAACGGTTGCCATGAATGAAATAACAGGGAACGTCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCa                                                                                                                                                                                                                                                    <  2:473568/272‑1 (MQ=255)

                     aGCACCTTTTCTTCCGGCAATAACGTCATGCCACTTTCACGGGCAAAGCGTTTGCCGAGCAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGTGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCa                                                                                                                                                                                                                                                    <  2:131347/275‑1 (MQ=255)

                                    ggCAATAACGTCATGCCACTTTCACGGGCAAAGCGTTTGCCGAGCAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGa                                                                                                                                                                                                                                             >  4:303221/1‑267 (MQ=255)

                                                       tttCACGGGCAAAGCGTTTGCCGAGCAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAAAGCAGAAAACCGGCAGTTATCGCCGGTTCCTCCACGCAGACATGAAGCTCTGCCATAAAGGGTGTCGcacagatta                                                                                                                                                                                                              >  1:104501/1‑271 (MQ=255)

                                                         tCACGGGCAAAGCGTTTGCCGAGCAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGa                                                                                                                                                                                                                 >  2:244542/1‑274 (MQ=255)

                                                            cGGGCAAAGCGTTTGCCGAGCAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTa                                                                                                                                                                                                              >  2:275437/1‑274 (MQ=255)

                                                                               gCAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCttt                                                                                                                                                                                          >  2:392610/1‑275 (MQ=255)

                                                                               gCAAAAAATCACGGTTTCGATGAATGAAATAACAGGGAACGACGAGATCGGAACCCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGTGTCGTCGTCGCCAAACCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCtt                                                                                                                                                                                           <  4:244539/274‑1 (MQ=255)

                                                                                                       cgaaATAACCGGGCGCGCCGGAATCGGACACCGCATTGATGGACCCCGCCATTTTGCGATGGAGTGGGTGGAGAACGCCGCCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGAGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCa                                                                                                                                                                  <  4:392606/274‑1 (MQ=255)

                                                                                                                   acacccctgcaTCCGACACCCCTTTTTCCGCCGCCTCCCTCTTTGTATGGAGTGGGTTGGGATCGTCGTGCCCAATCCNTGTTTCACACAGATCGCCAAGAATATCCAGCGCGTCGGCCTGGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGTTTCTTCCACGCAGAGGTGAAGATCTGCAATAAAGCGTGTCGCCACGATTACTCGCATACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCAt                                                                                                                                                     <  3:226388/266‑1 (MQ=255)

                                                                                                                                 gaaacCGCTTTGATGGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGAcc                                                                                                                                         <  2:323800/271‑1 (MQ=255)

                                                                                                                                 gACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCTCTAAAAAACGCAGAAAACCGGCGGTGAt                                                                                                                                                                                                                                                                      <  6:1550373/149‑1 (MQ=255)

                                                                                                                                 gACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGAt                                                                                                                                                                                                                                                                      <  7:249094/149‑1 (MQ=255)

                                                                                                                                 gACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGAt                                                                                                                                                                                                                                                                      >  5:1467566/1‑149 (MQ=255)

                                                                                                                                    acCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGc                                                                                                                                                                                                                                                                   >  7:1075647/1‑149 (MQ=255)

                                                                                                                                     aggCTTTGATCGCCGCCGCCAGCTTGCGATGGAAGTGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGt                                                                                                                                      <  2:74049/271‑1 (MQ=255)

                                                                                                                                       gCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATc                                                                                                                                                                                                                                                                     >  7:314905/1‑144 (MQ=255)

                                                                                                                                       gCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCgg                                                                                                                                                                                                                                                                >  5:1220479/1‑149 (MQ=255)

                                                                                                                                         tttGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGtt                                                                                                                                                                                                                                                              <  7:973227/149‑1 (MQ=255)

                                                                                                                                                gccgccgccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGttt                                                                                                                                                                                                                                                             <  6:1127573/143‑2 (MQ=255)

                                                                                                                                                gccgccgccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGttt                                                                                                                                                                                                                                                             >  5:1127573/1‑142 (MQ=255)

                                                                                                                                                   gccgccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGAtt                                                                                                                      <  2:63854/275‑1 (MQ=255)

                                                                                                                                                       ccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAgaga                                                                                                                                                                                                                                                <  6:8115/149‑1 (MQ=255)

                                                                                                                                                        cATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGAt                                                                                                                                                                                                                                               <  6:1467566/149‑1 (MQ=255)

                                                                                                                                                              gCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGAt                                                                                                                                                                                                                                         <  8:1643240/149‑1 (MQ=255)

                                                                                                                                                              gCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGAt                                                                                                                                                                                                                                         <  8:156489/149‑1 (MQ=255)

                                                                                                                                                                      gTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATa                                                                                                                                                                                                                                 >  8:1197475/1‑149 (MQ=255)

                                                                                                                                                                        gggTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATaaa                                                                                                                                                                                                                               <  5:210607/149‑1 (MQ=255)

                                                                                                                                                                         ggTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAag                                                                                                                                                                                                                              <  6:1220479/149‑1 (MQ=255)

                                                                                                                                                                         ggTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAag                                                                                                                                                                                                                              <  5:837009/149‑1 (MQ=255)

                                                                                                                                                                          gTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAaga                                                                                                                                                                                                                             <  5:1616030/149‑1 (MQ=255)

                                                                                                                                                                           ttGGCATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAagag                                                                                                                                                                                                                            <  5:77401/149‑1 (MQ=255)

                                                                                                                                                                            tGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAgt                                                                                                                                                                                                                           <  8:1109556/149‑1 (MQ=255)

                                                                                                                                                                              gggtcgtcgtcgCCAATCCATACTTAAAACAGATCGCCAAGAATATACAGCGCGTGGGCCTTGCGGGCTTCACCCGCTAAAAAACGCAGAAAACCGGCGGTGATCACCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAgtgt                                                                                                                                                                                                                         >  8:1048725/4‑149 (MQ=255)

                                                                                                                                                                               gATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTCAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGAt                                                                                                                                                                                                                                         >  3:56361/1‑132 (MQ=255)

                                                                                                                                                                               gATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGAt                                                                                                                                                                                                                                         <  3:142970/132‑1 (MQ=255)

                                                                                                                                                                               gATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGAt                                                                                                                                                                                                                                         >  1:142970/1‑132 (MQ=255)

                                                                                                                                                                                 tcgtcgtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGc                                                                                                                                                                                                                      >  6:398519/1‑149 (MQ=255)

                                                                                                                                                                                  cgtcgtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGcc                                                                                                                                                                                                                     <  5:999170/149‑1 (MQ=255)

                                                                                                                                                                                   gtcgtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAg                                                                                                                                                                                                                                           >  2:436775/1‑126 (MQ=255)

                                                                                                                                                                                    tcgtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCACGTCGGCCCTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCAc                                                                                                                                                                                                                   <  7:1455114/149‑1 (MQ=255)

                                                                                                                                                                                      gtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGa                                                                                                                                                                                                                 <  6:288145/149‑1 (MQ=255)

                                                                                                                                                                                      gtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGa                                                                                                                                                                                                                 >  5:698288/1‑149 (MQ=255)

                                                                                                                                                                                      ctcgCCAATCCCGGCTGCAAACAGATCGACACGAAGATGCAGCGCTTCGGCCGTGCGGGCGTCCACCGCTAAGAAACGCAGAAAACCGGCGGTGATCGCCGGTTGTTCCACGCAGGGATGAAGTCACGCAATAAAGGGTGTCGGCACGATTACTCGCTAACCGTCATGCTTGCAATGCTACCGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACAGGTTGCAACACCGTTGATTTTGTCTACCACGCCCATGCCGTCAACCACTTCAGc                                                                                   <  4:122348/274‑1 (MQ=255)

                                                                                                                                                                                       tcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGAt                                                                                                                                                                                                                >  8:74652/1‑149 (MQ=255)

                                                                                                                                                                                         gCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTa                                                                                                                                                                                                              >  8:27718/1‑149 (MQ=255)

                                                                                                                                                                                         gCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTa                                                                                                                                                                                                              <  6:683949/149‑1 (MQ=255)

                                                                                                                                                                                         gCCAATCCAAGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTa                                                                                                                                                                                                              <  7:1424478/149‑1 (MQ=255)

                                                                                                                                                                                            aaTCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTc                                                                                                                                                                                                           <  6:644673/149‑1 (MQ=255)

                                                                                                                                                                                            aaTCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGAGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGc                                                                                                                                                                   >  1:55256/1‑189 (MQ=255)

                                                                                                                                                                                                cATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTa                                                                                                                                                                                                       >  5:1093248/1‑149 (MQ=255)

                                                                                                                                                                                                 aTGCTTCAAACAGATCGCCCAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCACGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTaa                                                                                                                                                                                                      >  8:1165183/1‑149 (MQ=255)

                                                                                                                                                                                                 aTGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTTCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTaa                                                                                                                                                                                                      <  8:501058/149‑1 (MQ=255)

                                                                                                                                                                                                 aTGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTaa                                                                                                                                                                                                      >  7:787335/1‑149 (MQ=255)

                                                                                                                                                                                                 aTGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTaa                                                                                                                                                                                                      <  8:886857/149‑1 (MQ=255)

                                                                                                                                                                                                 aTGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTaa                                                                                                                                                                                                      >  5:159113/1‑149 (MQ=255)

                                                                                                                                                                                                        aaaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCAc                                                                                                                                                                                               >  5:349833/1‑149 (MQ=255)

                                                                                                                                                                                                        aaaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCAc                                                                                                                                                                                               >  6:232129/1‑149 (MQ=255)

                                                                                                                                                                                                        aaaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCAc                                                                                                                                                                                               <  5:675437/149‑1 (MQ=255)

                                                                                                                                                                                                              aTCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGtt                                                                                                                                     >  4:559541/1‑201 (MQ=255)

                                                                                                                                                                                                                  ccAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATg                                                                                                                                                                                     >  6:49177/1‑149 (MQ=255)

                                                                                                                                                                                                                   cAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGa                                                                                                                                                                                    >  5:1275110/1‑149 (MQ=255)

                                                                                                                                                                                                                    aaGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGAt                                                                                                                                                                                   >  7:1235033/1‑149 (MQ=255)

                                                                                                                                                                                                                    aaGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGttt                                                                                                                                                                                                                                                             >  5:416418/1‑74 (MQ=38)

                                                                                                                                                                                                                    aaGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGttt                                                                                                                                                                                                                                                             >  5:911520/1‑74 (MQ=38)

                                                                                                                                                                                                                    aaGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGttt                                                                                                                                                                                                                                                             <  6:416418/75‑2 (MQ=38)

                                                                                                                                                                                                                      gAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCGATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGc                                                                                                                                                                                             >  6:1371326/1‑137 (MQ=255)

                                                                                                                                                                                                                      gAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGc                                                                                                                                                                                             <  5:1371326/137‑1 (MQ=255)

                                                                                                                                                                                                                      gAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATaa                                                                                                                                                                                 <  6:1620124/149‑1 (MQ=255)

                                                                                                                                                                                                                       aaTATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAAc                                                                                                                                                                                <  6:1395789/149‑1 (MQ=255)

                                                                                                                                                                                                                         tataCAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGt                                                                                                                                                                              <  8:1332019/149‑1 (MQ=255)

                                                                                                                                                                                                                            aCAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGtctt                                                                                                                                                                           <  5:1646493/149‑1 (MQ=255)

                                                                                                                                                                                                                               gcgcGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTcttctt                                                                                                                                                                        >  6:207097/1‑149 (MQ=255)

                                                                                                                                                                                                                                 gcgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACCCGCAGTCGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCa                                    >  1:95646/1‑279 (MQ=255)

                                                                                                                                                                                                                                  cgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTg                                                                                                                                                                     <  7:131107/149‑1 (MQ=255)

                                                                                                                                                                                                                                   gTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTgg                                                                                                                                                                    >  8:556778/1‑149 (MQ=255)

                                                                                                                                                                                                                                    tCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCa                                                                                                                                                                  <  2:429622/150‑1 (MQ=255)

                                                                                                                                                                                                                                          ttGCGGGCTTCCCCCGGTAAAAAACGCAAAAAACCGGCGGTGATCGCAGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATACCGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTGGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCg                            <  1:291951/278‑1 (MQ=255)

                                                                                                                                                                                                                                              gggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGt                                                                                                                                                         >  8:393308/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  ttCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACg                                                                                                                                                                                                                  >  4:75154/1‑88 (MQ=255)

                                                                                                                                                                                                                                                  ttCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCAt                                                                                                                                                     >  7:1597662/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  ttCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCAt                                                                                                                                                     >  6:874413/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    cccccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATAc                                                                                                                                                   <  5:49177/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      cccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGc                                                                                                                                                       >  7:1045567/1‑143 (MQ=255)

                                                                                                                                                                                                                                                      cccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCg                                                                                                                                                 >  7:1165793/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      cccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCg                                                                                                                                                 >  6:621847/1‑149 (MQ=255)

                                                                                                                                                                                                                                                         gCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTa                                                                                                                                              >  8:1515324/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           tAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGt               <  2:548961/274‑1 (MQ=255)

                                                                                                                                                                                                                                                            aaaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGa                                                                                                                                           >  6:483120/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 aCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGt                                                                                                                                      <  7:1143096/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 aCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGt                                                                                                                                      <  7:556778/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     gaagaaCCGGCGATCACCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAAt                                                                                                                                                                                      <  8:1667509/93‑1 (MQ=38)

                                                                                                                                                                                                                                                                     gaagaaCCGGCGATCACCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAAt                                                                                                                                                                                      >  7:1667509/5‑97 (MQ=39)

                                                                                                                                                                                                                                                                         aaCCGGGGGTGACTGCCGGTTCTTCCACCCAGAGAGGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACtg  <  4:485635/273‑1 (MQ=255)

                                                                                                                                                                                                                                                                          aCCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGAc                                                                                                                                          >  7:271090/1‑136 (MQ=255)

                                                                                                                                                                                                                                                                          aCCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGAc                                                                                                                                          <  8:271090/136‑1 (MQ=255)

                                                                                                                                                                                                                                                                               cggTGATCGCCGGTTCTTCCAAGCAGGGAAGAAGGACTGGAATAAAGGGTGTTGCCAAGAATAATTCGTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGa                                                                                                                        <  7:941291/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                                   gATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGAtttt                                                                                                                    <  6:1093248/149‑1 (MQ=255)

ACCCGGCGACCATAAAGTTCGAGCACCTTTTCTTCCGGCAATAACGTCATGCCACTTTCACGGGCAAAGCGTTTGCCGAGCAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTG  >  NC_000913/553612‑554149

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`