Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,518,356	A→G	intergenic (‑121/‑16)	fecI ← / → ykiE	RNA polymerase sigma factor FecI/protein YkiE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,518,356	0	A	G	100.0%	209.8 / NA	58	intergenic (‑121/‑16)	fecI/ykiE	RNA polymerase sigma factor FecI/protein YkiE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/24);  total (34/24)

AATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAA  >  NC_000913/4518113‑4518591                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            aaTGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAaggta                                                                                                                                                                                                                                      <  2:289821/251‑1 (MQ=255)      cagcTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTGCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGGCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAggtaggta                                                                                                                                                                                                                                  <  1:307805/247‑1 (MQ=255)                                ttGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACAGTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCt                                                                                                                                                                                                        <  1:680559/250‑1 (MQ=255)                                      gTCAGCCAGCTTTTCAACGAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCAtgt                                                                                                                                                                                                 >  2:601254/1‑251 (MQ=255)                                         agccagcTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTc                                                                                                                                                                                              <  1:10470/251‑1 (MQ=255)                                         agccagcTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTc                                                                                                                                                                                              >  1:330279/1‑251 (MQ=255)                                            cagcTTTTCAACCAGCCGGGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATGACAATTGGAAATAAAATTGTTTCCAATAGACGTTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGt                                                                                                                                                                                           <  2:703772/251‑1 (MQ=255)                                              gcTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGAGGGAGCCTCTGCTTCTGGCATGTGTCGGTCa                                                                                                                                                                                         <  1:561662/251‑1 (MQ=255)                                                         aGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCCTATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGTGCCTCTGCTTCTGGCATGTGGCGTTCAGAATGACTCat                                                                                                                                                                              >  2:4564/1‑251 (MQ=255)                                                         aGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCat                                                                                                                                                                              <  1:258036/251‑1 (MQ=255)                                                                  gtgtgCCATAAAGCGACTTGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAAGTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTTGGTATAAAAGGGGATGGAGCATCTGCTTCTGGCATGTGTCGGACAGAATGACTCATGATGTGGTc                                                                                                                                                                     <  2:450631/251‑1 (MQ=255)                                                                  gtgtgCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGGTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTc                                                                                                                                                                     <  1:256269/251‑1 (MQ=255)                                                                    gtgCCATAAAGCGACTAGAACGTTAATGAAGCTGTGGGAGTTGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACAGTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCt                                                                                                                                                                    <  2:223845/250‑1 (MQ=255)                                                                         aTAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTAct                                                                                                                                                              >  1:414147‑M1/1‑249 (MQ=255)                                                                             aGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGGATAAAATGGGATGGAGCCTCTGCTTCTGGCATGGGTCGGTCCGAAGGACTCATGATGTGGTCTGCTActca                                                                                                                                                            >  2:64050‑M1/1‑245 (MQ=255)                                                                                gACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtg                                                                                                                                                       >  1:213466‑M1/1‑242 (MQ=255)                                                                                      aaCGTTAAGGAAGCTGTGGTAGTGGCGCGGTCCGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTAGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggaggaattg                                                                                                                                                 >  2:701885‑M1/1‑236 (MQ=255)                                                                                         gTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgt                                                                                                                                              >  1:662823‑M1/1‑233 (MQ=255)                                                                                                 aGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGGTGTGGTCTGCTActcaaggtggaattgtgtcgcagtat                                                                                                                                      <  2:587691‑M1/251‑27 (MQ=255)                                                                                                  gCTGTGGTAGTGGCGCGGGCAGACATGCGGAGTTCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACAGGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgggtcgcagtata                                                                                                                                     <  2:711595‑M1/251‑28 (MQ=255)                                                                                                     gtgGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTTAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtatac                                                                                                                                    <  2:376954‑M1/249‑29 (MQ=255)                                                                                                     gtgGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGATGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtatacc                                                                                                                                   <  1:470875‑M1/250‑30 (MQ=255)                                                                                                     gtgGTAGTGGCGCGGTCAAACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtatac                                                                                                                                    <  2:163925‑M1/249‑29 (MQ=255)                                                                                                            tGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaa                                                                                                                            >  1:715625‑M1/1‑214 (MQ=255)                                                                                                                             cGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctt                                                                                                          >  1:250303‑M1/1‑197 (MQ=255)                                                                                                                              ggAGTGAATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTTCAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGGAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggttgaattgggtcgcagtataccgccagaaacggggattaatactcttt                                                                                                         <  2:77939‑M1/251‑56 (MQ=255)                                                                                                                                          aagataattATCAAGTAGTCATAATAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcacggtggaattgtgtcgcagtataccgccagaaacggggataaatgctctttacctgttacccg                                                                                             >  2:383000‑M1/5‑184 (MQ=255)                                                                                                                                             ttAATTATCACGTAGTCACATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgaagtataccgccagaaacggggattaatactctttacctgttacccggct                                                                                          >  2:598081‑M1/1‑181 (MQ=255)                                                                                                                                                   ctcACGTAGTCAGATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTATAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacgggtattaatactctttacctgttacccggctacggt                                                                                     >  2:707354‑M1/2‑175 (MQ=255)                                                                                                                                                   aTCACGTAGTCATATTAATATGAGAATGGTNATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGCGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggt                                                                                     <  2:339629‑M1/250‑76 (MQ=255)                                                                                                                                                     cACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagacacggggattaatactctttacctgttacccggctacggtagc                                                                                  >  2:229395‑M1/1‑173 (MQ=255)                                                                                                                                                       cGTAGTCATATTAATATGAGAATGGAAATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagcaacggggattaatactctttacctgttacccggctacggtatct                                                                                 >  1:195019‑M1/1‑171 (MQ=255)                                                                                                                                                        gTAGTCATATTAATATGAGAATGGTTATCATTAAAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCAActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttaccaggctacggtagctat                                                                               >  1:490296‑M1/1‑170 (MQ=255)                                                                                                                                                          aGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAAAATGTTGTTTTTCTAAGTGTTGTAAGGAAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtctaccgccagaaacggggattaatcctntttacctgttacccggctacggtagctatcc                                                                             >  2:332719‑M1/1‑168 (MQ=255)                                                                                                                                                              atatTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagt                                                                         >  2:506730‑M1/1‑164 (MQ=255)                                                                                                                                                                atTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttg                                                                       >  1:317786‑M1/1‑162 (MQ=255)                                                                                                                                                                atTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagtt                                                                        <  2:96331‑M1/250‑89 (MQ=255)                                                                                                                                                                    atatGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTACGGTAGGTATAAAATGGGATGATGCCTCTGCTTCTGGCATGTGTCGGTCAGCATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtctaccgccagaaacggggattaatactctttacctgttacccggcgacggtagctatcccagttgaagt                                                                   >  2:169514‑M1/1‑158 (MQ=255)                                                                                                                                                                       tGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatc                                                                 >  1:46226‑M1/1‑155 (MQ=255)                                                                                                                                                                         agaATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacc                                                              <  1:588422‑M1/251‑99 (MQ=255)                                                                                                                                                                              ggTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccaca                                                          <  1:601254‑M1/250‑103 (MQ=255)                                                                                                                                                                                 tATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTACGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgagtcgcagtataccgccagaaacggggattaatacactttacctgttacccggctacggtagctatcccagttgaagataacccccaaact                                                      >  1:426895‑M1/1‑145 (MQ=255)                                                                                                                                                                                 tATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaact                                                      >  2:104425‑M1/1‑145 (MQ=255)                                                                                                                                                                                 tATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaact                                                      >  1:255247‑M1/1‑145 (MQ=255)                                                                                                                                                                                  aTCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTATGGTAGGCATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggcggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatccc                                                                            >  1:307317‑M1/1‑144 (MQ=255)                                                                                                                                                                                  aTCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTATTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtagaccgccagaaacggggattaatactctttacctgttaccctgctacggtagctatccc                                                                            <  2:307317‑M1/228‑85 (MQ=255)                                                                                                                                                                                    cATTACAATTGGAAATAAAATTGTTTACCCTAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGAATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaagggggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagctaacccacaaactgtt                                                   >  1:566779‑M1/1‑142 (MQ=255)                                                                                                                                                                                     aTTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgtt                                                   >  1:677760‑M1/1‑141 (MQ=255)                                                                                                                                                                                         cAATTGGAAATAAAATTGTTTCCAATAGACATTCTTAACATGTTGTTTTTCTAAGTGTTGTCAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCCTGTGTCGGACAGAATGACTCATGATGTGGTCTGCTCctcaaggtggaattgtgtcgcagtataccgccagaaacggggattaataatctttacctgatacccggctacggtagctatcacagttgcagatcacccacaaactgttacca                                               >  2:608724‑M1/1‑137 (MQ=255)                                                                                                                                                                                                 aaTAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcat                                       >  1:301456‑M1/1‑129 (MQ=255)                                                                                                                                                                                                    caaatTGTTTCCAATAGACATTTTTACCATGTTGTTTTTCTCAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActccaggtggaattgtgtcgcagtatacagccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggctcatacgg                                   >  2:306804‑M1/2‑126 (MQ=255)                                                                                                                                                                                                       aTTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtcc                                <  1:395787‑M1/251‑129 (MQ=255)                                                                                                                                                                                                            ttCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcatggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcac                           >  2:235336‑M1/1‑118 (MQ=255)                                                                                                                                                                                                                   aGACATTTTTAACATGTTGTTTTTCTCAGTGTTGTAAGGTAGGTATAAAAAGGGATGGAGCCTCTGCTTCTGGAATGTGTCGGTCAGAATGACTCATGCTGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaaatgttacccaggcgcatacggtccatcacgcgttcg                    >  2:99823‑M1/1‑111 (MQ=255)                                                                                                                                                                                                                      cATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgccc                 <  1:525053‑M1/251‑144 (MQ=255)                                                                                                                                                                                                                           ttAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGGATAAAATGGAATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaagggggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccccaaactgttacccaggcgcatacggtccatcacgcggtcgcccagca             >  2:445092‑M1/1‑103 (MQ=255)                                                                                                                                                                                                                               cATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttg        <  2:414147‑M1/251‑153 (MQ=255)                                                                                                                                                                                                                                     gttTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatc  >  1:18451‑M1/1‑93 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            AATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAA  >  NC_000913/4518113‑4518591

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