BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,425,971	C→A	5.1%	intergenic (‑55/+35)	hisQ ← / ← hisJ	histidine ABC transporter permease/histidine ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,425,971	0	C	A	5.1%	693.4 / 20.2	198	intergenic (‑55/+35)	hisQ/hisJ	histidine ABC transporter permease/histidine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand): ref base C (102/86); new base A (7/3); total (109/89)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-02

Only 100 of 201 total aligned reads displayed.

CCGGCACGCCACGAATCAGCGTGGTGTACCCTTCGAAAATCAGCCCCGAAAGCCGATTTTGCGAGAGCTTACCGCCAGCGCCAATTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCA  >  NC_000913/2425736‑2426240

ccGGCACGCCACGAATCAGCGTGGTGTACCCTTCGAAAATCAGCCCCGAAAGCCGATTTTGCGAGAGCTTACCGCCAGCGCCAATTAAACCGATGATTACAGCGAGCACTACAGCGCTGATAGCCAGCCCCAGCGTGACGAGCGCACCCTGTAAACTAACACCTGAAAACCCATCCAACATGCTGCCTGTCCTGTCGTGTGTGGCTAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGGAAGGCGCGGGCTTTCTTTCACCGATCGCCATTAGc                                                                                                                                                                                                                                        >  4:219041/1‑275 (MQ=255)

                                           ccccGAAAGCCGATTTTGCGAGAGCTTACCGCCAGCGCCAATTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTCGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTCGCCACCATAAACATCAAAATCGAAGTACCTTTTCGCTACTTtct                                                                                                                                                                                              >  1:482/1‑274 (MQ=255)

                                                                                           gatgatTACAGCGAGCACTACAGAGCTGATAGCCAGCTACAGCGTGACGAGCGCACCCTGTAAAATGACACCTGAAAACCCATACAACATGCTGCCTGCCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGa                                                                                                                                                                                                                                                                           >  7:906724/1‑149 (MQ=255)

                                                                                            atgatTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGaa                                                                                                                                                                                                                                                                          <  7:1374433/149‑1 (MQ=255)

                                                                                            atgatTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGaa                                                                                                                                                                                                                                                                          <  6:105779/149‑1 (MQ=255)

                                                                                             tgatTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAg                                                                                                                                                                                                                                                                         <  8:961659/149‑1 (MQ=255)

                                                                                                  aCAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGc                                                                                                                                                                                                                                                                    <  8:776769/149‑1 (MQ=255)

                                                                                                    aGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCgg                                                                                                                                                                                                                                                                  <  5:954770/149‑1 (MQ=255)

                                                                                                         gCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGcttt                                                                                                                                                                                                                                                             >  8:371062/1‑149 (MQ=255)

                                                                                                          cACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTAAATTACCGCCTCTTGTCTCCCTCCGGTTNCCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGCTTAGCCACCATAAACATCAAAATCGAAGTACTTTNTCGCTAATTTCCCGTAAGTAACGTCAGCGCGCATTTCGGCAAAGGNTNTTTTAAGTGCTTCGAGCAGTTCGtttt                                                                                                                               >  3:34126/1‑272 (MQ=255)

                                                                                                           aCTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttct                                                                                                                                                                                                                                                           <  5:17549/149‑1 (MQ=255)

                                                                                                           aCTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttct                                                                                                                                                                                                                                                           <  6:413471/149‑1 (MQ=255)

                                                                                                           aCTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttct                                                                                                                                                                                                                                                           <  7:511808/149‑1 (MQ=255)

                                                                                                            cTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttctt                                                                                                                                                                                                                                                          >  6:1480910/1‑149 (MQ=255)

                                                                                                             tACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttcttt                                                                                                                                                                                                                                                         >  5:77981/1‑149 (MQ=255)

                                                                                                             tACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttcttt                                                                                                                                                                                                                                                         <  7:752773/149‑1 (MQ=255)

                                                                                                                agagCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCAc                                                                                                                                                                                                                                                      >  5:797078/1‑149 (MQ=255)

                                                                                                                     tGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATg                                                                                                                                                                                                                                                 <  8:1147275/149‑1 (MQ=255)

                                                                                                                     tGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATg                                                                                                                                                                                                                                                 <  5:1382921/149‑1 (MQ=255)

                                                                                                                       aTAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGc                                                                                                                                                                                                                                               >  5:77803/1‑149 (MQ=255)

                                                                                                                       aTAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAACTCGAAGTAATTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCCCATTTCGTCACAGGCTTTGTTCAGTGCTTCGCACCGTTCGTTCTCTTCTTTACGCAg                                                                                                                  >  4:81960/1‑274 (MQ=255)

                                                                                                                       aTAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGATTTCTTTCACCGATTGCCACTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTCCCGTCAGCGCGCATTTCGGCAACGGCTTTGTTCATTGCTTCGAGCAGTTCTTTATCTTCTTTACGCAg                                                                                                                  >  2:283109/1‑274 (MQ=255)

                                                                                                                            cAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTa                                                                                                                                                                                                                                          >  5:1230335/1‑149 (MQ=255)

                                                                                                                            cAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTa                                                                                                                                                                                                                                          >  7:554925/1‑149 (MQ=255)

                                                                                                                            cAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTa                                                                                                                                                                                                                                          <  6:304009/149‑1 (MQ=255)

                                                                                                                                tCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGcca                                                                                                                                                                                                                                      >  7:211612/1‑149 (MQ=255)

                                                                                                                                   aGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGccacca                                                                                                                                                                                                                                   <  5:836074/149‑1 (MQ=255)

                                                                                                                                       tGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTa                                                                                                 >  2:62736/1‑275 (MQ=255)

                                                                                                                                         aCGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACa                                                                                                                                                                                                                             <  6:1422872/149‑1 (MQ=255)

                                                                                                                                         aCGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACa                                                                                                                                                                                                                             >  5:597279/1‑149 (MQ=255)

                                                                                                                                           gAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATc                                                                                                                                                                                                                           <  5:240916/149‑1 (MQ=255)

                                                                                                                                                cACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAAt                                                                                                                                                                                                                      <  8:54560/149‑1 (MQ=255)

                                                                                                                                                   ccTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGa                                                                                                                                                                                                                   >  6:750083/1‑149 (MQ=255)

                                                                                                                                                    cTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                                  <  6:77981/149‑1 (MQ=255)

                                                                                                                                                    cTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                                  <  7:459589/149‑1 (MQ=255)

                                                                                                                                                     tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAg                                                                                                                                                                                                                 >  5:586489/1‑149 (MQ=255)

                                                                                                                                                      gTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGt                                                                                                                                                                                                                >  7:1357297/1‑149 (MQ=255)

                                                                                                                                                        aaaaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTAc                                                                                                                                                                                                              <  8:144473/149‑1 (MQ=255)

                                                                                                                                                          aaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACtt                                                                                                                                                                                                            <  7:534173/149‑1 (MQ=255)

                                                                                                                                                          aaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACtt                                                                                                                                                                                                            >  6:1017022/1‑149 (MQ=255)

                                                                                                                                                            tAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACtttt                                                                                                                                                                                                          >  7:612986/1‑149 (MQ=255)

                                                                                                                                                               cacCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCCAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACAGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGTCCCATTCCGGTCCCTACCCCAAACCGTTTTTCATCtttc                                                                          >  4:232286/1‑273 (MQ=255)

                                                                                                                                                                acCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGc                                                                                                                                                                                                      <  6:1081809/149‑1 (MQ=255)

                                                                                                                                                                 ccTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGAtt                                                                                                                                                                                                                                           >  5:662315/1‑111 (MQ=255)

                                                                                                                                                                 ccTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGAtt                                                                                                                                                                                                                                           <  6:662315/111‑1 (MQ=255)

                                                                                                                                                                  cTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAggg                                                                                                                                                                                                                                                                       <  6:262356/82‑1 (MQ=255)

                                                                                                                                                                  cTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                                  >  5:1293783/1‑135 (MQ=255)

                                                                                                                                                                        aCCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATg                                                                                                                                                                                                                                                 >  8:58211/1‑98 (MQ=255)

                                                                                                                                                                        aCCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATg                                                                                                                                                                                                                                                 <  7:58211/98‑1 (MQ=255)

                                                                                                                                                                        aCCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtct                                                                                                                                                                                              >  8:1213635/1‑149 (MQ=255)

                                                                                                                                                                         cccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                                                                                                                                                             <  6:1261873/149‑1 (MQ=255)

                                                                                                                                                                          ccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGAt                                                                                                                                                                                                                                                  <  6:933062/95‑1 (MQ=255)

                                                                                                                                                                           cATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                                  <  6:472356/126‑1 (MQ=255)

                                                                                                                                                                           cATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                                  >  5:472356/1‑126 (MQ=255)

                                                                                                                                                                           cATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGt                                                                                                                                                                                           >  8:23798/1‑149 (MQ=255)

                                                                                                                                                                             tACAACCTGCTGCCTGTCCTGTCGTGTGTGGTAAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCg                                                           <  1:34126/275‑1 (MQ=255)

                                                                                                                                                                                      ctgcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCa                                                                                                                                                                                <  5:717994/149‑1 (MQ=255)

                                                                                                                                                                                        gcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTGTCATCTTTACCGAACGGGCCCCCGGATTTggactc                                                >  2:305253/1‑270 (MQ=255)

                                                                                                                                                                                        gcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATC‑AAATCGaa                                                                                                                                                                                                                  >  7:1154481/1‑112 (MQ=255)

                                                                                                                                                                                            gtcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgcgc                                                                                                                                                                          <  8:1357297/149‑1 (MQ=255)

                                                                                                                                                                                             tcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCa                                                                                                                                                                         >  8:70163/1‑149 (MQ=255)

                                                                                                                                                                                             tcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCAGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCa                                                                                                                                                                         >  5:1165006/1‑149 (MQ=255)

                                                                                                                                                                                               ctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCAtt                                                                                                                                                                       >  5:90065/1‑149 (MQ=255)

                                                                                                                                                                                               ctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCAtt                                                                                                                                                                       >  7:1417664/1‑149 (MQ=255)

                                                                                                                                                                                                 gtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTc                                                                                                                                                                     <  7:99879/149‑1 (MQ=255)

                                                                                                                                                                                                      gtgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCaa                                                                                                                                                                >  5:261981/1‑149 (MQ=255)

                                                                                                                                                                                                       tgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCaaa                                                                                                                                                               <  6:597279/149‑1 (MQ=255)

                                                                                                                                                                                                       tgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGGATTCTTTACCGACAGtttg                                 >  4:109772/1‑275 (MQ=255)

                                                                                                                                                                                                       tgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCAAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCCAAAGgtt                                  >  4:302405/1‑274 (MQ=255)

                                                                                                                                                                                                       tgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCAGGCCGGTCCCTACGCCAAACCGTTTTTCATCTTTAAACACCGGGCCACGGACTTTTGAATCtttcacaaaaggtt                                  >  4:243575/1‑263 (MQ=255)

                                                                                                                                                                                                        gtggtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAg                                                                                                                                                              <  6:1298876/149‑1 (MQ=255)

                                                                                                                                                                                                           gtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAg                                                                                                                                                              >  7:1139083/1‑146 (MQ=255)

                                                                                                                                                                                                             gAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCttt                                                                                                                                                         >  8:285676/1‑149 (MQ=255)

                                                                                                                                                                                                             gAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCttt                                                                                                                                                         >  6:119812/1‑149 (MQ=255)

                                                                                                                                                                                                              aaTTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCgg                                                                                                                                                                   <  5:647611/138‑1 (MQ=255)

                                                                                                                                                                                                                ttACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGtt                                                                                                                                                      >  5:1313981/1‑149 (MQ=255)

                                                                                                                                                                                                                 tACCGCCTCTTGGCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACAATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGGACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTGAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGaaa                       <  2:219041/275‑1 (MQ=255)

                                                                                                                                                                                                                   ccGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAg                                                                                                                                                   <  8:670435/149‑1 (MQ=255)

                                                                                                                                                                                                                    cGCCTCTTGTCTCCCTCCGGTTACCCGTAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                                                                                                                                                             <  8:1049388/106‑1 (MQ=255)

                                                                                                                                                                                                                    cGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                                                                                                                                                             >  7:1049388/1‑106 (MQ=255)

                                                                                                                                                                                                                       tcgtttttcccccccGTTTACCNNAAGGTNGTGGGCTTTTTTTCACCGTTGGTGATTATCCGCCATAAACAGCAGAATCTAAGGACTTTTTCGCTAATTTCTCGTAAGTCCCTTCATCGCGCTTTTCGGCAAGCGCTTTGTTCAGTGCTTCCCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCg                 <  3:482/266‑1 (MQ=255)

                                                                                                                                                                                                                       ctctTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCt                                                                                                                                               >  7:1010720/1‑149 (MQ=255)

                                                                                                                                                                                                                         ctTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTc                                                                                                                                             >  6:812818/1‑149 (MQ=255)

                                                                                                                                                                                                                           ccgctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCg                                                                                                                                                                                            <  7:567796/97‑1 (MQ=38)

                                                                                                                                                                                                                           ccgctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCg                                                                                                                                                                                            >  8:567796/4‑100 (MQ=38)

                                                                                                                                                                                                                            cgctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTcgcg                                                                                                                                          >  7:1203352/3‑149 (MQ=255)

                                                                                                                                                                                                                              ctccctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCa                                                                                                                                        >  6:619286/1‑149 (MQ=255)

                                                                                                                                                                                                                               tccctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCtt                                                                                                                                              >  8:109516/1‑142 (MQ=255)

                                                                                                                                                                                                                                 ccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                                                                                                                                                             >  5:271023/1‑93 (MQ=255)

                                                                                                                                                                                                                                 cctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGtt                                                                                                                                     >  7:482250/1‑149 (MQ=255)

                                                                                                                                                                                                                                 cctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGtt                                                                                                                                     >  6:13377/1‑149 (MQ=255)

                                                                                                                                                                                                                                   tccGGTTACCCGAAGGGAGCGGGCTTTGTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCg                                                                                                                                   >  7:975494/1‑149 (MQ=255)

                                                                                                                                                                                                                                   tccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGAcc     <  4:62736/275‑1 (MQ=255)

                                                                                                                                                                                                                                   tccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGAcc     <  2:109772/275‑1 (MQ=255)

                                                                                                                                                                                                                                    ccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGt                                                                                                                                  <  6:84119/149‑1 (MQ=255)

                                                                                                                                                                                                                                       gTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATACACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCCTAAGTACCGTCAGCGCGCATTTCGGCAATGGCTTTGTTCAGTGCTTCGCGCTGTTCGTTATCTTCTTTACGCATGCCCATTCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACATCCAGGCCCCCGAAGTTGTACTCTTTACAGACCGGTTGTTTGAGGCAACATTCCCTGGcgaccacctga  >  4:185861/1‑264 (MQ=255)

                                                                                                                                                                                                                                          aCCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtctt                                                                                                                            <  7:1455695/149‑1 (MQ=255)

                                                                                                                                                                                                                                          aCCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtctt                                                                                                                            >  6:671792/1‑149 (MQ=255)

                                                                                                                                                                                                                                           cccGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtcttc                                                                                                                           >  7:335371/1‑149 (MQ=255)

                                                                                                                                                                                                                                           cccGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAt                                                                                                                               <  2:138145/145‑1 (MQ=255)

CCGGCACGCCACGAATCAGCGTGGTGTACCCTTCGAAAATCAGCCCCGAAAGCCGATTTTGCGAGAGCTTACCGCCAGCGCCAATTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCA  >  NC_000913/2425736‑2426240

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