Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,867,455	G→A	Q33* (CAG→TAG)	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,455	0	G	A	100.0%	162.6 / NA	49	Q33* (CAG→TAG)	rpoS	RNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/23);  total (26/23)

GATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCC  >  NC_000913/2867218‑2867697                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   gATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGcc                                                                                                                                                                                                                                       <  1:411134/251‑1 (MQ=255)  aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGCGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCt                                                                                                                                                                                                                                      <  2:183249/251‑1 (MQ=255)                  gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                                                                                                                                                                                                      >  1:476120/1‑251 (MQ=255)                              accaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtctc                                                                                                                                                                                                          >  2:222631/1‑251 (MQ=255)                                   aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaa                                                                                                                                                                                                     <  1:111340/251‑1 (MQ=255)                                     gCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAt                                                                                                                                                                                                   >  2:77180/1‑251 (MQ=255)                                       aaGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAACGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCCAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTc                                                                                                                                                                                                 >  2:114785/1‑251 (MQ=255)                                       aaGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTACGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTc                                                                                                                                                                                                 >  1:103697/1‑251 (MQ=255)                                        aGTTACTCTCGATCATCCGGCAGCGAGAGGGGAACTCTCCACGCAGTGCGCGACGCGAAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATCACCACTCTCACCAAGGTAAAGCTGCGTCGCGACCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTACTCTTCGGCCAACTCGTTATAACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAACACCTCAAATCCGTTCTCATCAAATTcc                                                                                                                                                                                                >  1:416563/1‑251 (MQ=255)                                             ctctCGATCATCCGGCGGCGCGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAACTCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAAATCCGTTCTCATAAAATTCCGCATc                                                                                                                                                                                           >  2:381430/1‑251 (MQ=255)                                                 cGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCa                                                                                                                                                                                       <  1:482164/251‑1 (MQ=255)                                                 cGATCAACCGGCGGCGAGAGGCGACATCTCCATGCAGTGCGCGACGCGCAAAATAAACGTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCagca                                                                                                                                                                                          <  2:499712/248‑4 (MQ=255)                                                  gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAt                                                                                                                                                                                      <  1:307309/251‑1 (MQ=255)                                                     catcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTa                                                                                                                                                                                   <  1:110022/251‑1 (MQ=255)                                                         cggcggCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCCTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                               <  1:2394/251‑1 (MQ=255)                                                             ggcgAGCGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGACAAAACCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                                                                                                                                                            >  2:368805/1‑250 (MQ=255)                                                                gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACt                                                                                                                                                                        >  1:181712/1‑251 (MQ=255)                                                                              ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTAtt                                                                                                                                                             <  1:381430/248‑1 (MQ=255)                                                                              ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTg                                                                                                                                                          >  1:123502/1‑251 (MQ=255)                                                                              ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTGCAGCGGAtt                                                                                                                                                             <  2:388645/248‑1 (MQ=255)                                                                              ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCAGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGtcttctt                                                                                                                                                          >  2:137101/1‑250 (MQ=255)                                                                                 cGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACt                                                                                                                                                       <  2:143084/251‑1 (MQ=255)                                                                                     gtgaacgaaccgcACAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAAACAATCTCCCCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGCCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTACGGCCTTTTCGTCAAAAACATCAACACCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTATGATTCAt                                                                                                                                                    >  2:99575/11‑250 (MQ=255)                                                                                         gcgACGCGCAAAATAAACTTCTTCATCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCGCTTCTGCCAAATCGTTATCACTGTGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTTAACTCCGTTCTCAACAAATTCCGCATCTTCATTTAAAGCATGAACTTTCATCGTATTCTGACTCATAAGGt                                                                                                                                               <  2:409326/251‑1 (MQ=255)                                                                                            aCGCGCAAAATAAACTTCTTCTTCCGCCGTTAACAGTGGTGCATAACCAATCTCACCAAGGTAAAGCTGAGTCGCATCCACCACACGCTGTGTGGCTCCCTGCGATAACAGTTACTCTTCGGCCAAATCGTTATCACTGGGATCTAATTCTACTAATACCTTTTCGTCAAAAACCTCAACTCCGTTCTCATAAAATTCCGCATCTTCATTTAAATCATGCCCTTTCAGCGTATTCTGACTCTTACGGTGGc                                                                                                                                            >  2:35831/1‑251 (MQ=255)                                                                                              gcacAAAATAAACTTCTTCTTCGGCCGTTAACTGAGGTGAATATCCGATCTCACCAAGGTAAAGCTGATGCGCGTCCTACACCCGCCATGTGGTTTCCTGCGATAACAGTTCCACTTCGGCCAAATCTGTATCACTGGGTTCCTATTCTACCAAGGCCGTGTAGTCAAAATCCATAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAAGCATGAACGGCCAGCGTGTTCTGACGCATAAGGGGGCta                                                                                                                                          <  2:416563/248‑2 (MQ=255)                                                                                                     aTAAACTTCTTCTTCGGCCGTTAACAGTGGCGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGAGTCTTCTGACTCATAAGGTGGCTCCTACCCg                                                                                                                                   >  2:372038/1‑251 (MQ=255)                                                                                                       aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGt                                                                                                                                  <  1:434305/250‑1 (MQ=255)                                                                                                                           taCAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCACAAACCTCACCTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTc                                                                                                              >  2:122495/2‑250 (MQ=255)                                                                                                                            aCAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCCACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTTGCTCCTACCCGTGATCCCTTGACGGAACATTCa                                                                                                             >  1:340021/1‑250 (MQ=255)                                                                                                                              aGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGc                                                                                                          <  1:77180/251‑1 (MQ=255)                                                                                                                                                            aGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAATCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATTAACTTTCAGTGTAGTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCGTGACGGAACATTCAAGCAAAAGCTTGGTTCCGCCGATTTATCGCTga                                                                            <  2:63115/251‑2 (MQ=255)                                                                                                                                                                 aGTCGCGTCCCACTCATGCTGTGTGGCGCCCTGCTATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTGCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCaa                                                                       <  2:470356/251‑1 (MQ=255)                                                                                                                                                                      cgTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAAc                                                                  >  1:30311/1‑251 (MQ=255)                                                                                                                                                                       gTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGAGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACg                                                                 >  2:512281/1‑251 (MQ=255)                                                                                                                                                                          cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCCAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGc                                                                >  2:110187/1‑249 (MQ=255)                                                                                                                                                                            acacacGCTGTGTGACTACCTGCGATAACCGTCTCTCTTCCGCCAAATCGTTATCACCGGGTTCCTATTCTTCTAAGCCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGACCTTTCAGCGTATTCTGACTCATAAGGTGGCTACTACCCGTGATCCCTATACGGAACATTCTAGCAAAAGCCTCGTTCCGCAGATTTATCGCTGCGGCAAATAACGCAGCg                                                            >  2:491763/1‑251 (MQ=255)                                                                                                                                                                               cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATCCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGt                                                         >  2:308418/1‑251 (MQ=255)                                                                                                                                                                                 cGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGttt                                                       >  1:291157/1‑251 (MQ=255)                                                                                                                                                                                       gtgGCTCCCTTCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGa                                                  <  1:201020/250‑1 (MQ=255)                                                                                                                                                                                                 gCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGTGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTg                                        <  1:363208/250‑1 (MQ=255)                                                                                                                                                                                                    aTAACATTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCACCTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTGACGGATTTCCCCTTGTAAc                                    <  2:47277/251‑1 (MQ=255)                                                                                                                                                                                                    aTAACAGTTCCCCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAATCCTCAACTCCGTGCTCATCAAATTCCGCATCTTCATTTAAATCATGAACGTTCAGCGTATTCTGACTCATAAGGGGGTTCCTTCCCGGTATCCCTTGACGGAACATTCATTCGGAAGCCTGGTTCCGCCGAGTTATCGCTGCGGCAAATAACGTAGCGGGGTTACGGATTTCCGCTTGTaa                                     <  2:591938/250‑1 (MQ=255)                                                                                                                                                                                                                    ggCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCa                     >  2:198206/1‑250 (MQ=255)                                                                                                                                                                                                                         aaacGTTATCACTGGTTTCCTATTCTACTAATGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATAATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATGCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCgtg               <  2:327383/248‑1 (MQ=255)                                                                                                                                                                                                                         aaTCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCgtg               >  1:440373/1‑251 (MQ=255)                                                                                                                                                                                                                           tCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGtt             >  2:161179/1‑251 (MQ=255)                                                                                                                                                                                                                                 tCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTg       >  1:445610/1‑251 (MQ=255)                                                                                                                                                                                                                                      gggTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATGCCGCATCTGCATTTAAATCAGGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTGCCCGGTATCCCTTGCCGGAACATTCAAGCAAAAGCCTTGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTcc  <  1:354306/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   GATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCC  >  NC_000913/2867218‑2867697

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