BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,344,036	G→C	5.6%	T104S (ACT→AGT)	yciT ←	global regulator of transcription; DeoR family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,344,036	0	G	C	5.6%	624.3 / 22.9	179	T104S (ACT→AGT)	yciT	global regulator of transcription; DeoR family
Reads supporting (aligned to +/- strand): ref base G (82/87); new base C (8/2); total (90/89)
Fisher's exact test for biased strand distribution p-value = 9.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-03

Only 100 of 183 total aligned reads displayed.

TTTTCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGACCCTGGT  >  NC_000913/1343769‑1344306

ttttCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCgg                                                                                                                                                                                                                                                                           >  2:316554/1‑273 (MQ=255)

                                          cgACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTTCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑ATGGCCAACCATACTTTCGCTTTTTTTCTGATCCACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCCTACCGTGATGATAGTGCACTTTTTCTTCCACTCGCCCAAAGACCGGGCGATTAGGGCATTTCGTCTGCCATTTTGGATAAACCTGGTTTc                                                                                                                                                                                                                               >  4:490927/1‑275 (MQ=255)

                                            aCCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGc                                                                                                                                                                                                                             <  4:455690/275‑1 (MQ=255)

                                                                   gccagccaTCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCa                                                                                                                                                                                                      >  1:81289/1‑275 (MQ=255)

                                                                                       aaTGGTTTGCTGAAAGGTACCTGTTGGATGCACGTGC‑TGTTAATAGGGCCAACCTTACTTTCGGTTTTTTTCTGATACACGCCACCGAGAAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCtt                                                                                                                                                                                  <  4:559498/275‑1 (MQ=255)

                                                                                       aaTGCTTTTCTGAATTGCACCTTTTGGATGCACTGCCGTGTCAA‑AGGGCAAACCAGACTTTCGCTTTTTTTCTGTTACACGCCACCGAGCAAAATAACTTCACAAGGCCCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCTAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCtt                                                                                                                                                                                  <  4:533809/275‑1 (MQ=255)

                                                                                         tgtttttCTTAAATTTACCTTTTGGATGCACTGTCGTGTCAA‑ATGGCCAACCCTACTTTCGTTTTTTTTCTGCTACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCa                                                                                                                                                                                <  4:491575/272‑1 (MQ=255)

                                                                                             ttGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGtt                                                                                                                                                                                                                                 >  4:641819/1‑222 (MQ=255)

                                                                                              tGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACGCGGCGAGTTCACGCTTCAGCGta                                                                                                                                                                           >  2:396656/1‑275 (MQ=255)

                                                                                                        cACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGtt                                                                                                                                                                                                                                 >  2:591165/1‑211 (MQ=255)

                                                                                                              ttGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATa                                                                                                                                                                                                                                          >  2:451169/1‑196 (MQ=255)

                                                                                                                ggATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCATCGTATAGTTGCTCATCATatg                                                                                                                                                          >  1:12745/1‑272 (MQ=255)

                                                                                                                        tGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGt                                                                                                                                                                                                                                                                               >  5:43672/1‑149 (MQ=255)

                                                                                                                        tGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGt                                                                                                                                                                                                                                                                               >  5:120539/1‑149 (MQ=255)

                                                                                                                        tGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGAc                                                                                                                                                 >  4:61717/1‑275 (MQ=255)

                                                                                                                           cGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATa                                                                                                                                                                                                                                          >  2:523651/1‑183 (MQ=255)

                                                                                                                             tgtCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTACTCGCCCGGACTCTggg                                                                                                                                                                                                                                                                          <  7:827969/149‑1 (MQ=255)

                                                                                                                               ccAA‑AGGGCCCACCATATTTTCGCTTTTTTTCTGATACACGCCACTGAGCAAAATAACTTCTCAAGGTTCGTCTTTCAGCAAATGCGCGACGTAGCTGCTGACCGTGCTGATAGTGCCATTTTTCTTCTGCACGCCCCGAGTCTGGGCGAGTAGGGAATTGCTGCTGCCATTTTCGACAAAGATGGTTTCGCCCGGTTGAAACAGCGACGGGGCAAACTCGGCGAGTTCACGATTCCGAGTATTGTTGCACAACATCGGGGTTTCGagggaattc                                                                                                                                          >  4:593511/2‑267 (MQ=255)

                                                                                                                                    aGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGAgca                                                                                                                                                                                                                                                                    <  8:29721/149‑3 (MQ=255)

                                                                                                                                      ggCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAgg                                                                                                                                                                                                                                                                  <  6:455523/149‑1 (MQ=255)

                                                                                                                                      ggCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAgg                                                                                                                                                                                                                                                                  <  5:261386/149‑1 (MQ=255)

                                                                                                                                      ggCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTACTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAgg                                                                                                                                                                                                                                                                  <  8:6886/149‑1 (MQ=255)

                                                                                                                                         cAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCa                                                                                                                                                                                                                                                               >  8:104649/1‑149 (MQ=255)

                                                                                                                                           ccctttttttctttttttttCANNAACACNACCCCGGCAAAAATACCTGCACCAGGGGCGCGTTGCAGCAAATGGCCGATTTGGCTGCTGACCGTGATTATGGTGTCCTTTTGCTTATGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGa                                                                                                                               <  3:9082/262‑1 (MQ=255)

                                                                                                                                            ccATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtg                                                                                                                                                                                                                                                            <  5:309311/149‑1 (MQ=255)

                                                                                                                                                  tttCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgctg                                                                                                                                                                                                                                                      >  5:796107/1‑149 (MQ=255)

                                                                                                                                                  tttCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCa                                                                                                                        <  2:61717/275‑1 (MQ=255)

                                                                                                                                                       cTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGTTTTCGCCCGGTTGAACCAGCGACGCGCCAAACCCGGCGAGTTCACGCCTCAGCGTTTACTTGCTCAACATAGGGGTTTCGACCTCCTCAATATCAAGCGAAACACCAAAccc                                                                                                                   >  4:621565/1‑275 (MQ=255)

                                                                                                                                                        tttttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATGGCTACTGCCAttt                                                                                                                                                                                                                                                >  6:850423/1‑149 (MQ=255)

                                                                                                                                                        tttttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGAAGATCGTGACATTTTACTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAttt                                                                                                                                                                                                                                                >  6:789976/1‑149 (MQ=255)

                                                                                                                                                          tttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTc                                                                                                                                                                                                                                              <  6:125961/149‑1 (MQ=255)

                                                                                                                                                           ttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCg                                                                                                                                                                                                                                             >  8:458545/1‑149 (MQ=255)

                                                                                                                                                                tGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCgg                                                                                                                                                                                                                                                                           >  6:365187/1‑114 (MQ=255)

                                                                                                                                                                  aTACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGt                                                                                                                                                                                                                                  >  4:328183/1‑153 (MQ=255)

                                                                                                                                                                   tACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGAt                                                                                                                                                                                                                                     <  5:445421/149‑1 (MQ=255)

                                                                                                                                                                    acacGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATg                                                                                                                                                                                                                                    <  6:909521/149‑1 (MQ=255)

                                                                                                                                                                     cacGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATgg                                                                                                                                                                                                                                   <  8:172912/149‑1 (MQ=255)

                                                                                                                                                                     cacGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATgg                                                                                                                                                                                                                                   <  8:318338/149‑1 (MQ=255)

                                                                                                                                                                      acGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGt                                                                                                                                                                                                                                  <  6:284090/149‑1 (MQ=255)

                                                                                                                                                                      acGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCg                                                                                                                                                                                                                              <  2:328183/153‑1 (MQ=255)

                                                                                                                                                                         ccACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                                                                                                                                                                                               <  7:695660/149‑1 (MQ=255)

                                                                                                                                                                         ccACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                                                                                                                                                                                               <  6:562139/149‑1 (MQ=255)

                                                                                                                                                                         ccACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                                                                                                                                                                                               <  6:168834/149‑1 (MQ=255)

                                                                                                                                                                           aCCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGc                                                                                                                                                                                                                             >  7:1017817/1‑149 (MQ=255)

                                                                                                                                                                            ccGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                                                                                                                                                                                               >  8:908600/1‑146 (MQ=255)

                                                                                                                                                                                 cAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGtt                                                                                                                                                                                                                       <  7:708495/149‑1 (MQ=255)

                                                                                                                                                                                   aaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGa                                                                                                                                                                                                                     <  6:326785/149‑1 (MQ=255)

                                                                                                                                                                                   aaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGa                                                                                                                                                                                                                     >  7:117744/1‑149 (MQ=255)

                                                                                                                                                                                   aaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGa                                                                                                                                                                                                                     <  6:736022/149‑1 (MQ=255)

                                                                                                                                                                                    aaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGt                                                                                                                                                                       >  4:246584/1‑194 (MQ=255)

                                                                                                                                                                                       aaCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCa                                                                                                                                                                                                                 >  7:177563/1‑149 (MQ=255)

                                                                                                                                                                                        aCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCt                                                                                                                                                                   <  2:246584/194‑1 (MQ=255)

                                                                                                                                                                                         cTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGc                                                                                                                                                                                                               <  5:765543/149‑1 (MQ=255)

                                                                                                                                                                                         cTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATAATACTGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTCACACCCTTTTTTGaag                                                                                 >  1:189710/1‑273 (MQ=255)

                                                                                                                                                                                          ttCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCg                                                                                                                                                                                                              >  7:615363/1‑149 (MQ=255)

                                                                                                                                                                                             acaAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATgg                                                                                                                                                                                                                                   >  5:64037/1‑125 (MQ=255)

                                                                                                                                                                                             acaAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATgg                                                                                                                                                                                                                                   <  6:64037/125‑1 (MQ=255)

                                                                                                                                                                                                    gcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTa                                                                                                                                                                                                                                                                    >  7:453015/1‑85 (MQ=255)

                                                                                                                                                                                                    gcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGc                                                                                                                                                                                                                             <  7:225478/124‑1 (MQ=255)

                                                                                                                                                                                                    gcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGc                                                                                                                                                                                                                             >  8:225478/1‑124 (MQ=255)

                                                                                                                                                                                                    gcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCaaa                                                                                                                                                                                                    >  8:884334/1‑149 (MQ=255)

                                                                                                                                                                                                     cgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATa                                                                                                                                                                                                                                          >  4:447907/1‑110 (MQ=255)

                                                                                                                                                                                                            cAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGa                                                                                                                                                                                            <  6:386968/149‑1 (MQ=255)

                                                                                                                                                                                                             aGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAg                                                                                                                                                                                           >  6:258556/1‑149 (MQ=255)

                                                                                                                                                                                                              gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAt                                                                                                                                                                                                                                           >  5:71839/1‑100 (MQ=255)

                                                                                                                                                                                                              gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAt                                                                                                                                                                                                                                           <  6:71839/100‑1 (MQ=255)

                                                                                                                                                                                                              gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGt                                                                                                                                                                                          >  8:416903/1‑149 (MQ=255)

                                                                                                                                                                                                              gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGt                                                                                                                                                                                          <  7:369016/149‑1 (MQ=255)

                                                                                                                                                                                                                aaaTGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTc                                                                                                                                                                                        >  5:157722/1‑149 (MQ=255)

                                                                                                                                                                                                                      gcgATGTAGCTGCTGACCGTGATGATAGTGACCTTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCAGCAAACTCGGCGAGTTCACGCtt                                                                                                                                                                                  <  5:913145/149‑1 (MQ=255)

                                                                                                                                                                                                                        gATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCa                                                                                                                                                                                >  5:1020815/1‑149 (MQ=255)

                                                                                                                                                                                                                               ctgctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGtata                                                                                                                                                                         <  7:458545/149‑1 (MQ=255)

                                                                                                                                                                                                                                     aCCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCt                                                                                                                                                                   >  7:789376/1‑149 (MQ=255)

                                                                                                                                                                                                                                        gTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTg                                                                           <  2:633351/234‑1 (MQ=255)

                                                                                                                                                                                                                                           atgatAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTcatcat                                                                                                                                                             >  7:950645/1‑149 (MQ=255)

                                                                                                                                                                                                                                              atAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGTCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACg                                                                                                                                                          >  8:599966/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGttt                                                                                                                                                     <  6:43672/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    aCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTTCTGCCATTTTCGATAAAGATGGGTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTAAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCgg                       <  2:621565/274‑1 (MQ=255)

                                                                                                                                                                                                                                                    aCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTc                                                                                                                                                    <  8:254222/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    aCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTCTCGCCCGGTTGAACCAGCGACGCGGCAAACCCGGCGAGTTCCCGCTTCAGCGTATAGTTGctcctc                                                                                                                                                               >  4:132674/1‑138 (MQ=255)

                                                                                                                                                                                                                                                      aTTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGa                                                                                                                                                  >  6:839289/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       tttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGAc                                                                                                                                                 <  8:442221/149‑1 (MQ=255)

                                                                                                                                                                                                                                                        ttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCGCGCTTCAGCGTATAGTTGCTCATCATAc                                                                                                                                                           <  2:132674/138‑1 (MQ=255)

                                                                                                                                                                                                                                                        ttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGa                                                                                   <  2:87550/210‑1 (MQ=255)

                                                                                                                                                                                                                                                            gccctacTCGCCCGGACTCTGGGCGAGCAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcat                                                                                                                                            >  6:880422/7‑149 (MQ=255)

                                                                                                                                                                                                                                                            cttctGCTCGCCCAGAGTCCTGGCGAGTAGGGCATTGCTGCTGCCATTATCGATAAAGATGGTTTCGACCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcat                                                                                                                                            <  5:850423/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             ttctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcatc                                                                                                                                           >  8:118125/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             ttctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcatc                                                                                                                                           >  7:255267/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             ttctGCTCGCCCAGAGTCC‑GGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcatca                                                                                                                                          >  8:858166/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              cctacTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcatca                                                                                                                                          >  5:349296/5‑149 (MQ=255)

                                                                                                                                                                                                                                                                tacTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACt                                                                                                                                        >  7:252985/3‑149 (MQ=255)

                                                                                                                                                                                                                                                                    cGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAGACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGAccc      <  2:593511/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                     gCCCGGACTCTGGGCGAGCAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTTGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCaa                                                                                                                                   >  6:541576/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     gCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCaa                                                                                                                                   <  5:1054005/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                       ccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGc                                                                                                                                 >  5:69860/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       ccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGc                                                                                                                                 <  7:394078/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                        cAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGACCCTGGt  <  2:490927/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                         agagTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGaa                                                                                                                                                                                                                    >  7:907029/1‑64 (MQ=255)

                                                                                                                                                                                                                                                                           agTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGaaa                                                                                                                             <  6:327539/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                            gTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAAc                                                                                                                            <  6:144909/149‑1 (MQ=255)

TTTTCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGACCCTGGT  >  NC_000913/1343769‑1344306

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