BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,718,281	C→T	6.2%	D324D (GAC→GAT) ‡	ghrB →	glyoxylate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,718,281	0	C	T	6.2%	506.7 / 13.9	160	D324D (GAC→GAT) ‡	ghrB	glyoxylate reductase
Reads supporting (aligned to +/- strand): ref base C (77/72); new base T (7/3); total (85/75)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

Only 100 of 164 total aligned reads displayed.

CTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTG  >  NC_000913/3718132‑3718429

ctctCAATGGCCAACGTCGTCGCAGTACCGCATTTTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGAc                                                                                                                                                      <  2:610199/150‑1 (MQ=255)

ctctCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGAc                                                                                                                                                      <  2:1672594/150‑1 (MQ=255)

ctctCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGAc                                                                                                                                                      >  2:711026/1‑150 (MQ=255)

ctctCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCACACGTCGCGGAc                                                                                                                                                      <  2:1582855/150‑1 (MQ=255)

 tctcAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACt                                                                                                                                                     <  2:321717/150‑1 (MQ=255)

   tcAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTaa                                                                                                                                                   >  2:1383189/1‑150 (MQ=255)

     aaTGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGc                                                                                                                                                 <  1:491251/150‑1 (MQ=255)

     aaTGGCCAACGTCGTCGCAGTACCACATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGc                                                                                                                                                 >  2:352482/1‑150 (MQ=255)

        ggCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCcgc                                                                                                                                              >  2:575788/1‑150 (MQ=255)

          ccAACGGCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACCCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGTAAAGGTTGAGAAGATCTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGa                                                                                                                                            <  2:1265993/150‑1 (MQ=255)

           cAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGAc                                                                                                                                           >  1:320260/1‑150 (MQ=255)

              cgtcgtCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGCTTAGTCCGCGAcgtg                                                                                                                                        <  2:387802/150‑4 (MQ=255)

                tcgtcgCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGt                                                                                                                                      <  2:1756838/150‑1 (MQ=255)

                tcgtcgCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGt                                                                                                                                      >  1:1930518/1‑150 (MQ=255)

                 cgtcgCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTa                                                                                                                                     >  2:724225/1‑150 (MQ=255)

                  gtcgCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTag                                                                                                                                    <  1:2155864/150‑1 (MQ=255)

                  gtcgCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTag                                                                                                                                    >  2:866042/1‑150 (MQ=255)

                         tACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAg                                                                                                                             <  2:1225585/150‑1 (MQ=255)

                          aCCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGc                                                                                                                            >  2:1200333/1‑150 (MQ=255)

                           cctcATCTTGGATCTGCCACCCATGAGAGTCGTTATCGCATGGCCGCCTGTGCCCGGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTCTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGcc                                                                                                                           <  1:920273/147‑1 (MQ=255)

                            cGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGccc                                                                                                                          >  1:735191/1‑150 (MQ=255)

                            cGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGccc                                                                                                                          >  2:2049562/1‑150 (MQ=255)

                               atatTGGATCTGCCACCGATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTTGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGAt                                                                                                                       <  1:2449139/150‑1 (MQ=255)

                                  ttGGATCTGCCACCCATGAGACGCGTTATGGCATGGTCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAAt                                                                                                                    >  1:397796/1‑150 (MQ=255)

                                  ttGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAAt                                                                                                                    >  1:1515294/1‑150 (MQ=255)

                                    ggATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCg                                                                                                                  <  2:2351482/150‑1 (MQ=255)

                                     gATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGAt                                                                                                                       >  2:514484/1‑144 (MQ=255)

                                     gATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGc                                                                                                                 >  1:1115118/1‑150 (MQ=255)

                                        cTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCg                                                                                                              <  1:1963615/150‑1 (MQ=255)

                                         tGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCgg                                                                                                             >  1:318419/1‑150 (MQ=255)

                                         tGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCgg                                                                                                             >  2:2382609/1‑150 (MQ=255)

                                               ccATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCtttt                                                                                                       <  1:1342842/150‑1 (MQ=255)

                                                 aTGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGCTTAGTCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTAc                                                                                                     <  1:1371268/150‑1 (MQ=255)

                                                       cgcgTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTa                                                                                               <  1:785201/150‑1 (MQ=255)

                                                       cgcgTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTa                                                                                               >  1:1304475/1‑150 (MQ=255)

                                                        gcgTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTAt                                                                                              >  2:1899520/1‑150 (MQ=255)

                                                          gTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTg                                                                                            >  1:221828/1‑150 (MQ=255)

                                                           ttATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTgg                                                                                           <  2:318419/150‑1 (MQ=255)

                                                           ttATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTgg                                                                                           <  2:1061097/150‑1 (MQ=255)

                                                               ggcatggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTg                                                                                       >  1:660450/1‑150 (MQ=255)

                                                                gcatggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGc                                                                                      <  2:1747719/150‑1 (MQ=255)

                                                                gcatggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGc                                                                                      <  1:2210987/150‑1 (MQ=255)

                                                                 catggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCa                                                                                     >  2:1486393/1‑150 (MQ=255)

                                                                  atggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAg                                                                                    <  2:1355334/150‑1 (MQ=255)

                                                                  atggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAg                                                                                    <  2:269645/150‑1 (MQ=255)

                                                                  atggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAg                                                                                    <  2:372759/150‑1 (MQ=255)

                                                                     gccgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTaa                                                                                 >  2:981919/1‑150 (MQ=255)

                                                                      ccgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAAc                                                                                <  2:802641/150‑1 (MQ=255)

                                                                      ccgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAAc                                                                                <  1:65986/150‑1 (MQ=255)

                                                                      ccgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAAc                                                                                <  1:1394076/150‑1 (MQ=255)

                                                                       cgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAAct                                                                               >  1:1055959/1‑150 (MQ=255)

                                                                       cgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAAct                                                                               <  2:922886/150‑1 (MQ=255)

                                                                         ccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCAGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAActgc                                                                             >  1:1114708/1‑150 (MQ=255)

                                                                          cTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAActgct                                                                            <  1:2049562/150‑1 (MQ=255)

                                                                          cTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTC‑TTATTGGGTTGCAGTAActgctg                                                                           <  2:759883/150‑1 (MQ=255)

                                                                            gtgCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGt                                                                          <  2:188369/150‑1 (MQ=255)

                                                                              gCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAg                                                                        >  2:892626/1‑150 (MQ=255)

                                                                              gCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAg                                                                        >  1:1508919/1‑150 (MQ=255)

                                                                               ccGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGt                                                                       >  1:1304189/1‑150 (MQ=255)

                                                                                cGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGCGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTc                                                                      <  1:28764/150‑1 (MQ=255)

                                                                                      tAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGcc                                                                >  1:1181013/1‑150 (MQ=255)

                                                                                        aTTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTg                                                              >  1:889993/1‑150 (MQ=255)

                                                                                        aTTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCAGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCAGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTg                                                              >  1:2084354/1‑150 (MQ=255)

                                                                                          ttgattgaTGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGAt                                                            <  2:1853229/150‑1 (MQ=255)

                                                                                          ttgattgaTGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGAt                                                            >  2:297052/1‑150 (MQ=255)

                                                                                           tgattgatGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGAtt                                                           <  1:1329609/150‑1 (MQ=255)

                                                                                            gattgatGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTa                                                          <  2:1579070/150‑1 (MQ=255)

                                                                                             attgatGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTaa                                                         >  1:454936/1‑150 (MQ=255)

                                                                                                    cGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTg                                                  >  2:104951/1‑150 (MQ=255)

                                                                                                     gTTACAAGTAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGGTGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTTCTGTAGTCCAGGCCTGATTAAACGCCTGa                                                 <  2:2084354/150‑1 (MQ=255)

                                                                                                     gTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGa                                                 >  2:1004431/1‑150 (MQ=255)

                                                                                                     gTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGa                                                 >  1:1684448/1‑150 (MQ=255)

                                                                                                       tACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATg                                               >  1:440280/1‑150 (MQ=255)

                                                                                                        aCAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGc                                              <  1:575788/150‑1 (MQ=255)

                                                                                                           tggAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCAGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCtgt                                           >  2:2213935/2‑150 (MQ=255)

                                                                                                           aGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCtgt                                           >  2:2222615/1‑150 (MQ=255)

                                                                                                           aGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCtgt                                           >  2:191687/1‑150 (MQ=255)

                                                                                                            ggAAAGGTTGAGAAGAACTGTGTGAATCAGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCtgtg                                          >  1:2127472/1‑150 (MQ=255)

                                                                                                                 ggTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGt                                     <  1:741565/150‑1 (MQ=255)

                                                                                                                  gTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTa                                    >  1:524136/1‑150 (MQ=255)

                                                                                                                   ttGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTaa                                   >  1:1972203/1‑150 (MQ=255)

                                                                                                                     gagaAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATg                                 >  1:2145460/1‑150 (MQ=255)

                                                                                                                      agaagaACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATgg                                <  1:711026/150‑1 (MQ=255)

                                                                                                                      agaagaACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATgg                                <  2:1914770/150‑1 (MQ=255)

                                                                                                                       gaagaaCTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGgc                               >  1:2169962/1‑150 (MQ=255)

                                                                                                                       gaagaaCTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGgc                               >  2:2453895/1‑150 (MQ=255)

                                                                                                                         agaaCTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGgcgc                             <  2:1228987/150‑1 (MQ=255)

                                                                                                                         agaaCTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGgcgc                             <  2:221828/150‑1 (MQ=255)

                                                                                                                          gaaCTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCa                            >  2:2095729/1‑150 (MQ=255)

                                                                                                                               gtgtgAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAg                       >  1:197134/1‑150 (MQ=255)

                                                                                                                                tgtgAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGt                      >  2:455108/1‑150 (MQ=255)

                                                                                                                                   gAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTg                   >  2:929569/1‑150 (MQ=255)

                                                                                                                                    aaTCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGt                  <  2:1930716/150‑1 (MQ=255)

                                                                                                                                         cgcagTCGCGGCTTAGTCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTatat             >  2:215415/5‑150 (MQ=255)

                                                                                                                                          gcagTCGCGGCTTAGTCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATAtt            >  1:1570916/4‑150 (MQ=255)

                                                                                                                                          gcagTCGCGGCTTAGTCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATAtt            >  1:2277478/4‑150 (MQ=255)

                                                                                                                                           cagTCGCGGCTTAGTCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATAtta           >  2:217063/3‑149 (MQ=255)

                                                                                                                                           cagTCGCGGCTTAGTCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTc           >  2:1212006/3‑150 (MQ=255)

                                                                                                                                              tCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACt        >  1:886002/1‑150 (MQ=255)

                                                                                                                                                gcgGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTg      >  1:1103844/1‑150 (MQ=255)

CTCTCAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTG  >  NC_000913/3718132‑3718429

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