BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,470,952	T→C	7.1%	intergenic (‑39/+34)	ridA ← / ← pyrI	enamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,470,952	0	T	C	7.1%	427.1 / 15.4	134	intergenic (‑39/+34)	ridA/pyrI	enamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand): ref base T (62/62); new base C (6/4); total (68/66)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.49e-01

Only 100 of 137 total aligned reads displayed.

GGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGG  >  NC_000913/4470689‑4471206

  tAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAg                                                                                                                                                                                                                                                   <  2:700897/275‑1 (MQ=255)

  gaGTTTTCACGGTTGCGGCCTCTTTCAGGGCCGGGACTTCGCCCATCGGTTTTTCGTTATCCAGCGACTTACGTGCCTGTGCAGCTTCTTCTGCCGGTACTTCGCCCGTTTGCGGAGTTACCTGGATCTGACCGGAGGTGATGATCATATTGCCCAGTTCAACGCCCTATACGTAAGGACCGATAGCTTCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAg                                                                                                                                                                                                                                                   <  4:679252/274‑1 (MQ=255)

        tAACGATGTCGCCCACGTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                                                                                                                                                                             <  1:140426/275‑1 (MQ=255)

               gTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTCCGTAAGGACCGATCGCTGCCGGTGCCTTTTCCGTCTCGATAGTTTTGACCCTGAATTCTCCTTTATTACAGTGGTAAAAAGGCTGGCGCATTATAGGGAGCCGGACTTTGTTTAccc                                                                                                                                                                                                                                      >  4:699476/1‑274 (MQ=255)

                                             gaTCGCTTTTACGTTATCCAGCGACTGACGCGCCTGGGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTa                                                                                                                                                                                                         <  2:699478/274‑1 (MQ=255)

                                                             tCCAGCGACTGACGTGCCTGGGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttttt                                                                                                                                                                                         <  1:70266/274‑1 (MQ=255)

                                                                      tGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGc                                                                                                                                                                                                                      <  1:40430/236‑1 (MQ=255)

                                                                                 tgCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg                                                                                                                                                                      <  2:100759/273‑1 (MQ=255)

                                                                                      cGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt                                                                                                                                                                          >  4:72639/1‑264 (MQ=255)

                                                                                          gTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg                                                                                                                                                                      <  2:72640/264‑1 (MQ=255)

                                                                                           tCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAt                                                                                                                                                           >  2:128455/1‑274 (MQ=255)

                                                                                                 ggTACTTCGCCCGTTTTCGGATTTACCGGTATCTGACCGGANGTGATGCTCATATTGCCCAGTTCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAg                                                                                                                                                                                                                                                   <  3:130182/180‑1 (MQ=255)

                                                                                                       tCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGCCTTTTATTACCAACCGCAATTAATTGGCCAGCACCCCATTATGTGAAAACTCTTTTTCACAGTTTTTGCTTTTGATCGCGATCTCTTTGGCGCTTTTTcc                                                                                                                                              >  4:427411/1‑274 (MQ=255)

                                                                                                             gTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCCCCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTGCCACGGc                                                                                                                                        >  2:312121/1‑275 (MQ=255)

                                                                                                             gTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                                                                                                                         >  2:652835/1‑274 (MQ=255)

                                                                                                                 ttggATTTACCGGGATCTGACCGTGGGTGATGATCATATTTCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTTCATTTTCCGTCGCGATAGTTTTGCTCATGATGTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAg                                                                                                                                    <  4:407882/273‑1 (MQ=255)

                                                                                                                     aTTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAtt                                                                                                                                                                                                                                                               >  5:199175/1‑148 (MQ=255)

                                                                                                                        tACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGTTCATTATAg                                                                                                                                                                                                                                                           <  8:111825/149‑1 (MQ=255)

                                                                                                                          ccGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAgg                                                                                                                                                                                                                                                          <  6:849503/148‑1 (MQ=255)

                                                                                                                           cGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGa                                                                                                                                                                                                                                                        <  8:565559/149‑1 (MQ=255)

                                                                                                                           cGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGa                                                                                                                                                                                                                                                        >  8:153683/1‑149 (MQ=255)

                                                                                                                            gggATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGa                                                                                                                                                                                                                                                        >  5:218983/1‑148 (MQ=255)

                                                                                                                            gggATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGa                                                                                                                                                                                                                                                        >  5:539971/1‑148 (MQ=255)

                                                                                                                              gATCTGACCGGAGGTGCTGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTAGCTCATTATAGGGAGc                                                                                                                                                                                                                                                      >  6:63132/1‑148 (MQ=255)

                                                                                                                               cgcTGACCGGAGGTTATGATCAGATTGCCCAGATCAACGCCCTGTACGTAAGAACAGATAGCTGCCGGTGCATTTTCCGTCTCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGCCGCATTATAGGGAGCCAGAATTTTATTACCAACCGCAATTAATTGACCAGGACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                                                                                                                      <  1:101861/273‑1 (MQ=255)

                                                                                                                                  tGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGa                                                                                                                                                                                                                                                  <  6:893760/148‑1 (MQ=255)

                                                                                                                                  tGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGa                                                                                                                                                                                                                                                  <  5:855145/148‑1 (MQ=255)

                                                                                                                                    aCCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtt                                                                                                                                                                                                                                               <  7:334153/149‑1 (MQ=255)

                                                                                                                                      cGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACttt                                                                                                                                                                                                                                              >  5:115800/1‑148 (MQ=255)

                                                                                                                                        gAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTtta                                                                                                                                                                                                                                            <  6:19688/148‑1 (MQ=255)

                                                                                                                                         aGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaat                                                                                                                                                                                                                          >  2:754741/1‑165 (MQ=255)

                                                                                                                                                 atCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                                                                                                                                                                                   >  5:359317/1‑148 (MQ=255)

                                                                                                                                                   cATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGc                                                                                                                                                                                                                                 >  6:869609/1‑148 (MQ=255)

                                                                                                                                                    atatTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCa                                                                                                                                                                                                                                <  6:416142/148‑1 (MQ=255)

                                                                                                                                                      atTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaat                                                                                                                                                                                                                              >  5:150925/1‑148 (MQ=255)

                                                                                                                                                           ccAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                                                                                                                                                                         <  8:731956/148‑1 (MQ=255)

                                                                                                                                                           ccAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                                                                                                                                                                         >  7:731956/1‑148 (MQ=255)

                                                                                                                                                            cAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                                                                                                                                                                       >  8:478770/1‑149 (MQ=255)

                                                                                                                                                            cAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTg                                                                                                                                                                                                                        >  6:37113/1‑148 (MQ=255)

                                                                                                                                                             aGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGTAATTAATTGGCCAGCACCACATTATGGGAAAACTATTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGTAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCgg                                                                                        >  2:439125/1‑275 (MQ=255)

                                                                                                                                                             aGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGc                                                                                                                                                                                                                      >  7:879614/1‑149 (MQ=255)

                                                                                                                                                              ggtcaaagcgcgttACGGAAGGACCGATAGCTTCCGGGTCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCggg                                                                                       <  4:312120/262‑1 (MQ=255)

                                                                                                                                                              gATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGc                                                                                                                                                                                                                      >  5:618491/1‑148 (MQ=255)

                                                                                                                                                               aTCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAg                                                                                                                                                                                                                                                   <  8:646151/118‑1 (MQ=255)

                                                                                                                                                               aTCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAg                                                                                                                                                                                                                                                   >  7:646151/1‑118 (MQ=255)

                                                                                                                                                                   aCGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGAc                                                                                  <  2:74071/275‑1 (MQ=255)

                                                                                                                                                                              gTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATggg                                                                                                                                                                                                     >  7:876160/1‑149 (MQ=255)

                                                                                                                                                                                aaGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                                                                                                                                                                         >  5:914016/1‑127 (MQ=255)

                                                                                                                                                                                  ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaat                                                                                                                                                                                                                          <  7:82664/124‑1 (MQ=255)

                                                                                                                                                                                  ggACCGATAGCAGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaat                                                                                                                                                                                                                          >  8:82664/1‑124 (MQ=255)

                                                                                                                                                                                   gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                                                                                                                                                                         >  2:561838/1‑124 (MQ=255)

                                                                                                                                                                                     ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAc                                                                                                                                                                                                                                    <  2:677763/111‑1 (MQ=255)

                                                                                                                                                                                     ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActc                                                                                                                                                                                              >  7:330810/1‑149 (MQ=255)

                                                                                                                                                                                            cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                                                                                                                                                                        <  5:739553/148‑1 (MQ=255)

                                                                                                                                                                                            cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                                                                                                                                                                        <  6:41751/148‑1 (MQ=255)

                                                                                                                                                                                             tGCCGGTGCATTTTCCGTCGCGATAGCTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCTAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCTCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCa                                                        >  2:62133/1‑275 (MQ=255)

                                                                                                                                                                                             tGCCGGTGCATTTTCCGTCACGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCa                                                        <  1:100458/275‑1 (MQ=255)

                                                                                                                                                                                              gCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcaca                                                                                                                                                                                     <  7:560500/149‑1 (MQ=255)

                                                                                                                                                                                                      aTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                                                                                                                                                                              >  5:397311/1‑148 (MQ=255)

                                                                                                                                                                                                       ttttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCa                                                                                                                                                                            >  7:419616/1‑149 (MQ=255)

                                                                                                                                                                                                        tttccgtgggcgtAGTTTTGCTCCTGTGTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATGTTAGGGATCCATACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGAGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGgcgc                                             <  4:567588/263‑1 (MQ=255)

                                                                                                                                                                                                            cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg                                                                                                                                                                        <  6:199175/148‑1 (MQ=255)

                                                                                                                                                                                                            cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGa                                                                                                                                                                       <  7:108572/149‑1 (MQ=255)

                                                                                                                                                                                                             gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGc                                                                                                                                        >  2:425796/1‑179 (MQ=255)

                                                                                                                                                                                                                 cgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAg                                                                                                                                    <  4:425794/179‑1 (MQ=255)

                                                                                                                                                                                                                    tatttttGCTCATGATTTCTCCTTTATTACAGCGGGAAAAAGTCTTGCTCATTATATGGAGCCATACTTTTATTACCAACCGCAATTAATTGGCCCGCACCACATTAAGGGAAAACTCGTTTTCACCGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACccac                                  <  1:156420/271‑1 (MQ=255)

                                                                                                                                                                                                                       ttttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATc                                                                                                                                                             >  6:100551/1‑148 (MQ=255)

                                                                                                                                                                                                                         ttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAtt                                                                                                                                                          <  7:16807/149‑1 (MQ=255)

                                                                                                                                                                                                                         ttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAt                                                                                                                                                           <  5:849480/148‑1 (MQ=255)

                                                                                                                                                                                                                          tGCTCATGATTTTTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCACCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAGGATGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCa                           <  2:346430/275‑1 (MQ=255)

                                                                                                                                                                                                                          tGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAtt                                                                                                                                                          <  6:115800/148‑1 (MQ=255)

                                                                                                                                                                                                                             tCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACCGTATTTGCATTTGAGCGCGATATCATTGgc                                                                                                                                                       >  5:590600/1‑148 (MQ=255)

                                                                                                                                                                                                                              cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGgcgc                                                                                                                                                     <  7:153683/149‑1 (MQ=255)

                                                                                                                                                                                                                                  aTTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtt                                                                                                                                                  >  6:344368/1‑148 (MQ=255)

                                                                                                                                                                                                                                   tttCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtttt                                                                                                                                                <  8:876160/149‑1 (MQ=255)

                                                                                                                                                                                                                                      ctcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTtta                                                                                                                                                                                                                                            >  4:482946/1‑54 (MQ=255)

                                                                                                                                                                                                                                       tcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGc                                                                                                                                             <  6:618491/148‑1 (MQ=255)

                                                                                                                                                                                                                                          tttATTACAGCGGTAAAAAGTCTTGCTCATTATAGGGAGCCAGACTTTTATTAc                                                                                                                                                                                                                                        <  2:482948/54‑1 (MQ=255)

                                                                                                                                                                                                                                           ttattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                                                                                                                         >  5:493221/1‑148 (MQ=255)

                                                                                                                                                                                                                                                aCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGc                                                                                                                                   <  7:745770/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGTGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGCTGCGCTCCGGCCGACTTGGGTGCTCTTTACCCACCATTTCATAGTCGTCGATCCGTTTAACCGTgg  >  2:555160/1‑275 (MQ=255)

                                                                                                                                                                                                                                                   gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcaca                                                                                                                                                                                     >  2:138913/1‑96 (MQ=255)

                                                                                                                                                                                                                                                   gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAt                                                                                                              >  3:218053/1‑167 (MQ=255)

                                                                                                                                                                                                                                                   gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGc                                                                                      >  1:108848/1‑191 (MQ=255)

                                                                                                                                                                                                                                                   gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTgg  >  2:448696/1‑275 (MQ=255)

                                                                                                                                                                                                                                                   gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGACCGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTTTTCGGGCAGACCAGCACATTGTCGATGCGCCCCGGCAGACGTGGGCGCGATTTACCCCCCCCTTACTAGTTGTCGATACGGTTAACCGtc   >  4:144305/1‑273 (MQ=255)

                                                                                                                                                                                                                                                     ggTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGAt                                                                                                                              <  8:255940/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       tAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAt                                                                                                                                                                                 <  4:138912/96‑1 (MQ=255)

                                                                                                                                                                                                                                                       tAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCt                                                                                                          <  1:218054/167‑1 (MQ=255)

                                                                                                                                                                                                                                                         aaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGaa                                                                                                                           <  6:218983/148‑1 (MQ=255)

                                                                                                                                                                                                                                                         aaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGaa                                                                                                                           <  6:493221/148‑1 (MQ=255)

                                                                                                                                                                                                                                                          aaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGaaa                                                                                                                          <  5:37113/148‑1 (MQ=255)

                                                                                                                                                                                                                                                          aaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGaaa                                                                                                                          >  6:887584/1‑148 (MQ=255)

                                                                                                                                                                                                                                                           aaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAAcc                                                                                                                        >  7:495050/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              tCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                                                                                                                      >  5:107337/1‑148 (MQ=255)

                                                                                                                                                                                                                                                              tCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGt                                                                                                                     >  8:340420/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               cTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGAATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGt                                                                                                                     >  5:76325/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                    tCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTggg                                                <  2:588482/212‑1 (MQ=255)

                                                                                                                                                                                                                                                                     cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCa                                                                                                               >  6:714001/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                       ttATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCa                                                                                                               <  6:610613/146‑1 (MQ=255)

GGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGG  >  NC_000913/4470689‑4471206

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