BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,092,779	Δ1 bp	8.3%	coding (902/951 nt)	gshB →	glutathione synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,092,779	0	T	.	8.3%	349.0 / 25.5	120	coding (902/951 nt)	gshB	glutathione synthetase
Reads supporting (aligned to +/- strand): ref base T (65/44); new base . (4/6); total (70/50)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 126 total aligned reads displayed.

CTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGA  >  NC_000913/3092630‑3092926

cTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                                     >  2:911614/1‑150 (MQ=255)

cTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                                     >  2:465307/1‑150 (MQ=255)

cTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCAACGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                                     >  2:1518053/1‑150 (MQ=255)

       aTCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGa                                                                                                                                              <  2:1331665/150‑1 (MQ=255)

        tCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAACGGGCTGATTTTTGTTGGGCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaa                                                                                                                                             <  1:696797/150‑1 (MQ=255)

          gCCCGTCAGATCGGGCCGACGCTTAAAGAAAAAGTGCTGATTTTTGTTTGTATGGATATAATCTGCGAACGTCTGACTTAAATTAAATTCACCAGCCCAACCTGTATTCGTTAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaac                                                                                                                                           <  1:1321763/150‑3 (MQ=255)

          gCCCGTCAGATCGGGCCGACGCTGAAAGAACAAGGCCTGACTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACTTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaac                                                                                                                                           <  1:1020256/150‑3 (MQ=255)

          gCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGGTTCCGGTGATCGA‑CACCGGaaac                                                                                                                                           <  1:1549067/150‑3 (MQ=255)

          gCCCGGCATATCGGGCAGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaac                                                                                                                                           <  1:333898/150‑3 (MQ=255)

               tCAGAACCGGCCGACGCTGAAAGAACAAGGGCTGATTTTTGATGCTCGGGATATCCTCGGCGACCGTCAGACTGAAATTAGCGTCACCAGCCCAACCTCTATTCTTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAAt                                                                                                                                      <  1:321560/150‑1 (MQ=255)

                cAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATg                                                                                                                                     <  2:53104/150‑1 (MQ=255)

                cAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATg                                                                                                                                     >  2:977970/1‑150 (MQ=255)

                  gATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGAAATCATCGGCGACCGTCTGACTGAAATTAAGGTCACCCGCGCAACCTGTATTCGTGTGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAAAGTTAATGGa                                                                                                                                   >  2:319715/1‑150 (MQ=255)

                       ggCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCa                                                                                                                              >  1:134190/1‑150 (MQ=255)

                               cTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGTATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGTATGCCATCTAAGCa                                                                                                                      >  1:1634661/1‑150 (MQ=255)

                               cTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCa                                                                                                                      >  1:1662936/1‑150 (MQ=255)

                               cTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCa                                                                                                                      >  1:1619967/1‑150 (MQ=255)

                               cTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTCACTGAAATTAACGTCACCAGCCCAACCTGAATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCa                                                                                                                      >  2:968537/1‑150 (MQ=255)

                                 gaaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACg                                                                                                                    >  2:306667/1‑150 (MQ=255)

                                 gaaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCAGAATGTTAATGGATGCCATCGAAGCACg                                                                                                                    >  2:1122037/1‑150 (MQ=255)

                                   aagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGtt                                                                                                                  >  2:5622/1‑150 (MQ=255)

                                    agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGttt                                                                                                                 <  1:454421/150‑1 (MQ=255)

                                    agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGttt                                                                                                                 >  1:1708692/1‑150 (MQ=255)

                                     gaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTa                                                                                                                >  2:73861/1‑150 (MQ=255)

                                      aaaaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAc                                                                                                               >  2:804710/1‑150 (MQ=255)

                                        aaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcag                                                                                                             <  2:1226624/150‑1 (MQ=255)

                                        aaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcag                                                                                                             <  2:1226747/150‑1 (MQ=255)

                                           gggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTTTTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcag                                                                                                          <  1:871905/150‑1 (MQ=255)

                                           gggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcag                                                                                                          >  1:1042963/1‑150 (MQ=255)

                                           gggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcag                                                                                                          >  2:566761/1‑150 (MQ=255)

                                              cTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagcag                                                                                                       >  1:1614512/1‑150 (MQ=255)

                                                    tttGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAAccc                                                                                                 <  1:109336/150‑1 (MQ=255)

                                                         tgGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCtt                                                                                            >  2:1159987/1‑150 (MQ=255)

                                                          ggTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATTGATGCCATCGAAGAACGTTTACAGCAGCAGTAACCCACCTTa                                                                                           >  2:910470/1‑150 (MQ=255)

                                                          ggTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTa                                                                                           >  2:75964/1‑150 (MQ=255)

                                                          ggTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTa                                                                                           >  2:1123089/1‑150 (MQ=255)

                                                               ggATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGAAGCAGTAACCCACCTTAGCgag                                                                                      >  2:337683/1‑150 (MQ=255)

                                                               ggATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGACCGACCTGTATTAGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAACACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTatacat                                                                                      >  1:1474235/1‑145 (MQ=255)

                                                                   atcatcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAgg                                                                                  <  1:1565014/150‑1 (MQ=255)

                                                                     catcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTAGTGAGAGTGAAGAAGAGTTTCCGTAG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTATACAGCAGCAGTAACCCACCTTAGCGAGAAGGAt                                                                                >  2:1113332/1‑150 (MQ=255)

                                                                      atcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAtc                                                                               <  2:1662936/150‑1 (MQ=255)

                                                                        cGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAtctc                                                                             >  2:396698/1‑150 (MQ=255)

                                                                        cGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAtctc                                                                             <  2:27628/150‑1 (MQ=255)

                                                                         ggCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCg                                                                            <  1:1122037/150‑1 (MQ=255)

                                                                           cGACCGTCTGACTGAAATTAACGTCACCAG‑ACAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTg                                                                         <  2:500295/150‑1 (MQ=255)

                                                                            taCCGTCTTACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCATCAGTAACCCACCTTAGCGAGAAGGATCTCGTTg                                                                         <  2:444019/149‑1 (MQ=255)

                                                                            gACCGTGTGACTGAAATTAACGTCACCAGCCCAACCTGTAGTCGTGTGTTTGAAGGAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAAGCCACCTTAGCGAGAAGGATCTCGTTg                                                                         <  2:464734/150‑1 (MQ=255)

                                                                                gTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGAc                                                                     <  2:466191/150‑1 (MQ=255)

                                                                                gTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGAc                                                                     <  1:5622/150‑1 (MQ=255)

                                                                                gTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGAc                                                                     <  1:1433026/150‑1 (MQ=255)

                                                                                  ctgactgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGActc                                                                   >  2:1379073/1‑150 (MQ=255)

                                                                                    gactgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTTATGGATGCCATGGAAGCACGTTTACAGCAGCAGTAACCCACCGTAGCGAGAAGGTTCTCGTTGAGACGCTg                                                                 >  2:285435/1‑150 (MQ=255)

                                                                                     actgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                                <  1:1513487/150‑1 (MQ=255)

                                                                                       tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGt                                                              <  1:461198/150‑1 (MQ=255)

                                                                                         aaaTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGa                                                            <  1:1123089/150‑1 (MQ=255)

                                                                                             tAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAgc                                                        >  2:1626751/1‑150 (MQ=255)

                                                                                               gcGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACCGCAGCAGTAACCCACCTTAGCGAGAAGGATGACGTTGAGACTCTGAGTGACAgcgc                                                      >  1:1065722/2‑150 (MQ=255)

                                                                                                 gTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGccc                                                    >  1:1218132/1‑150 (MQ=255)

                                                                                                  tCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCCCCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCct                                                   <  1:306667/150‑1 (MQ=255)

                                                                                                    aCCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACCGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCcttc                                                 >  1:254066/1‑150 (MQ=255)

                                                                                                     ccAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCcttct                                                >  1:737227/1‑150 (MQ=255)

                                                                                                      cAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCcttctt                                               <  1:1705633/150‑1 (MQ=255)

                                                                                                       aGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCttt                                              >  2:376894/1‑150 (MQ=255)

                                                                                                        gCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGTCAGCGCCCTTCTTTc                                             >  1:1299937/1‑150 (MQ=255)

                                                                                                          ccAACCTGTATTTGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCa                                           <  2:1538874/150‑1 (MQ=255)

                                                                                                            aaCCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACg                                         <  2:162876/150‑1 (MQ=255)

                                                                                                               cTGTTTTCGTGATATTGAAGGAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGT‑ACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATa                                     <  1:1622499/150‑1 (MQ=255)

                                                                                                               cTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCAt                                      <  1:600744/150‑1 (MQ=255)

                                                                                                                 gTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGAGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGAATAc                                    >  2:1523026/1‑150 (MQ=255)

                                                                                                                   aTTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCCCCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTg                                  >  1:1287640/1‑150 (MQ=255)

                                                                                                                      cGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgc                               <  1:114975/150‑1 (MQ=255)

                                                                                                                      cGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgc                               >  1:483347/1‑150 (MQ=255)

                                                                                                                      cGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgc                               >  2:1495587/1‑150 (MQ=255)

                                                                                                                       gTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgcg                              >  1:14491/1‑150 (MQ=255)

                                                                                                                        tGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgcgc                             <  1:566761/150‑1 (MQ=255)

                                                                                                                         gagaTTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCCCCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCTTACTGGgcggg                            >  2:179390/1‑148 (MQ=255)

                                                                                                                         gagaTTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCCCt                            >  2:179386/1‑150 (MQ=255)

                                                                                                                               gAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCt                            <  1:1539841/144‑1 (MQ=255)

                                                                                                                               gAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCt                      >  2:1304301/1‑150 (MQ=255)

                                                                                                                               gAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCt                      <  1:790498/150‑1 (MQ=255)

                                                                                                                               gAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCt                      >  2:1573222/1‑150 (MQ=255)

                                                                                                                               gAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACC‑TAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtt                     >  2:695057/1‑150 (MQ=255)

                                                                                                                                aaGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtt                     >  1:1474640/1‑150 (MQ=255)

                                                                                                                                 aGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCGACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttt                    >  1:880518/1‑150 (MQ=255)

                                                                                                                                  gCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtttt                   >  1:180617/1‑150 (MQ=255)

                                                                                                                                  gCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtttt                   <  2:1042963/150‑1 (MQ=255)

                                                                                                                                   cAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttttt                  >  2:914362/1‑150 (MQ=255)

                                                                                                                                    agagTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTATAAAGCACCAGTAACCCACATTAGCGAGAAGGATATCGTTGAGACTCTGAGTGACCGCTCCCTTCTATACACGCATACTGGGCGCAGTTGCttattt                 >  1:339597/1‑150 (MQ=255)

                                                                                                                                     gagTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTg                <  2:681054/150‑1 (MQ=255)

                                                                                                                                     gagTTTCCGGCG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTg                <  2:1106146/150‑1 (MQ=255)

                                                                                                                                     cagTTTCCGGTG‑GCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAGCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTg                <  2:472935/149‑1 (MQ=255)

                                                                                                                                      agTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGa               >  1:603290/1‑150 (MQ=255)

                                                                                                                                      acattCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCCCCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTGTCCCCGCATACTGGGCGCTGGTGCTTTTTTGa               >  1:347502/4‑150 (MQ=255)

                                                                                                                                       cattCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa              >  2:979165/3‑150 (MQ=255)

                                                                                                                                        attCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAAc             >  1:768515/2‑150 (MQ=255)

                                                                                                                                         ttCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAAcc            >  2:363400/1‑150 (MQ=255)

                                                                                                                                          tCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCa           >  1:306061/1‑150 (MQ=255)

                                                                                                                                             ggTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATAATGGGCGCTGTTGAGTTTTTGCACCAGGa        >  1:495455/1‑150 (MQ=255)

                                                                                                                                                    gATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGa  >  2:1528457/1‑150 (MQ=255)

                                                                                                                                                    gATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGa  >  2:1199034/1‑150 (MQ=255)

CTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGA  >  NC_000913/3092630‑3092926

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