Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,867,454	2 bp→AA	coding (97‑98/993 nt)	rpoS ←	RNA polymerase sigma factor RpoS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,867,454	0	T	A	100.0%	231.7 / NA	60	Q33M (CAG→TTG)	rpoS	RNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (24/36);  total (24/36)
*	NC_000913	2,867,455	0	G	A	100.0%	224.4 / NA	60	Q33M (CAG→TTG)	rpoS	RNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/36);  total (24/36)

GACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACT  >  NC_000913/2867212‑2867691                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                              gACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                                                                                                                                                                                                        <  1:16368/251‑1 (MQ=255) gACCACGATTGCCATAACGGCGGGCAACTTTTACCCCCCGACTCAAGTTAATCTCTATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAAGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGGGAATAACCAAACTCACCAAGGGAAAGGTGAGACGCGTCCAACACACCCTGTGTGGCTCCCTGCGAATAA‑AGATCCTCTTCGGCCAAATCGGTATCACTGGGTTCCAATTctact                                                                                                                                                                                                                                        >  1:370986/1‑251 (MQ=255)           gCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACTGTGCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                                                                                                                                                                                                               <  1:143779/250‑1 (MQ=255)           gCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                                                                                                                                                                                                               <  2:492514/250‑1 (MQ=255)                  cggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACCGTTCCTCTTCGGCCAAATCGTTATCCCTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                                                                                                                                                                                                        >  2:443819/1‑250 (MQ=255)                      gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAAcc                                                                                                                                                                                                                   >  1:1052/1‑251 (MQ=255)                        gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                                                                                                                                                                                                 >  2:153126/1‑251 (MQ=255)                                  ccaccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTtc                                                                                                                                                                                                      <  2:573817/251‑1 (MQ=255)                                                  aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCAt                                                                                                                                                                                      <  1:311958/251‑1 (MQ=255)                                                  aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCAt                                                                                                                                                                                      <  2:264718/251‑1 (MQ=255)                                                    tttcGATCATCCGCAGGCGAGAGGCCAAAGCTCCACGCCGTGCGCGACGAGCACAATAAACTTCTACTTCGGCAGATAACCGTGGCGACTGACCATTCTCACCCAGGTAAAGCTGAGTCGCCTCCATCCCCCGCTGTGTGACTCCCTACGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAACTCTGCTCATGCCTTTTTGTCAAAA‑CCCCAACTTCCCTTCTACTCAAATTCCGCAtac                                                                                                                                                                                    >  2:185711/3‑249 (MQ=255)                                                         atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCAtt                                                                                                                                                                               >  1:386642/1‑251 (MQ=255)                                                         atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGGCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACC‑CCGTTCTCATCAAATTCCGCATCTTCAtt                                                                                                                                                                               <  1:540185/251‑1 (MQ=255)                                                                       agagGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCATCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                                                                                                                                                                                                                >  1:385024/1‑189 (MQ=255)                                                                         agGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTc                                                                                                                                                               <  1:10383/251‑1 (MQ=255)                                                                                        gCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTc                                                                                                                                                <  1:195832/251‑1 (MQ=255)                                                                                                cgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAt                                                                                                                                              <  1:128663/245‑1 (MQ=255)                                                                                                cgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAt                                                                                                                                              >  2:128663/1‑245 (MQ=255)                                                                                                      gcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCt                                                                                                                                  <  2:261372/251‑1 (MQ=255)                                                                                                       cAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                                                                                                                                                      >  2:566978/1‑214 (MQ=255)                                                                                                       cAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                                                                                                                                                      <  1:566978/214‑1 (MQ=255)                                                                                                          aaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTAccc                                                                                                                              <  2:184057/251‑1 (MQ=255)                                                                                                           aTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCg                                                                                                                             >  2:488733/1‑251 (MQ=255)                                                                                                             aaccttcgtcttCGGGCGTTAACAGTGGTGAATACCCATGCTCACCAAGGGAAAGCTGAGTCGGGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCAATGGGTTCCAATTCTACTAAGTCCTTTTCTTCAAAAACCTCAACT‑CCGTTCTTATCAAAGTCCGCATCTTCAGTTAAATCTTGAACTTTCATCGTATTCTGACTCAGAAGGTTGCttttggccgtg                                                                                                                           <  2:464824/248‑11 (MQ=255)                                                                                                             aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTg                                                                                                                           <  1:443732/251‑1 (MQ=255)                                                                                                              aaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTg                                                                                                                           >  1:320975/1‑250 (MQ=255)                                                                                                              aaCTTCTTCTTCGGCCGTTAACAGGGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCGCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGa                                                                                                                          <  2:115272/251‑1 (MQ=255)                                                                                                                  tcttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTTCCCGTGATcc                                                                                                                       <  2:166883/250‑1 (MQ=255)                                                                                                                    ttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCtt                                                                                                                    >  1:487674/1‑251 (MQ=255)                                                                                                                     tcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCtt                                                                                                                    <  2:278644/250‑1 (MQ=255)                                                                                                                       ttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGAc                                                                                                                 >  2:3813/1‑251 (MQ=255)                                                                                                                        tcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGAc                                                                                                                 <  1:548546/250‑1 (MQ=255)                                                                                                                         cGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACgg                                                                                                               >  1:480404/1‑251 (MQ=255)                                                                                                                         cGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACgg                                                                                                               >  1:447877/1‑251 (MQ=255)                                                                                                                           gCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGAAAAGCTGAGTCGCGACCAACACACGCTGTGTGGCCCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGCTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGa                                                                                                              >  1:216431/1‑250 (MQ=255)                                                                                                                                tAACAGTGGTGAATAACCAATCTCACCAAGGAAAACCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAttt                                                                                                        >  1:128934/1‑250 (MQ=255)                                                                                                                                 aaCAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCa                                                                                                       <  1:474990/251‑1 (MQ=255)                                                                                                                                     gtggggAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATGCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCGTCATTTAAATCATGAACTTTCAGCGTAGTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCa                                                                                                   <  2:370986/251‑1 (MQ=255)                                                                                                                                              aaCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTg                                                                                           <  1:338785/250‑1 (MQ=255)                                                                                                                                                ccAATCTCACCAAGGAAAAGCTGAGTCGCGTCGAACACACGCTGGGTGGATCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CAGTTCTCATCAAATTCCGCATCATCATTTAAATCATGAACTGTCGGCGGATTCTGACCCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACCTTCAAGCAAAAGCCTGGtt                                                                                        >  2:106910/1‑251 (MQ=255)                                                                                                                                                    tctcACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTccgc                                                                                    >  2:25146/1‑251 (MQ=255)                                                                                                                                                         ccAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAtt                                                                               >  1:355889/1‑251 (MQ=255)                                                                                                                                                                      gAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCgg                                                                    <  1:346116/249‑1 (MQ=255)                                                                                                                                                                            cgTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACG‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGATCTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAAc                                                            <  2:86250/251‑1 (MQ=255)                                                                                                                                                                             gTCCAACACAAGCTGTGTGGCTCCCTGCGA‑TACCAGATCCTCTTCGGCCAAATCGTTAGCACGGGCGTCCAATGCGACTAAGGCATCTTCGGCACAAGACTCAACT‑CCCTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAAATTTCAGCGTATTCTGCCTCATAAGGCGGCTCCTGCCCGTGATCCCTTGACGGAACATTCACGCAACAGCCTGGTTCCGCCGATTTATCGCTGCGGCCAATAAcc                                                           >  2:163639/1‑250 (MQ=255)                                                                                                                                                                                  acacacGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCg                                                      <  1:222896/251‑1 (MQ=255)                                                                                                                                                                                           gtgtgGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTAc                                               <  2:355889/249‑1 (MQ=255)                                                                                                                                                                                           gtgtgGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACg                                              <  1:248039/250‑1 (MQ=255)                                                                                                                                                                                               ggCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGAtt                                          <  2:18358/250‑1 (MQ=255)                                                                                                                                                                                                ccTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTCCCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTc                                        <  2:244733/250‑1 (MQ=255)                                                                                                                                                                                                  tCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTccc                                      >  2:221743/1‑251 (MQ=255)                                                                                                                                                                                                          a‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAAc                              >  2:500192/1‑251 (MQ=255)                                                                                                                                                                                                             aaCAGTTCCTCTTCGGCCAAAACGTTATCACTGGGGTCCAATTCTACTAACGCCTTTTCGTAAAAAACCTCAACA‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCTTGATCCCTTGCCGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGAGTTTACGGATTTCCCCTTGTaccg                             >  2:383021/1‑247 (MQ=255)                                                                                                                                                                                                                         tCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATg                 <  1:443819/250‑1 (MQ=255)                                                                                                                                                                                                                           ggCCAAATCGTTATCACTGGGTTCCAATTCTAGTAAGGCCTTTTCGTCAACAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATGTAAATCATGAACTTTCAGCGTATTCTGACTCATACGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGTTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCGTGTAACGAATTTCAAAATGCaa              <  2:443884/251‑1 (MQ=255)                                                                                                                                                                                                                           ggCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCaa              <  1:334591/251‑1 (MQ=255)                                                                                                                                                                                                                                   cGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCTTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAATCGTGtt       <  2:142437/250‑1 (MQ=255)                                                                                                                                                                                                                                      tttCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCGACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCACGCGTGTTGaa    <  2:216431/248‑1 (MQ=255)                                                                                                                                                                                                                                       aTCACTGGGTTCCAATTCTCCTAAGGCCTTTTCGTCAAAAACCTCAACT‑CTGTTCTCATAAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCAGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTGCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGatc   <  2:274977/250‑3 (MQ=255)                                                                                                                                                                                                                                       aTCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACt  <  1:92090/251‑1 (MQ=255)                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                              GACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGA‑TAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACT‑CCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACT  >  NC_000913/2867212‑2867691

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