Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,517,171	G→A	L46F (CTT→TTT)	damX ←	cell division protein DamX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,517,171	0	G	A	100.0%	202.5 / NA	67	L46F (CTT→TTT)	damX	cell division protein DamX
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/37);  total (30/37)

CTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAGTTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCTGATTGACAATCG  >  NC_000913/3516931‑3517416                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       cTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATatc                                                                                                                                                                                                                                               >  1:176230/1‑250 (MQ=255)              accgaaaccgatGATCAAAAGAAGTAGAACCAGAATGACGACTCCCATCATTATATACTGACGAGAAGCGGGTTTACTTGCTGCTCTTTTGCACTTACGCGGACGACGCTCTACTCGATCTTCATCCACTGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGGTCCTCATTGCGCTATTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGGCATCAAATTCAATTTCATCAAAATTGATCt                                                                                                                                                                                                                                 <  2:404439/251‑1 (MQ=255)                        atgatCTACAGCAGTAGAACCAGAATGCCGACGCACATCATCATATACTGGCGAAAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTAGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGTACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATCTCATCAAAATTGATCTGCGGTTCGcc                                                                                                                                                                                                                       >  2:587966/1‑251 (MQ=255)                           atcaaatgaagTAGAACCAGAATGCCGACGCCCATCATAATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTAAGGAGGACGACGCTCTACGCGCTCTGAATCCACGGTTTCATCTTCGGATTCGTCAATTTTTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTCCGCACGAGTCGGGCGGCGATCGTAAGTGTCTGCCAATTAAATATCATCAAAAGTGATCTGCGGTTCGccccg                                                                                                                                                                                                                    >  2:646894/10‑248 (MQ=255)                             ccaCAGCAGTAGACCCAGAATGCCGACGCCCTTCTTCATATACTGACGATAATCGGGTTTGCTGGCTTCTTTTTTGCGCCTACGCGGGCGACGCGCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACTGCGATTGTCAGTGTCATCAAATTCAATATCTTCAAAATTGATCTGCGGTTCGCCACGtt                                                                                                                                                                                                                  <  2:202797/249‑1 (MQ=255)                             cAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGtt                                                                                                                                                                                                                  <  1:251193/251‑1 (MQ=255)                               aCAGCAGTAGAACCCGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTc                                                                                                                                                                                                                 >  1:452392/1‑250 (MQ=255)                                                  gCCGACGCCCTACATCATATACTGACGAGAAGCGGGTTTACGTGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCTCGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCag                                                                                                                                                                                            <  2:491582/251‑1 (MQ=255)                                                    cGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCagaa                                                                                                                                                                                          <  1:46528/251‑1 (MQ=255)                                                     gACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGATATCCCCGATCTTCATCCACGGTTTCATATTCGGATTCGTCCATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGGCGGACGGCGATCGTCAGTGTCATCAAATTCAATCTCATCAAAATTGATCTGCGGATCGCCACGTTcctaccgctcaaaacgttacaaa                                                                                                                                                                                           >  2:277767/1‑228 (MQ=255)                                                       cGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACTCTCTACGCGCTCTTTATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGAt                                                                                                                                                                                       <  2:32030/251‑1 (MQ=255)                                                                 catATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAAc                                                                                                                                                                             <  2:526249/251‑1 (MQ=255)                                                                 catATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCGTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAAc                                                                                                                                                                             <  2:664155/251‑1 (MQ=255)                                                                         gACGAGAAGCTGGTTTACTGGCTCCTTTATTGCGCCTACGCGGACGTCGCTCTACGCGCTCTTCATCCTCGGTTTCATCTGCCGATACGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGTTCGTCAGTGTCATCAAAGTCAATATCATCAAAATTGACCGGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAg                                                                                                                                                                       <  1:472433/249‑1 (MQ=255)                                                                         gACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGt                                                                                                                                                                      >  1:473180/1‑250 (MQ=255)                                                                                  cGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCg                                                                                                                                                            >  1:383927/1‑251 (MQ=255)                                                                                   gggTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATAGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑AACGCTCAGAAAGTTCAGAAGATTGACGAGAACGACCAGTacgacgctagc                                                                                                                                                           >  2:479942/1‑246 (MQ=255)                                                                                    ggTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGCGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGCACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCt                                                                                                                                                          >  2:223992/1‑251 (MQ=255)                                                                                       ttactggttgtcttgttGCGCTTACGCTGACGACGCTCCACCATCTCGTCCGCCACGGTTGCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCTGTGCGCTCTTTTTGCGCACTAGTCGGACGGCGATCGGCAGTGTCATCAAATGCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACTTTCAGAATATTGACGAGATCGACCATTACGCCGATCGCTgg                                                                                                                                                        <  2:678861/235‑1 (MQ=255)                                                                                         aCTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACATTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGAt                                                                                                                                                     >  1:637374/1‑251 (MQ=255)                                                                                          cTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGACCGACAAGTACGACGCTCGCTGGGATc                                                                                                                                                    >  2:83679/1‑251 (MQ=255)                                                                                            ggCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCgg                                                                                                                                                  >  1:243597/1‑251 (MQ=255)                                                                                                gctTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGttt                                                                                                                                              >  1:259634/1‑251 (MQ=255)                                                                                                      ttGCGCTTTCGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGTTTCCTTATTGCGCTCTTTTTGCGCACGAGCCGGACGTCGATCGTCATTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTGCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGGACTACGATCGCTGGGATCGTGTTTCAGCt                                                                                                                                         <  2:370532/250‑1 (MQ=255)                                                                                                       tGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTc                                                                                                                                        <  1:291818/250‑1 (MQ=255)                                                                                                                gcgcacgacgCTCGTCGCGCTCTTCATCCACCGGTTCATCCTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCGCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCAGCAAATTCAATATCATCAAAATTGATCTGCGGGTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGGCGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTggg                                                                                                                              <  2:340930/247‑2 (MQ=255)                                                                                                                  ggACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGCTCGCTGGGATCGGGTTTAAGCTCGTCTTCTGGttt                                                                                                                            >  2:435129/1‑251 (MQ=255)                                                                                                                  cggcgctcgtccgatctctctTCATCCACGGTTTCATCTTCGGTTGCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATAACATCAAAATTGATCTGCGGTTCGCCACTTTCAG‑TACGCTCAGAACGTTCAGAATATTTACGAGATCGACCAGTACGTCGATCGGTGGGATCGGGTTTCAGCTCGCCTTCTGGttt                                                                                                                            <  2:628560/234‑1 (MQ=255)                                                                                                                   gacgacGCACTTCGCGCTCTTCAACCACGGCTTCAGCTTCGGATTCGTCAATTTCTTCCTCGATTTCCAGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTTTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGATCGTGCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGCTTg                                                                                                                           <  1:232364/251‑1 (MQ=255)                                                                                                                    acgacgCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGa                                                                                                                          <  1:82574/251‑1 (MQ=255)                                                                                                                      gacgCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGCCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAAt                                                                                                                        >  1:413875/1‑251 (MQ=255)                                                                                                                              aCGCGCTCTTCATCCACGGTTTCCTCTTCGTATTCGTCAATTTCTTCTGCGATTTCCGGTTCCGCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCa                                                                                                                 <  2:491895/250‑1 (MQ=255)                                                                                                                              aCGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCAt                                                                                                                <  1:313523/251‑1 (MQ=255)                                                                                                                               cgcgCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAAGTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTGCTGGTTTGAATTCATCCAtt                                                                                                               <  2:661650/251‑1 (MQ=255)                                                                                                                                   ctctTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCCCGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGttt                                                                                                                            >  1:628560/1‑234 (MQ=255)                                                                                                                                         aTCCTCGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCa                                                                                                     <  1:366819/251‑1 (MQ=255)                                                                                                                                          tCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGATTTCATCCTTTTAACACCCCAc                                                                                                    >  2:306589/1‑251 (MQ=255)                                                                                                                                           ccACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACt                                                                                                   >  1:532209/1‑251 (MQ=255)                                                                                                                                                gTTTCATCTTCTGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTTGTTTGAATTCATCCATTTAACACCCCCCTAAAAg                                                                                              <  1:596778/251‑1 (MQ=255)                                                                                                                                                  ttCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTCAAAGGt                                                                                            >  1:114810/1‑251 (MQ=255)                                                                                                                                                   tCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGt                                                                                            <  1:83679/250‑1 (MQ=255)                                                                                                                                                    cATCTTCGGATTCGGCAATGTCTTCTGCGATTTCCGGTTCCTCATTGCGCTCTGTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCGTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTa                                                                                          <  2:224943/251‑1 (MQ=255)                                                                                                                                                      tCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGATTAAt                                                                                        >  1:580477/1‑251 (MQ=255)                                                                                                                                                         tCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCGG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAAtt                                                                                                                       >  1:672554/1‑217 (MQ=255)                                                                                                                                                         tCGGACTCGTCAAGTTCTTCTTCGATGTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGATTCGGACGGCGATCGTCAGTGTCATCAAAGTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCGG‑GACGCTCAGAACGTGCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCGTCTGGTTTGAAtt                                                                                                                       <  2:672554/217‑1 (MQ=255)                                                                                                                                                          cGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCtt                                                                                    <  2:189949/251‑1 (MQ=255)                                                                                                                                                            gATTCGTCAATTTCGTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCt                                                                                     <  1:197527/248‑1 (MQ=255)                                                                                                                                                            gATTCGCCAATTTATTCTTCCATTTCCGGTTCCTCATAGCGCTCTTTTTGCGAACGAGTAGGACGGCGCTAGTCAGTGTCATCAAATTCACTATCATCAAAATTGCTCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATAGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAAAACCCCGCTAAAAGGTTAATGCttc                                                                                   >  2:322461/1‑249 (MQ=255)                                                                                                                                                                              tcGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTAAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAAc                                                                   <  1:430363/248‑1 (MQ=255)                                                                                                                                                                              tcGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCCGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCt                                                                 <  1:677102/250‑1 (MQ=255)                                                                                                                                                                              tcGATTTCCCGTTCCTCATTTCTCTCTTTTTGCGCACGATCCGGACGGCTATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCTGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACTAGAACGACCAGTACGACGATCGCTTGGATCGTGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTTCCTCGTTGCAATTAACCt                                                                 <  1:677120/250‑1 (MQ=255)                                                                                                                                                                                     ccGGTTCCTCATTGCGCGCTTTTTGCGAACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTGCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGTCCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACGCCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTaa                                                         <  2:195270/251‑1 (MQ=255)                                                                                                                                                                                      cGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCCGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGCATTCATCCATTTACCACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAAt                                                        >  2:475923/1‑251 (MQ=255)                                                                                                                                                                                        gTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCCTTTAACACCCCCCTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAAt                                                        >  2:52326/1‑249 (MQ=255)                                                                                                                                                                                              cATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTg                                                                                                                                                                                     >  1:149064/1‑118 (MQ=255)                                                                                                                                                                                              cATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTg                                                                                                                                                                                     <  2:149064/118‑1 (MQ=255)                                                                                                                                                                                               aTTGTGCTCTTTTTGCGCACGAGTCCGACGGCGATCGTTAGGGTCATCAAAATCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACTATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGcgt                                               <  2:458436/251‑3 (MQ=255)                                                                                                                                                                                                 tGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTg                                              <  1:198277/250‑1 (MQ=255)                                                                                                                                                                                                    gcTCTTTTTGCGCAAGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGAAATTAACCTGAAGCTAATAGACCGCTTGATa                                           >  1:510757/1‑250 (MQ=255)                                                                                                                                                                                                     ctctTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATa                                           >  1:525516/1‑249 (MQ=255)                                                                                                                                                                                                           ttGCGAACGAGTCGGACGGCGATCGTCAGGGTCATCAAATTCAATATCATCAAAAGTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGcc                                    <  2:243597/250‑1 (MQ=255)                                                                                                                                                                                                                   gAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTc                           >  1:241158/1‑251 (MQ=255)                                                                                                                                                                                                                        ggACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCAACTAAAAGGTTAATGCTTACAACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCttttt                      >  2:9264/1‑248 (MQ=255)                                                                                                                                                                                                                                     cAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCTGAtt          <  1:280543/250‑1 (MQ=255)                                                                                                                                                                                                                                         gtCAACAAATTCAATATCATCAAAATTTATCTGCGGGTCGCCACTTTTAGGTA‑GCTCAGAACGTTCAGAATTCTGACGAGAACGAACATTGCTACTATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAAGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCTGATTGACaa     <  2:383927/251‑1 (MQ=255)                                                                                                                                                                                                                                          tcatcaAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCTGATTGACAAt    >  1:70624/1‑251 (MQ=255)                                                                                                                                                                                                                                            atcaAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCTGATTGACaagcg  <  2:413875/251‑4 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       CTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAGTTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAG‑TACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCTGATTGACAATCG  >  NC_000913/3516931‑3517416

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