BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,880,852	T→A	6.1%	intergenic (‑93/+34)	leuE ← / ← dmlR	leucine efflux protein/transcriptional activator of dmlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,880,852	0	T	A	6.1%	641.4 / 20.0	183	intergenic (‑93/+34)	leuE/dmlR	leucine efflux protein/transcriptional activator of dmlA
Reads supporting (aligned to +/- strand): ref base T (101/70); new base A (6/6); total (107/76)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 187 total aligned reads displayed.

CCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGTACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGA  >  NC_000913/1880602‑1881108

                                 ccTTTCATACCGCTACTGACGCTATTTTTGAGTACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAAACAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTAGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAACCGATAATTGAATCATTAATCAGAAATGAGCCAGTGAAAAAt                                                                                                                                                                                                          >  1:63416/1‑274 (MQ=255)

                                            tcTACTGACGCTATTTTTGAGTACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCt                                                                                                                                                                                              <  1:47100/274‑1 (MQ=255)

                                                         tttttGAGTACAAACAGGGAATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTg                                                                                                                                                                                 <  2:79402/275‑1 (MQ=255)

                                                                                      tGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACt                                                                                                                                                     >  2:9005/1‑274 (MQ=255)

                                                                                                            gCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATCAACCTCTTTAATCATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGCGAAACTCAGNTTTTAACCTGATATCAACCCGACAATTGAATCATTAATCAGGCATGCTCCAGTGCAAAATTCGGAGAGTGCCCGNCAAAATACTGGCGAAAAAATATTACCGTAANNCGATCTTTCGCTGACGTCGCCAGcatg                                                                                                                                >  3:42116/1‑270 (MQ=255)

                                                                                                             ccccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGa                                                                                                                                                                                                                                                           >  6:1471822/1‑149 (MQ=255)

                                                                                                              cccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGat                                                                                                                                                                                                                                                          >  8:538167/1‑149 (MQ=255)

                                                                                                              cccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGat                                                                                                                                                                                                                                                          <  6:1261077/149‑1 (MQ=255)

                                                                                                               ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGGATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGatt                                                                                                                                                                                                                                                         <  7:1260614/149‑2 (MQ=255)

                                                                                                               ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGata                                                                                                                                                                                                                                                         >  6:893291/1‑149 (MQ=255)

                                                                                                                 aaCCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATc                                                                                                                                                                                                                                                       >  5:180141/1‑149 (MQ=255)

                                                                                                                     aGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACGTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACGTGatttagatc                                                                                                                                                                                                                                                   >  7:1009869/1‑142 (MQ=255)

                                                                                                                     aGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAcc                                                                                                                                                                                                                                                   >  7:269154/1‑149 (MQ=255)

                                                                                                                     aGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAcc                                                                                                                                                                                                                                                   >  6:1132085/1‑149 (MQ=255)

                                                                                                                      gATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAccc                                                                                                                                                                                                                                                  <  8:1175922/149‑1 (MQ=255)

                                                                                                                         aGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAc                                                                                                                                                                                                                                                    >  7:200757/1‑144 (MQ=255)

                                                                                                                         aGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGAt                                                                                                                                                                                                                                               <  7:935984/149‑1 (MQ=255)

                                                                                                                            tCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAAt                                                                                                                                                                                                                                            >  8:1002437/1‑149 (MQ=255)

                                                                                                                            tCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAAt                                                                                                                                                                                                                                            >  5:529764/1‑149 (MQ=255)

                                                                                                                            tCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAAt                                                                                                                                                                                                                                            >  5:327269/1‑149 (MQ=255)

                                                                                                                              cAGCAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTg                                                                                                                                                                                                                                          >  8:924843/1‑149 (MQ=255)

                                                                                                                               aGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAAt                                                                                                                                                                                                                                <  4:129278/158‑1 (MQ=255)

                                                                                                                                gTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGaa                                                                                                                                                                                                                                        <  7:358205/149‑1 (MQ=255)

                                                                                                                                   atccacaaCCCCGTATTCAACAAACACATNAACCTTTTTAAGTATATTGCGATCTNAGCCAATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCNATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAGATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCg                                                                                                       <  3:63416/269‑1 (MQ=255)

                                                                                                                                    ttCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCAt                                                                                                                                                                                                                                    >  6:279044/1‑149 (MQ=255)

                                                                                                                                    ttCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCAt                                                                                                                                                                                                                                    >  5:279538/1‑149 (MQ=255)

                                                                                                                                    ttCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCAt                                                                                                                                                                                                                                    >  6:636852/1‑149 (MQ=255)

                                                                                                                                      cAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTa                                                                                                                                                                                                                                  <  5:1132246/149‑1 (MQ=255)

                                                                                                                                      cAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTa                                                                                                                                                                                                                                  >  7:765722/1‑149 (MQ=255)

                                                                                                                                      cAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTa                                                                                                                                                                                                                                  >  5:695225/1‑149 (MQ=255)

                                                                                                                                        gAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAAt                                                                                                                                                                                                                                >  8:396522/1‑149 (MQ=255)

                                                                                                                                              cGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGa                                                                                                                                                                                                                                                           >  4:109256/1‑116 (MQ=255)

                                                                                                                                                  ttCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATc                                                                                                                                                                                                                                                       <  2:109256/116‑1 (MQ=255)

                                                                                                                                                      gCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGa                                                                                                                                                                                                                                                           <  5:1015042/108‑1 (MQ=255)

                                                                                                                                                      gCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGa                                                                                                                                                                                                                                                           >  6:1015042/1‑108 (MQ=255)

                                                                                                                                                         aaCACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGa                                                                                                                                                                                                               >  8:878394/1‑149 (MQ=255)

                                                                                                                                                         aaCACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGa                                                                                                                                                                                                               >  5:458248/1‑149 (MQ=255)

                                                                                                                                                          acacATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGaa                                                                                                                                                                                                              <  6:725552/149‑1 (MQ=255)

                                                                                                                                                            acaTTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGaaaa                                                                                                                                                                                                            >  6:80742/1‑149 (MQ=255)

                                                                                                                                                             caTTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGaaaaa                                                                                                                                                                                                           >  8:210822/1‑149 (MQ=255)

                                                                                                                                                             caTTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGaaaaa                                                                                                                                                                                                           <  8:269154/149‑1 (MQ=255)

                                                                                                                                                              aTTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAAt                                                                                                                                                                                                          >  7:736530/1‑149 (MQ=255)

                                                                                                                                                              aTTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAAt                                                                                                                                                                                                          >  6:1468362/1‑149 (MQ=255)

                                                                                                                                                                  aCCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCgg                                                                                                                                                                                                      <  6:607401/149‑1 (MQ=255)

                                                                                                                                                                    ctctTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGt                                                                                                                                                                                                    >  6:1066614/1‑149 (MQ=255)

                                                                                                                                                                    ctctTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGt                                                                                                                                                                                                    >  8:1455414/1‑149 (MQ=255)

                                                                                                                                                                      ctTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGATAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTTAAACATAAACGGCACAGACGTTCGCTGGCTTGTAATATTCCGGTAGCCATTGCAcgcaa                                                                    >  2:216322/1‑271 (MQ=255)

                                                                                                                                                                       tttAATTATCTTTCGATCATGCGCGCGCAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGCTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAGTTCGGGTAGt                                                                                                                                                                                                 >  8:1104829/1‑149 (MQ=255)

                                                                                                                                                                       tttAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATc                                                                                                                                                                                                                                      <  5:890227/112‑1 (MQ=255)

                                                                                                                                                                       tttAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGt                                                                                                                                                                                                 >  5:1100827/1‑149 (MQ=255)

                                                                                                                                                                       tttAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGt                                                                                                                                                                                                 >  6:485456/1‑149 (MQ=255)

                                                                                                                                                                           aTTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTc                                                                                                                                                                                             >  8:71718/1‑149 (MQ=255)

                                                                                                                                                                            ttATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCg                                                                                                                                                                                            >  5:155225/1‑149 (MQ=255)

                                                                                                                                                                                 tttCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAg                                                                                                                                                                                                                             >  5:921302/1‑111 (MQ=255)

                                                                                                                                                                                 tttCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAg                                                                                                                                                                                                                             <  6:921302/111‑1 (MQ=255)

                                                                                                                                                                                     gaccATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAAGAc                                                                                                                                                                                   >  7:970082/4‑149 (MQ=255)

                                                                                                                                                                                       tCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTg                                                                                                                                                                                 <  7:1294019/149‑1 (MQ=255)

                                                                                                                                                                                       tCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTg                                                                                                                                                                                 >  8:84064/1‑149 (MQ=255)

                                                                                                                                                                                        cATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTgg                                                                                                                                                                                <  7:396522/149‑1 (MQ=255)

                                                                                                                                                                                         aTGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGc                                                                                                                                                                               <  7:217897/149‑1 (MQ=255)

                                                                                                                                                                                             gcgATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTaa                                                                                                                                                                           <  8:799274/149‑1 (MQ=255)

                                                                                                                                                                                                 ttAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAAt                                                                                                                                                                       >  6:1430225/1‑149 (MQ=255)

                                                                                                                                                                                                 ttAAAGGTAAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAAt                                                                                                                                                                       >  8:126009/1‑149 (MQ=255)

                                                                                                                                                                                                  tAAAGGTGAAATTGGTAACCAACCTGTAGCGGGTTAAAAAGGAGAAAATCAGGGTTTAACCGGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACAt                                        <  4:9005/275‑1 (MQ=255)

                                                                                                                                                                                                            tatGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGc                                                                                                                                                                               <  7:1016610/130‑1 (MQ=255)

                                                                                                                                                                                                            tatGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGc                                                                                                                                                                               >  8:1016610/1‑130 (MQ=255)

                                                                                                                                                                                                               gCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCg                                                                                                                                                         <  7:71718/149‑1 (MQ=255)

                                                                                                                                                                                                                cTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGc                                                                                                                                                        >  5:1047731/1‑149 (MQ=255)

                                                                                                                                                                                                                 tAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCa                                                                                                                                                       <  8:1310785/149‑1 (MQ=255)

                                                                                                                                                                                                                 tAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCa                                                                                                                                                       >  8:430489/1‑149 (MQ=255)

                                                                                                                                                                                                                 tAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCa                                                                                                                                                       >  6:1109229/1‑149 (MQ=255)

                                                                                                                                                                                                                  aaCCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACa                                                                                                                          >  4:121220/1‑177 (MQ=255)

                                                                                                                                                                                                                     cAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACttt                                                                                                                                                   >  6:1347512/1‑149 (MQ=255)

                                                                                                                                                                                                                     cAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACttt                                                                                                                                                   >  5:387261/1‑149 (MQ=255)

                                                                                                                                                                                                                      aaTCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATaaa                                                                                                                      <  2:121220/177‑1 (MQ=255)

                                                                                                                                                                                                                       aTCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGaa                                                                                                                                                                                                                                        <  7:1361220/62‑1 (MQ=255)

                                                                                                                                                                                                                       aTCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGaa                                                                                                                                                                                                                                        >  8:1361220/1‑62 (MQ=255)

                                                                                                                                                                                                                        tCTGTAGCGGCTTAGAAAGGAGAAAATTAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGc                                                                                                                                                <  6:1110632/149‑1 (MQ=255)

                                                                                                                                                                                                                        tCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGTATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGc                                                                                                                                                <  8:549057/149‑1 (MQ=255)

                                                                                                                                                                                                                        tCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCAATTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAAt                  >  2:157907/1‑275 (MQ=255)

                                                                                                                                                                                                                         cTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCt                                                                                                                                               <  5:1471822/149‑1 (MQ=255)

                                                                                                                                                                                                                                   ttAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCa                                                                                                                                     >  6:416565/1‑149 (MQ=255)

                                                                                                                                                                                                                                    tAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCgc                                                                                                      >  4:327162/1‑179 (MQ=255)

                                                                                                                                                                                                                                     aGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGc                                                                                                                                   >  6:292560/1‑149 (MQ=255)

                                                                                                                                                                                                                                      gAAAGATGAAAAGCAGGTTTTAACCGTATAACAACCCGAAAATTGAATCATTAATCAGGCAGGCTCCAGTTAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCga    <  1:36486/275‑1 (MQ=255)

                                                                                                                                                                                                                                        aaGGAGAAATCCAGGTTTTAACTGTATATCAACCCGATAATTGAATCATTAAACAGGCATGCGCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCACCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCgaga  <  2:23284/275‑1 (MQ=255)

                                                                                                                                                                                                                                        aaGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCgaga  <  1:16060/275‑1 (MQ=255)

                                                                                                                                                                                                                                           gagaAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGaa                                                                                                                             >  5:1441176/1‑149 (MQ=255)

                                                                                                                                                                                                                                            agaAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGc                                                                                                                                                >  4:364348/1‑129 (MQ=255)

                                                                                                                                                                                                                                            agaAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGaaa                                                                                                                            <  8:736530/149‑1 (MQ=255)

                                                                                                                                                                                                                                             gaAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAAc                                                                                                                           <  8:416469/149‑1 (MQ=255)

                                                                                                                                                                                                                                               aaaTCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCAGACGTCGCCAGCCTTGAAACAt                                                                                                                         <  6:180141/149‑1 (MQ=255)

                                                                                                                                                                                                                                                aaTCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGAc                                                                                                                                            <  2:364349/129‑1 (MQ=255)

                                                                                                                                                                                                                                                 aTCAGGTTAAAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATaa                                                                                                                       >  6:1177008/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 aTCAGGTCTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATaa                                                                                                                       >  6:302434/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  tCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATaaa                                                                                                                      >  6:1164155/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   cAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAAc                                                                                                                     >  7:559666/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      gTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGc                                                                                                                  <  6:517135/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      gTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTAGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGc                                                                                                                  >  8:341821/1‑149 (MQ=255)

                                                                                                                                                                                                                                                         ttAACCTGATATCAACCCGATAATGGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGTTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCa                                                                                                               <  7:894120/149‑1 (MQ=255)

CCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGTACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGA  >  NC_000913/1880602‑1881108

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