Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,500,389	(T)8→7	coding (61/981 nt)	ydcK ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,500,382	0	T	.	100.0%	368.0 / NA	77	coding (68/981 nt)	ydcK	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (42/35);  total (42/35)

CCGCTCCGCCGAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACATTTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGAAAATTCCATAAAAT  >  NC_000913/1500233‑1500528                                                                                                                                                      |                                                                                                                                                   ccGCTCCGCCGAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAtt                                                                                                                                                   <  1:16558/151‑1 (MQ=255)  cGCTCCGCCGAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAttt                                                                                                                                                  <  1:684553/151‑1 (MQ=255)   gCTCCGCCGAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAtttt                                                                                                                                                 <  2:1157396/151‑1 (MQ=255)    cTCCGCCGAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAtttt                                                                                                                                                 <  1:85049/150‑1 (MQ=255)     tCCGCCGAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAtttttt                                                                                                                                               >  2:988127/1‑151 (MQ=255)      ccgccgAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAttttttt                                                                                                                                              <  1:1267154/151‑1 (MQ=255)      ccgccgAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAttttttt                                                                                                                                              >  1:858532/1‑151 (MQ=255)      ccgccgAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAttttttt                                                                                                                                              <  2:516389/151‑1 (MQ=255)      ccgccgAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAttttt                                                                                                                                                <  2:448498/149‑1 (MQ=255)       cgccgAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAttttttt                                                                                                                                              >  1:277271/1‑150 (MQ=255)       cgccgAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAttttttt                                                                                                                                              >  1:93547/1‑150 (MQ=255)        gccgAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAtttttttag                                                                                                                                            <  2:1347090/151‑3 (MQ=255)          cgAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAtttttttagtg                                                                                                                                          <  2:1183282/151‑5 (MQ=255)          cgAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAtttttttagtg                                                                                                                                          >  2:109320/1‑147 (MQ=255)          cgAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAtttttttagtg                                                                                                                                          <  1:1144600/151‑5 (MQ=255)           gAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGc                                                                                                                                        >  2:965040/1‑151 (MQ=255)            aaCGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGcc                                                                                                                                       >  1:942930/1‑151 (MQ=255)              cGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCAt                                                                                                                                     <  2:1288911/151‑1 (MQ=255)                cAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATc                                                                                                                                    <  2:355730/150‑1 (MQ=255)                cAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCt                                                                                                                                   >  1:24633/1‑151 (MQ=255)                cAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCt                                                                                                                                   <  2:794303/151‑1 (MQ=255)                  aTAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTc                                                                                                                                 >  2:969927/1‑151 (MQ=255)                  aTAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTc                                                                                                                                 >  1:621130/1‑151 (MQ=255)                   tAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTcc                                                                                                                                >  2:1065784/1‑151 (MQ=255)                    aGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCt                                                                                                                               <  1:583319/151‑1 (MQ=255)                           tGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATa                                                                                                                           >  1:755747/1‑148 (MQ=255)                            ggTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACt                                                                                                                        <  1:1386653/150‑1 (MQ=255)                            ggTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTa                                                                                                                       <  1:916207/151‑1 (MQ=255)                                gTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAgg                                                                                                                   >  2:113703/1‑151 (MQ=255)                                 tAAACCCAGCAATTTCCTTCTTGTGCT‑ACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAg                                                                                                                    <  1:620876/148‑1 (MQ=255)                                        aGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGt                                                                                                           >  2:9019/1‑151 (MQ=255)                                             tttCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGttcttc                                                                                                      >  2:1332734/1‑151 (MQ=255)                                              ttCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGttcttc                                                                                                      >  2:109158/1‑150 (MQ=255)                                              ttCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGttcttc                                                                                                      <  1:691636/150‑1 (MQ=255)                                               tCCTTCTTGTTCTAACACCTTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCa                                                                                                     <  1:1387566/150‑1 (MQ=255)                                                  tacctgtgCTAACACCGTTTCGCGTTCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTCCGGTCTTAACACA‑TTTTCTTAGTGCCATCTTCCTGCTAACTAAAGGCTCGCTGCTCTTCACt                                                                                                   >  1:1288304/5‑149 (MQ=255)                                                   tcttGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACtt                                                                                                  <  2:435528/149‑1 (MQ=255)                                                           tAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGTTCGCTGTTCTTCACTTAGACGat                                                                                           >  1:866724/1‑148 (MQ=255)                                                           tAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATAtt                                                                                        >  1:1290441/1‑151 (MQ=255)                                                           tAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATAtt                                                                                        >  1:224288/1‑151 (MQ=255)                                                            aaCACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATAttt                                                                                       <  1:906170/151‑1 (MQ=255)                                                             accccGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCCCATCGTTAAAATCGCTCATGGCGATAATCGGCCGTAATACCACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAAGCTCGCTGTTCTTCACTTATACGatat                                                                                         <  1:1084144/145‑1 (MQ=255)                                                              cacCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTgc                                                                                     <  2:737329/151‑1 (MQ=255)                                                               acCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTgcg                                                                                    >  1:923442/1‑151 (MQ=255)                                                               acCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTgcg                                                                                    >  2:730087/1‑151 (MQ=255)                                                                  gTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCAt                                                                                 >  2:727540/1‑151 (MQ=255)                                                                    ttCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGTCCATATTTGCGCAt                                                                                 >  1:551345/1‑149 (MQ=255)                                                                    ttCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATa                                                                                <  1:434672/150‑1 (MQ=255)                                                                     tCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAg                                                                               <  1:755406/150‑1 (MQ=255)                                                                        cgatcgatCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGt                                                                            <  2:172936/150‑1 (MQ=255)                                                                         gatcgatcCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGttt                                                                          >  2:434449/1‑151 (MQ=255)                                                                         gatcgatcCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGttt                                                                          <  2:513444/151‑1 (MQ=255)                                                                            cgatcCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTcc                                                                       >  1:262485/1‑151 (MQ=255)                                                                            cgatcCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTcc                                                                       >  2:638152/1‑151 (MQ=255)                                                                                ccAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCcaca                                                                   <  1:849659/151‑1 (MQ=255)                                                                                 cAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCcac                                                                    >  2:348943/1‑149 (MQ=255)                                                                                        cAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACtt                                                           <  1:738924/151‑1 (MQ=255)                                                                                        cAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACtt                                                           >  1:260523/1‑151 (MQ=255)                                                                                           cGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACt                                                        <  2:438554/151‑1 (MQ=255)                                                                                           cGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACt                                                        <  1:983300/151‑1 (MQ=255)                                                                                             gTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACtat                                                      >  2:1038136/1‑151 (MQ=255)                                                                                              ttCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACtat                                                      >  1:1183120/1‑150 (MQ=255)                                                                                                ccTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATaa                                                    <  2:168142/150‑1 (MQ=255)                                                                                                ccTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAAc                                                   >  2:1127661/1‑151 (MQ=255)                                                                                                  tGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCg                                                 >  1:866719/1‑151 (MQ=255)                                                                                                  tGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCg                                                 >  2:604926/1‑151 (MQ=255)                                                                                                      aTCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAg                                              <  2:641275/150‑1 (MQ=255)                                                                                                       tCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGc                                             >  1:272587/1‑150 (MQ=255)                                                                                                              gTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTg                                     <  2:1309012/151‑1 (MQ=255)                                                                                                                tAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGc                                    <  2:83456/150‑1 (MQ=255)                                                                                                                 aaaaTCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCg                                   <  2:1305471/150‑1 (MQ=255)                                                                                                                    aTCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATCTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGctc                               >  1:1358810/1‑151 (MQ=255)                                                                                                                      cGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGctct                              <  2:613942/150‑1 (MQ=255)                                                                                                                      cGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGctct                              <  2:415185/150‑1 (MQ=255)                                                                                                                      cGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTg                             >  1:1106236/1‑151 (MQ=255)                                                                                                                          aaTGGCGATAATCTGCCGTAATAACACA‑TTTTTTAAGACCCATCTTCCGGAAAACAAAAGGACCGCTGTTCTTCACTTAAACGATTTTTGGGCATAGCGTTTTCCCACAGGCACCTTACTATAACCGGACAAAATCCTGCGGCCCTGGc                           >  2:1025791/1‑149 (MQ=255)                                                                                                                           aTGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTa                         >  1:1399077/1‑150 (MQ=255)                                                                                                                              gCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCt                     <  2:200061/151‑1 (MQ=255)                                                                                                                                 aTAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCg                   >  2:196167/1‑150 (MQ=255)                                                                                                                                      cTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGaaa              >  2:182370/1‑150 (MQ=255)                                                                                                                                        gCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGaaaa             <  2:296761/149‑1 (MQ=255)                                                                                                                                        gCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGaaaa             <  2:689305/149‑1 (MQ=255)                                                                                                                                         ccGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGAAAATTc          >  2:536652/1‑151 (MQ=255)                                                                                                                                          cGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGAAAATTc          >  2:142305/1‑150 (MQ=255)                                                                                                                                            taataaCACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGAAAATTCCAt       >  1:252078/1‑151 (MQ=255)                                                                                                                                            taataaCACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGAAAATTCCAt       >  2:122097/1‑151 (MQ=255)                                                                                                                                             aataaCACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGAAAATTCCAt       >  2:170002/1‑150 (MQ=255)                                                                                                                                             aataaCACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGAAAATTCCAt       <  1:1399975/150‑1 (MQ=255)                                                                                                                                             aataaCACA‑TTTTTTTAGTGCCATCTGCTTGATAACTAAAGGAACGCTGTTCTTCACTTAGACTATATTTGCGCATAGCGTTTTCCCACAGGCGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGATATGCTCGAGAAAATTCCa        >  1:12976/1‑149 (MQ=255)                                                                                                                                              ataaCACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGAAAATTCCATaa     <  2:621161/151‑1 (MQ=255)                                                                                                                                                aaCACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGAAAATTCCATaaa    <  1:1287599/150‑1 (MQ=255)                                                                                                                                                 acacA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGAAAATTCCATAAAAt  <  1:907621/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CCGCTCCGCCGAACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACATTTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGAAAATTCCATAAAAT  >  NC_000913/1500233‑1500528

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