Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,276,133	A→G	K221K (AAA→AAG)	yjcC →	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,276,133	0	A	G	100.0%	323.4 / NA	91	K221K (AAA→AAG)	yjcC	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (45/46);  total (45/46)

TGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCT  >  NC_000913/4275985‑4276280                                                                                                                                                     |                                                                                                                                                    tGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCg                                                                                                                                                   >  2:805012/1‑151 (MQ=255)    aaTGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGccc                                                                                                                                                >  1:850134/1‑151 (MQ=255)     aTGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCa                                                                                                                                               <  2:850148/151‑1 (MQ=255)      tGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCa                                                                                                                                               >  2:670865/1‑150 (MQ=255)       ggggTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAAt                                                                                                                                             <  2:421157/151‑1 (MQ=255)          gtgtgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGc                                                                                                                                          >  2:201352/1‑151 (MQ=255)           tgtgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGc                                                                                                                                          >  1:125526/1‑150 (MQ=255)            gtgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAg                                                                                                                                        <  2:305552/151‑1 (MQ=255)                aTGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCAt                                                                                                                                    >  2:71523/1‑151 (MQ=255)                  gATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCAttg                                                                                                                                  >  1:463708/1‑151 (MQ=255)                  gATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCAttg                                                                                                                                  >  2:582280/1‑151 (MQ=255)                  gATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCAttg                                                                                                                                  >  1:207657/1‑151 (MQ=255)                      cGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGc                                                                                                                              >  2:58125/1‑151 (MQ=255)                      cGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGc                                                                                                                              <  2:1087112/151‑1 (MQ=255)                          gACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACt                                                                                                                          <  2:709505/151‑1 (MQ=255)                            cGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTc                                                                                                                        >  2:658999/1‑151 (MQ=255)                                     ttttCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAAc                                                                                                                <  2:292771/150‑1 (MQ=255)                                      tttCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACg                                                                                                               <  2:481172/150‑1 (MQ=255)                                           cGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCtt                                                                                                           <  1:917881/149‑1 (MQ=255)                                           cGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTtat                                                                                                         <  1:832192/151‑1 (MQ=255)                                            gTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTtat                                                                                                         >  1:177653/1‑150 (MQ=255)                                              tAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATaa                                                                                                       <  1:1110130/150‑1 (MQ=255)                                                 gCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAAccc                                                                                                    >  1:346240/1‑150 (MQ=255)                                                   aaaGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCAtt                                                                                                 <  2:342384/151‑1 (MQ=255)                                                   aaaGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCAtt                                                                                                 <  2:1163721/151‑1 (MQ=255)                                                        ggCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTtat                                                                                            <  1:281084/151‑1 (MQ=255)                                                         gCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTtat                                                                                            >  2:891078/1‑150 (MQ=255)                                                             cGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAatc                                                                                        <  2:249032/150‑1 (MQ=255)                                                             cGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAatc                                                                                        >  1:741484/1‑150 (MQ=255)                                                             cGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAat                                                                                         <  1:1146144/149‑1 (MQ=255)                                                              ggCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAatcat                                                                                      <  1:1128261/151‑1 (MQ=255)                                                              ggCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAatcat                                                                                      >  2:25344/1‑151 (MQ=255)                                                               gcagcaACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAatcatc                                                                                     <  2:271033/151‑1 (MQ=255)                                                                cagcaACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCt                                                                                    <  1:638117/151‑1 (MQ=255)                                                                    aaCGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTAtt                                                                                >  2:678943/1‑151 (MQ=255)                                                                      cGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTAtttt                                                                              >  2:578308/1‑151 (MQ=255)                                                                       gTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTAttttt                                                                             >  1:467044/1‑151 (MQ=255)                                                                         tttCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTg                                                                            <  1:95713/150‑1 (MQ=255)                                                                              ccGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGtt                                                                        <  1:784995/149‑1 (MQ=255)                                                                              ccGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGt                                                                         <  2:3637/148‑1 (MQ=255)                                                                                   gATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCccg                                                                   <  2:177655/149‑1 (MQ=255)                                                                                     ttCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCccgccg                                                                <  2:394717/150‑1 (MQ=255)                                                                                     ttCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCccgccg                                                                <  1:706409/150‑1 (MQ=255)                                                                                     ttCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCccgccg                                                                >  1:394714/1‑150 (MQ=255)                                                                                         tttGAAGGATTTACCCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTAt                                                            <  2:306317/150‑1 (MQ=255)                                                                                         tttGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTAtt                                                           >  2:748638/1‑151 (MQ=255)                                                                                         tttGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTAtt                                                           >  1:540254/1‑151 (MQ=255)                                                                                         tttAAAGGTTTTAACCGTAAAAAAAAATGGTTTTTTATTTGCGACAGTTTTTTCGAAAAAGCCCCCATTTGCAGCCTTTGTTGCTATTTCATATCACCGTCTTATAACCCATTTTTATAATCATTTTATTTTTGCTTGGCCCGCCGGt                                                              <  1:606492/148‑1 (MQ=255)                                                                                          ttGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTAttt                                                          >  1:424133/1‑151 (MQ=255)                                                                                                  tttAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGt                                                  >  1:639093/1‑151 (MQ=255)                                                                                                  tttAACCGTAAAAAAAAATGGTTATTTTTATGCAACGGTTTATTCGACAAAGCGCCCAATGGCAGCCATTGTTTCGACTTCATATAACCGTCTTATAACCCATTTTTATAATCATTTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAg                                                   <  2:360407/150‑1 (MQ=255)                                                                                                   ttAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGt                                                  <  2:838235/150‑1 (MQ=255)                                                                                                   ttAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGt                                                  <  1:540169/150‑1 (MQ=255)                                                                                                   ttAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTc                                                 >  1:874497/1‑151 (MQ=255)                                                                                                      aCCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTg                                              <  2:169701/151‑1 (MQ=255)                                                                                                      aCCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTAATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCt                                                <  1:448461/149‑1 (MQ=255)                                                                                                       ccGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTg                                              <  2:547949/150‑1 (MQ=255)                                                                                                       ccGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGt                                             <  2:893662/151‑1 (MQ=255)                                                                                                        cGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCTTTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGtt                                            <  1:805001/151‑1 (MQ=255)                                                                                                         gTACAAATAAATGGCTATTTATATGCGACATTTTTTTCGACAAAGCTCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCTTTGCCCGCCGGTATTTTGGGGAGTCTTGtt                                            <  2:407914/150‑1 (MQ=255)                                                                                                              aaGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCt                                       >  1:354077/1‑150 (MQ=255)                                                                                                              aaGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCt                                       >  1:337188/1‑150 (MQ=255)                                                                                                                 aaatgattgtttATATGGGACAGTTTATTCGACATGGCGCCCAGTTGCAGCCTTTGTTCCGAATTCATATCAACGTTTTATAACCCTTTTTTATAATCATTTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTAt                                     <  1:283372/140‑1 (MQ=255)                                                                                                                 aaaTGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTAtt                                    >  2:22553/1‑150 (MQ=255)                                                                                                                 aaaTGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTAtt                                    >  2:669420/1‑150 (MQ=255)                                                                                                                 aaaTGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTAtt                                    >  2:598647/1‑150 (MQ=255)                                                                                                                 aaaTGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTa                                   <  2:1011241/151‑1 (MQ=255)                                                                                                                    tGGCTATTTATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTAct                                 <  1:245329/150‑1 (MQ=255)                                                                                                                        tATTTATATGCGACATTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCTTCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTg                              <  2:654887/149‑1 (MQ=255)                                                                                                                          tttATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCt                           >  1:1100411/1‑150 (MQ=255)                                                                                                                           ttATATGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTAc                         >  1:418710/1‑151 (MQ=255)                                                                                                                            tataTGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTAc                         <  2:1042289/150‑1 (MQ=255)                                                                                                                            tataTGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACg                        >  2:253070/1‑151 (MQ=255)                                                                                                                             atatGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGt                       >  1:662480/1‑151 (MQ=255)                                                                                                                             atatGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGt                       >  1:531532/1‑151 (MQ=255)                                                                                                                              tatGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGt                       >  1:1145967/1‑150 (MQ=255)                                                                                                                                tGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTa                      >  1:1045902/1‑149 (MQ=255)                                                                                                                                tGCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTAtt                    >  1:16607/1‑151 (MQ=255)                                                                                                                                 gCGACAGTTTATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTc                   >  2:771885/1‑151 (MQ=255)                                                                                                                                        tttATTCGACAAAGCGCCCAATTGCAGCCATTTTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACaaa             <  1:131727/150‑1 (MQ=255)                                                                                                                                        tttATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACaaa             >  2:600030/1‑150 (MQ=255)                                                                                                                                        tttATTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACaa              <  1:31541/149‑1 (MQ=255)                                                                                                                                           aTTCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACt          >  2:452632/1‑150 (MQ=255)                                                                                                                                             tCGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTAtt       >  1:438285/1‑151 (MQ=255)                                                                                                                                              cGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTAtt       <  2:460250/150‑1 (MQ=255)                                                                                                                                              cGACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTAtt       <  2:278984/150‑1 (MQ=255)                                                                                                                                              cGACAAAGCCCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATACCCCATTTTTATAATCATCTTATTTTTGCGTTCCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTAttt      <  2:1059058/151‑1 (MQ=255)                                                                                                                                               gACAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTAttt      >  1:932855/1‑150 (MQ=255)                                                                                                                                                aCAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTAAAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTAt        <  1:672166/147‑1 (MQ=255)                                                                                                                                                 cAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTAtc   >  1:842191/1‑151 (MQ=255)                                                                                                                                                 cAAAGCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTAtc   >  1:1125994/1‑151 (MQ=255)                                                                                                                                                  aaagCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTAtct  >  1:28935/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCT  >  NC_000913/4275985‑4276280

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