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%	208.7 / 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 (30/28);  total (30/28)

TCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTG  >  NC_000913/4518120‑4518599                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    tCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGtata                                                                                                                                                                                                                                        >  2:328555/1‑250 (MQ=255)        cAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCCTATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATggg                                                                                                                                                                                                                                 >  2:185050/1‑250 (MQ=255)         aaaaGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGAt                                                                                                                                                                                                                               >  2:438043/1‑251 (MQ=255)          aaaGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATg                                                                                                                                                                                                                              <  1:567562/251‑1 (MQ=255)           aaGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATgg                                                                                                                                                                                                                             <  1:40680/251‑1 (MQ=255)           aaGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATgg                                                                                                                                                                                                                             >  1:193563/1‑251 (MQ=255)                                                    ccGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAACGTTAATTATCACGTAGTCCTATTAATATGAGAATGGTTATCATTACAATTGGAACTAAAATTGATTCCAATAGATATTTTTACCATGTTGTTTTTTTAAGAGTTGTAAGGTAGGTCTAAAATGGGATGGAGCCTCTTCTTCTGGCATGTGTCGGTCAGAATGACTTatg                                                                                                                                                                                     >  2:18426/1‑250 (MQ=255)                                                     cGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCatgat                                                                                                                                                                                   <  2:250329/251‑1 (MQ=255)                                                           gtgtgCCATAAAGAGACTCGAACGTTAAGGAAGGTGTTGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCAAGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAAATGTTTCCAATAGACAGTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGTGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACACATGATGGGGt                                                                                                                                                                              <  2:271123/250‑1 (MQ=255)                                                                  aTAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTAct                                                                                                                                                                      <  1:469125‑M1/251‑3 (MQ=255)                                                                               aaCGTTAACGAAGCTGTGGTATTGGTGCGGTCAGACATGCGGGGTGCATCAACCGTTAATTATCACGTAGTCATATTAATCTGAGAATGGTTATCATTCCAATTGGAAATAAAATTGATTCAAATAGACATTTTTAACATGTTGTTTTTATAAGTGTTGTATGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTTGCATGTGTCGGTCAGCGTCACTCATAATGTgggtggctcctacagggggaattt                                                                                                                                                         >  2:653134/1‑229 (MQ=255)                                                                                              gtgGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccg                                                                                                                                          >  1:471207‑M1/1‑221 (MQ=255)                                                                                              gtgGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtatacc                                                                                                                                           <  1:438043‑M1/250‑30 (MQ=255)                                                                                                       gcgcGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacg                                                                                                                                 <  1:328555‑M1/251‑40 (MQ=255)                                                                                                         gcgGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggg                                                                                                                               >  1:147061‑M1/1‑210 (MQ=255)                                                                                                               aGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggatta                                                                                                                          >  1:535454‑M1/1‑204 (MQ=255)                                                                                                                 aCATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatt                                                                                                                       >  1:23361‑M1/1‑202 (MQ=255)                                                                                                                       ggAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTTTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactcttt                                                                                                                 <  2:734818‑M1/251‑56 (MQ=255)                                                                                                                         aGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttac                                                                                                               <  1:348656‑M1/251‑58 (MQ=255)                                                                                                                          gTGCATCAAAAGTTAATTATCGCGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCGActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacc                                                                                                              <  1:252216‑M1/251‑59 (MQ=255)                                                                                                                                 aaaaGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacc                                                                                                       <  1:105401‑M1/251‑66 (MQ=255)                                                                                                                                 aaaaGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacc                                                                                                       >  1:567374‑M1/1‑186 (MQ=255)                                                                                                                                   aaGTTAATTATCACGTGTTCATATTCATATGAGACTGGTTATCATTACAATTGGAAATAAAATTGTTTCCAAAAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTACGATAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGAGGTGGTCTGCTAcgcaaggtggtattgtgtcgcagtataccgccagaaacggggattaatactctttcactgttacccg                                                                                                     >  2:46052‑M1/1‑184 (MQ=255)                                                                                                                                   aaGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccg                                                                                                     <  1:414300‑M1/251‑68 (MQ=255)                                                                                                                                     gTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaata                                                                                                                       >  2:568472‑M1/1‑182 (MQ=255)                                                                                                                                     gTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaata                                                                                                                       <  1:568472‑M1/231‑50 (MQ=255)                                                                                                                                        aaTTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggcta                                                                                                 <  2:567374‑M1/250‑72 (MQ=255)                                                                                                                                        aaTTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggcta                                                                                                 <  1:676193‑M1/250‑72 (MQ=255)                                                                                                                                            aTCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggta                                                                                            <  1:310104‑M1/251‑77 (MQ=255)                                                                                                                                            aTCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGGTTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatcctctttacctgttacccggctacggta                                                                                            <  2:111647‑M1/251‑77 (MQ=255)                                                                                                                                                cGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTATGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagct                                                                                         >  1:217740‑M1/1‑171 (MQ=255)                                                                                                                                                  tAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGAATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGACGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatc                                                                                      >  2:480863‑M1/1‑169 (MQ=255)                                                                                                                                                       atatTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctattccagt                                                                                 <  1:228759‑M1/251‑88 (MQ=255)                                                                                                                                                         atTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttg                                                                               >  1:188127‑M1/1‑162 (MQ=255)                                                                                                                                                             atatGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttga                                                                              >  1:755912‑M1/1‑158 (MQ=255)                                                                                                                                                             atatGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCAGGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaaga                                                                           <  1:577419‑M1/251‑94 (MQ=255)                                                                                                                                                                  agaAGGGTTATCATTACAATTGGAAATAAAAGTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCGAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGGGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtatgccgccagaaacggggattaatactctttacctgttagccggctacggtagctatcccagttgaagatcacc                                                                      <  2:23361‑M1/251‑99 (MQ=255)                                                                                                                                                                   gaATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcaccc                                                                     >  1:21204‑M1/1‑152 (MQ=255)                                                                                                                                                                         ttATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAAAATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGAATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaac                                                               >  2:195291‑M1/1‑146 (MQ=255)                                                                                                                                                                              aTTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaatgtggaattgtgtcgcagtataccgccagaaacgtggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgtta                                                          <  2:609380‑M1/251‑111 (MQ=255)                                                                                                                                                                              aTTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttc                                                          >  1:220577‑M1/1‑141 (MQ=255)                                                                                                                                                                               ttACAATTGGAAATAAAATTGTTTCAAATAGCCATTTTTAACATGTTGATTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGCGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcacggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttac                                                         >  1:263551‑M1/1‑140 (MQ=255)                                                                                                                                                                                   aattgtgaataaaattGTTTCCAATAGACATNTTTAACATGTTGTTTTTCTAAGTGTTGTAATGTAGGTATAAAATGGGATGGGGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTTGTCTGCTActcaaggtggaattttgtcgcagtataccgccagaaacggggattaatactctttccctgttacccggctacgggagctatcccaggtgaaggtctcccacaaactgttgccctg                                                     <  2:338565‑M1/244‑116 (MQ=255)                                                                                                                                                                                   aaTTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttaccc                                                       >  1:642602‑M1/1‑136 (MQ=255)                                                                                                                                                                                           aTAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatac                                             <  1:239179‑M1/251‑124 (MQ=255)                                                                                                                                                                                                 ttGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacncacaaactgttacccaggcgcatacggtcc                                        >  2:338952‑M1/1‑122 (MQ=255)                                                                                                                                                                                                  tatttccACTAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGGCAGAATGACTCCTGATGTGGTCTGCTActcaaggtggaattgtgacgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcaaccacaaactgttacccaggcgcatactgtccat                                      >  2:559213‑M1/3‑121 (MQ=255)                                                                                                                                                                                                  tGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATTATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcccccacaaactgttacccaggcgcatacggtcca                                       <  2:193563‑M1/250‑130 (MQ=255)                                                                                                                                                                                                       ccAATAGACATTTTTATTATGTTGTTTTTCTAAGTGTTGTAAGGGAGGGATAAAATGGGATGGAGACTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTCctcaaggtggaattgcgtcgcagtataccgccagaatcggggattaatactctatccctcttacccggctacggtagctatcccagttcaagatcacccacaatctgttacccaggctcatacggtccatcactc                                 <  2:217740‑M1/251‑136 (MQ=255)                                                                                                                                                                                                       ccAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgc                                 >  1:668268‑M1/1‑116 (MQ=255)                                                                                                                                                                                                            aGACATTTTTATCATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttc                             <  1:559213‑M1/250‑140 (MQ=255)                                                                                                                                                                                                            aGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATAGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACCCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatctcagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcg                            >  1:261646‑M1/1‑111 (MQ=255)                                                                                                                                                                                                                        cATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttg                >  1:546516‑M1/1‑99 (MQ=255)                                                                                                                                                                                                                         aTGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGAGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttgg               >  1:283338‑M1/1‑98 (MQ=255)                                                                                                                                                                                                                              gttTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggcttaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcccccagcaacttggtcctc          >  2:680720‑M1/1‑93 (MQ=255)                                                                                                                                                                                                                              gttTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatc          <  1:185050‑M1/251‑159 (MQ=255)                                                                                                                                                                                                                                ttttCTAATTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCGTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgtttcccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatct         <  2:618948‑M1/250‑160 (MQ=255)                                                                                                                                                                                                                                       aGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccat  >  1:110961‑M1/1‑84 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    TCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTG  >  NC_000913/4518120‑4518599

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: ‑