BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	473,057	A→G	5.3%	D31G (GAC→GGC)	amtB →	ammonium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	473,057	0	A	G	5.3%	378.3 / 7.9	114	D31G (GAC→GGC)	amtB	ammonium transporter
Reads supporting (aligned to +/- strand): ref base A (56/52); new base G (3/3); total (59/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 115 total aligned reads displayed.

CCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTC  >  NC_000913/472908‑473205

ccGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGa                                                                                                                                                      <  1:1716729/150‑1 (MQ=255)

 cGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATTAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGAc                                                                                                                                                     >  1:79757/1‑150 (MQ=255)

 cGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGAc                                                                                                                                                     >  1:79765/1‑150 (MQ=255)

       aGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCg                                                                                                                                               <  2:2098554/150‑1 (MQ=255)

       aGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCg                                                                                                                                               <  2:2098546/150‑1 (MQ=255)

       aGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCg                                                                                                                                               <  1:2290355/150‑1 (MQ=255)

         ccGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGtt                                                                                                                                             >  1:1782281/1‑150 (MQ=255)

           gacgaAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTa                                                                                                                                           >  2:1080745/1‑150 (MQ=255)

           gacgaAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTa                                                                                                                                           >  2:2319619/1‑150 (MQ=255)

           gacgaAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTa                                                                                                                                           >  2:847753/1‑150 (MQ=255)

            acgaAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTat                                                                                                                                          >  1:1305878/1‑150 (MQ=255)

              gaAGCGGCGCTGTAATCTCTGGCACACAGCAACTGGAACGAAAAATGAAGATAGCGACGATAATATCTGGGCTTGCTTCTCTGGCGATGCTTCCGGGAGTGGTAATGGCTGCTCCTGCGGTGGCCGATAAAGCCGACAATGCGTTTatgt                                                                                                                                        >  2:820824/1‑149 (MQ=255)

               aaGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGGCCATGCGGTTatgat                                                                                                                                       >  1:2307992/1‑150 (MQ=255)

               aaGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTatgat                                                                                                                                       >  2:776356/1‑150 (MQ=255)

               aaGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTatgat                                                                                                                                       >  2:421987/1‑150 (MQ=255)

               aaGCGGCGCTGTAATCTATGGCACACAGCAACAGGAAAGAAAAATGAAGATAGCGAAGATAAAAAAAGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAAAGGCTGCACCTGCGGTGGCCGATAAAGCCGGCAATGCGTTTatgat                                                                                                                                       >  1:619807/1‑150 (MQ=255)

                      gcTGTAATCTCTGGTACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGt                                                                                                                                <  1:918921/150‑1 (MQ=255)

                       cTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTa                                                                                                                               <  1:1394940/150‑1 (MQ=255)

                       cTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTa                                                                                                                               <  2:1664758/150‑1 (MQ=255)

                       cTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTa                                                                                                                               <  1:1394764/150‑1 (MQ=255)

                               tctGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTg                                                                                                                       >  1:313011/1‑150 (MQ=255)

                               tcaggCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGATTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTg                                                                                                                       >  1:1810397/4‑150 (MQ=255)

                                ctGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTgg                                                                                                                      >  1:1393638/1‑150 (MQ=255)

                                 tGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCTATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGAGGCCGATAAAGCCGACACTGAGTTTATGATGATTTGTACTGCGCTGGt                                                                                                                     >  2:1596859/1‑150 (MQ=255)

                                 tGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGt                                                                                                                     >  1:1046418/1‑150 (MQ=255)

                                 tGGCACACAGCAACAGGAAAGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGt                                                                                                                     <  2:1733178/150‑1 (MQ=255)

                                  ggCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTg                                                                                                                    <  1:634838/150‑1 (MQ=255)

                                     acacaGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTg                                                                                                                 >  1:2403274/1‑150 (MQ=255)

                                     acacaGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTg                                                                                                                 >  2:477737/1‑150 (MQ=255)

                                        caGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGttt                                                                                                              <  1:1356964/150‑1 (MQ=255)

                                          gCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTATTGCGCTGGTGCTGTTTAt                                                                                                            <  1:1551064/150‑1 (MQ=255)

                                          gCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTAt                                                                                                            <  1:2389003/150‑1 (MQ=255)

                                            aaCAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTCTGa                                                                                                          >  2:1229296/1‑150 (MQ=255)

                                            aaCAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGa                                                                                                          >  2:22142/1‑150 (MQ=255)

                                                     gAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCgg                                                                                                 <  2:1974274/150‑1 (MQ=255)

                                                      aaaaaTGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCggg                                                                                                <  1:695639/150‑1 (MQ=255)

                                                           tGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTg                                                                                           <  1:1603623/150‑1 (MQ=255)

                                                           tGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTg                                                                                           <  1:1603618/150‑1 (MQ=255)

                                                           tGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTg                                                                                           <  1:1183875/150‑1 (MQ=255)

                                                           tGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTg                                                                                           <  1:1115035/150‑1 (MQ=255)

                                                            gAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGc                                                                                          >  1:404903/1‑150 (MQ=255)

                                                             aaGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGcc                                                                                         <  1:451413/150‑1 (MQ=255)

                                                              aGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGccc                                                                                        <  1:2372336/150‑1 (MQ=255)

                                                              aGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGccc                                                                                        <  1:2372327/150‑1 (MQ=255)

                                                              aGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGccc                                                                                        <  1:2372317/150‑1 (MQ=255)

                                                                aTAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTg                                                                                      >  2:1132807/1‑150 (MQ=255)

                                                                aTAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTg                                                                                      >  2:1132648/1‑150 (MQ=255)

                                                                   gCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGttt                                                                                   >  1:2222568/1‑150 (MQ=255)

                                                                     gacgaTAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTa                                                                                 <  2:372678/150‑1 (MQ=255)

                                                                     gacgaTAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTa                                                                                 <  2:373130/150‑1 (MQ=255)

                                                                        gaTAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACgg                                                                              >  1:1011415/1‑150 (MQ=255)

                                                                         aTAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGt                                                                             >  2:1855010/1‑150 (MQ=255)

                                                                         aTAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGt                                                                             >  2:1550223/1‑150 (MQ=255)

                                                                         aTAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGt                                                                             >  2:1550216/1‑150 (MQ=255)

                                                                           aaaaaCTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTgg                                                                           >  2:1727406/1‑150 (MQ=255)

                                                                           aaaaaCTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTgg                                                                           <  2:1606073/150‑1 (MQ=255)

                                                                           aaaaaCTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTgg                                                                           <  1:1759887/150‑1 (MQ=255)

                                                                                  gggCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGAtt                                                                    <  1:1840091/150‑1 (MQ=255)

                                                                                      ttgcttCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTcgcg                                                                <  1:477737/150‑1 (MQ=255)

                                                                                       tgcttCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCgg                                                               <  2:1657986/150‑1 (MQ=255)

                                                                                       tgcttCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCgg                                                               <  2:1566114/150‑1 (MQ=255)

                                                                                       tgcttCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCgg                                                               <  1:1776017/150‑1 (MQ=255)

                                                                                       tgcttCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCgg                                                               <  1:1775942/150‑1 (MQ=255)

                                                                                        gcttCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGc                                                              <  1:2208518/150‑1 (MQ=255)

                                                                                        gcttCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGc                                                              <  2:862124/150‑1 (MQ=255)

                                                                                              cTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGGCAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAAc                                                        <  2:619807/150‑1 (MQ=255)

                                                                                              cTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGGCAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAAc                                                        <  2:2307992/150‑1 (MQ=255)

                                                                                              cTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGGCAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAAc                                                        <  2:2346856/150‑1 (MQ=255)

                                                                                              cTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAAc                                                        <  1:421987/150‑1 (MQ=255)

                                                                                              cTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAAc                                                        <  1:776356/150‑1 (MQ=255)

                                                                                               tGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACg                                                       >  2:2155464/1‑150 (MQ=255)

                                                                                               tGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACg                                                       >  1:1256040/1‑150 (MQ=255)

                                                                                                   gATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCt                                                   >  1:1579383/1‑150 (MQ=255)

                                                                                                   gATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCt                                                   >  1:1579329/1‑150 (MQ=255)

                                                                                                     tGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGt                                                 >  1:816349/1‑150 (MQ=255)

                                                                                                     tGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGt                                                 >  2:843794/1‑150 (MQ=255)

                                                                                                      gCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTc                                                <  1:2113422/150‑1 (MQ=255)

                                                                                                        ttCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGa                                              <  1:121821/150‑1 (MQ=255)

                                                                                                        ttCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGa                                              >  2:1520167/1‑150 (MQ=255)

                                                                                                        ttCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGa                                              >  2:372945/1‑150 (MQ=255)

                                                                                                         tCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGAt                                             >  1:654287/1‑150 (MQ=255)

                                                                                                                 tGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGc                                     >  1:807229/1‑150 (MQ=255)

                                                                                                                   gTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAg                                   >  2:108175/1‑150 (MQ=255)

                                                                                                                      aTGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAggtg                                >  2:873605/1‑150 (MQ=255)

                                                                                                                           tGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAggtgacggt                           <  1:481106/150‑1 (MQ=255)

                                                                                                                           tGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAggtgacggt                           <  1:1408591/150‑1 (MQ=255)

                                                                                                                             cACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAggtgacggtga                         <  1:2388578/150‑1 (MQ=255)

                                                                                                                                   cGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTg                   <  2:2337585/150‑1 (MQ=255)

                                                                                                                                          cGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGt            >  2:1218751/1‑150 (MQ=255)

                                                                                                                                          cGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGt            >  1:1977880/1‑150 (MQ=255)

                                                                                                                                          cGATAAAGCCGACAATGCGTGTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGAAGGTGAAATTTGCACTGGt            >  2:2032421/1‑150 (MQ=255)

                                                                                                                                           gATAAAGCCGACAATGCGTTTATGATGA‑TTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCt          <  1:1606746/150‑1 (MQ=255)

                                                                                                                                            aTAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCt          <  1:104623/150‑1 (MQ=255)

                                                                                                                                             tAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTg         >  1:96913/1‑150 (MQ=255)

                                                                                                                                             tAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTg         >  1:1072123/1‑150 (MQ=255)

                                                                                                                                              aaaGCCGACTATGCGTTTATGATGATTTGTACTTCGCTTGTGCTGTTTATGACTTTTCCGGTGATTGCCCTGTTTTTCGGTTTGTTGATTCGCGGCAAAAAAGTGCTGTCGATGCTGACGCAGTTGACGGTGACATTTGCAATGGTATGt        <  1:670278/150‑1 (MQ=255)

                                                                                                                                                aGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTAt      >  2:1320329/1‑150 (MQ=255)

                                                                                                                                                  ccGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATtc    >  1:28271/1‑150 (MQ=255)

                                                                                                                                                    gACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATtctc  >  1:587205/1‑150 (MQ=255)

                                                                                                                                                    gACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATtctc  >  1:1733961/1‑150 (MQ=255)

CCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTC  >  NC_000913/472908‑473205

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