BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,617	A→G	D307G (GAC→GGC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,617	0	A	G	100.0%	380.0 / NA	100	D307G (GAC→GGC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (44/56); total (44/56)

Only 100 of 104 total aligned reads displayed.

CGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACC  >  NC_000913/1756471‑1756762

cGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGgcgc                                                                                                                                                >  1:1109218/1‑150 (MQ=255)

cGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGgcgc                                                                                                                                                >  4:236839/1‑150 (MQ=255)

cGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGgcgc                                                                                                                                                >  2:1057121/1‑150 (MQ=255)

cGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCa                                                                                                                                               <  1:1036448/151‑1 (MQ=255)

 gCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGgcgc                                                                                                                                                >  1:1101029/1‑149 (MQ=255)

 gCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAg                                                                                                                                              <  3:300931/151‑1 (MQ=255)

  cccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAg                                                                                                                                              >  1:750932/1‑150 (MQ=255)

   ccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAGAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTg                                                                                                                                            <  3:219430/151‑1 (MQ=255)

   ccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTg                                                                                                                                            <  1:853545/151‑1 (MQ=255)

   ccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTg                                                                                                                                            >  4:158190/1‑151 (MQ=255)

    cAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGa                                                                                                                                           <  1:993174/151‑1 (MQ=255)

        agaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCt                                                                                                                                       >  1:887228/1‑151 (MQ=255)

             atgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCt                                                                                                                                   <  1:1084670/150‑1 (MQ=255)

             atgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTg                                                                                                                                  >  1:1100228/1‑151 (MQ=255)

             atgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTg                                                                                                                                  <  2:754202/151‑1 (MQ=255)

               gatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGt                                                                                                                                >  3:235659/1‑151 (MQ=255)

                atCGAAAAATGTATCCGTGCACGTAGAGTCGTTATCACTGCGACCCAGATGCTGGATTCCAGGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTAGTGGCGCAGTGATGCTGGCTGGt                                                                                                                                <  2:1012089/150‑1 (MQ=255)

                atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTg                                                                                                                               >  2:24915/1‑151 (MQ=255)

                 tCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTg                                                                                                                               >  4:233997/1‑150 (MQ=255)

                 tCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTg                                                                                                                               >  2:589132/1‑150 (MQ=255)

                  cGAAAAATGTATCCGTGCGCGTAAAATCGTTATCACTGCGACCCAGATGGTGGATTCCATGATCAAAACCCCACGCCCGACGCTCGCAGAAGCCGGTGACGTTGCAATCGCAATCCTCCCCGGTCCTGGCGCAGTGATGCTGTCTGGTGa                                                                                                                              <  2:117202/150‑1 (MQ=255)

                  cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGa                                                                                                                              >  4:271937/1‑150 (MQ=255)

                   gAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAAt                                                                                                                            >  3:75119/1‑151 (MQ=255)

                     aaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATcc                                                                                                                          <  2:283506/151‑1 (MQ=255)

                     aaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATcc                                                                                                                          >  2:1045971/1‑151 (MQ=255)

                     aaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATcc                                                                                                                          <  1:462223/151‑1 (MQ=255)

                          gTATCCGTGCGCGGAAAGTCGTTATCACGGCGACCCAGATGCTGGATTCCATGAGCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCAGCCTCGACGGTACTGGCGCAGTGATGCTGTCGGGTGAATCCGCaa                                                                                                                      <  2:1264979/150‑1 (MQ=255)

                          gTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCaa                                                                                                                      <  2:505599/150‑1 (MQ=255)

                            aTCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAg                                                                                                                   <  1:780934/151‑1 (MQ=255)

                              ccGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAgg                                                                                                                  <  2:37427/150‑1 (MQ=255)

                              ccGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGt                                                                                                                 <  2:1124807/151‑1 (MQ=255)

                              ccGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGt                                                                                                                 <  3:335354/151‑1 (MQ=255)

                               cGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTa                                                                                                                <  4:222882/151‑1 (MQ=255)

                                         aaGTCGTTATCACTGCGACCCAGATTCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCt                                                                                                      <  2:772769/151‑1 (MQ=255)

                                         aaGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGGCTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCt                                                                                                      <  1:547219/151‑1 (MQ=255)

                                         aaGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCt                                                                                                      <  3:33788/151‑1 (MQ=255)

                                          aGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTg                                                                                                     >  2:1148300/1‑151 (MQ=255)

                                          aGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTg                                                                                                     <  2:344538/151‑1 (MQ=255)

                                          aGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTg                                                                                                     <  3:427979/151‑1 (MQ=255)

                                           gTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTg                                                                                                     <  2:225299/150‑1 (MQ=255)

                                                tATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGAGGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAcgcg                                                                                               >  1:537699/1‑151 (MQ=255)

                                                 aTCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCgg                                                                                              >  4:201688/1‑151 (MQ=255)

                                                 aTCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCg                                                                                               <  3:150200/150‑1 (MQ=255)

                                                   cACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGtt                                                                                            >  1:207489/1‑151 (MQ=255)

                                                   cACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGtt                                                                                            <  2:1109012/151‑1 (MQ=255)

                                                     cTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTc                                                                                          >  2:404371/1‑151 (MQ=255)

                                                     cTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTc                                                                                          >  1:579094/1‑151 (MQ=255)

                                                      tGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCt                                                                                         >  2:1192258/1‑151 (MQ=255)

                                                       gCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTa                                                                                        <  2:116599/151‑1 (MQ=255)

                                                        cGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTAt                                                                                       >  3:330325/1‑151 (MQ=255)

                                                        cGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTAt                                                                                       <  2:700027/151‑1 (MQ=255)

                                                        cGACCCAGATGCTGGATCCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTAt                                                                                       >  1:1107458/1‑151 (MQ=255)

                                                         gACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATc                                                                                      >  4:379564/1‑151 (MQ=255)

                                                           cccAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCAt                                                                                    <  1:1044202/151‑1 (MQ=255)

                                                           cccAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCAt                                                                                    <  2:752293/151‑1 (MQ=255)

                                                            ccAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATg                                                                                   >  2:497764/1‑151 (MQ=255)

                                                             cAGATGCTGGATTCCATGATCAAAAACCCCCGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATc                                                                                      <  2:1030290/147‑1 (MQ=255)

                                                               gATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCg                                                                                <  1:170044/151‑1 (MQ=255)

                                                                aTGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCg                                                                                <  2:726870/150‑1 (MQ=255)

                                                                  gCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGAc                                                                              <  3:423978/150‑1 (MQ=255)

                                                                    tGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCAt                                                                           >  1:546582/1‑151 (MQ=255)

                                                                       aTTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTg                                                                        >  4:12768/1‑151 (MQ=255)

                                                                        ttCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGc                                                                       >  1:692540/1‑151 (MQ=255)

                                                                         tCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCg                                                                      >  4:361231/1‑151 (MQ=255)

                                                                          ccATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGa                                                                     <  1:441756/151‑1 (MQ=255)

                                                                               aTCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGt                                                                 >  2:1156951/1‑150 (MQ=255)

                                                                               aTCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGt                                                                 >  2:481843/1‑150 (MQ=255)

                                                                                 cAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTacc                                                              <  1:1094988/151‑1 (MQ=255)

                                                                                  aaaaaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTacc                                                              >  2:519883/1‑150 (MQ=255)

                                                                                     aaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTaccgacc                                                          <  1:34521/151‑1 (MQ=255)

                                                                                                gACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGaa                                                >  2:1258556/1‑150 (MQ=255)

                                                                                                 aCTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACa                                              <  2:612956/151‑1 (MQ=255)

                                                                                                  cTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGACGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACa                                              <  2:749183/150‑1 (MQ=255)

                                                                                                     gcgcAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCg                                          >  4:19076/1‑151 (MQ=255)

                                                                                                      cgcAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGt                                         >  1:1186395/1‑151 (MQ=255)

                                                                                                        cAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGtct                                       >  1:659509/1‑151 (MQ=255)

                                                                                                         aGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGtctc                                      <  4:383622/151‑1 (MQ=255)

                                                                                                          gAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCg                                     >  1:586838/1‑151 (MQ=255)

                                                                                                             gCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGt                                  <  3:379559/151‑1 (MQ=255)

                                                                                                              ccGGTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGtt                                 >  4:224164/1‑151 (MQ=255)

                                                                                                                ggTGACGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTca                               <  2:908287/151‑1 (MQ=255)

                                                                                                                     cGTTGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAAt                          <  1:700910/151‑1 (MQ=255)

                                                                                                                        tGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGa                        <  1:736059/150‑1 (MQ=255)

                                                                                                                        tGCAAACGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGAc                       >  1:898853/1‑151 (MQ=255)

                                                                                                                           aaaCGCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAAc                    >  4:165743/1‑151 (MQ=255)

                                                                                                                               gCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGt                 >  2:1236140/1‑150 (MQ=255)

                                                                                                                               gCCATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTa                <  4:338405/151‑1 (MQ=255)

                                                                                                                                ccATCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTaa               <  1:383715/151‑1 (MQ=255)

                                                                                                                                   tCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAAc             <  2:440327/150‑1 (MQ=255)

                                                                                                                                   tCCTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACt            >  1:1083068/1‑151 (MQ=255)

                                                                                                                                    ccTCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTg           <  3:142333/151‑1 (MQ=255)

                                                                                                                                      tCGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTgcg         >  1:778169/1‑151 (MQ=255)

                                                                                                                                       cGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTgcg         <  1:409738/150‑1 (MQ=255)

                                                                                                                                       cGACGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTgcg         <  1:145357/150‑1 (MQ=255)

                                                                                                                                         aCGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCAt      <  1:973488/151‑1 (MQ=255)

                                                                                                                                         aCGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCAt      >  1:543321/1‑151 (MQ=255)

                                                                                                                                         aCGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCAt      >  4:38237/1‑151 (MQ=255)

                                                                                                                                          cGGTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCAt      <  1:62350/150‑1 (MQ=255)

                                                                                                                                            gTACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTAc   <  2:1044680/151‑1 (MQ=255)

                                                                                                                                             tACTGGCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTAcc  >  1:79546/1‑151 (MQ=255)

CGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACC  >  NC_000913/1756471‑1756762

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