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%	276.3 / NA	72	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 (39/33);  total (39/33)

GCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCC  >  NC_000913/4518111‑4518601                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      gCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAagg                                                                                                                                                                                                                                                  <  2:876243/251‑1 (MQ=255)  cAATGTCATCTGCACCAAAAGCAGACTGGAGTTTGCGCGTCAGGCAGCTGTTCAACCAGCCGTGATGTGGGCCATTAAGCTACGCGAACGTTAAGGAAGCTGTGGTTGTGGCGCGGTCAGACCTGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTCCAATTGGAAATAAAATTGTGTCCAATAGACAGTTTTAACATGTTGTTTTTCTAAGTGTTGTAaggt                                                                                                                                                                                                                                                 <  2:439829/251‑1 (MQ=255)        cagcTGCATCAAAAGCAGACTGGATTTTGCTCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCAGAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTTGCGCGGTCAGACATGCGGAGTGCATCAAAATTTAATTATCACGTTGTCATATTAATAGGAGAATGGTTATCATTACAATTGGAAATAAAATTGGTACCAGTAGACATTTTTAACATGTTGTTTTTCTAAGTGGTGTAAggtaggtag                                                                                                                                                                                                                                           <  2:341255/248‑2 (MQ=255)        cATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTTGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGtat                                                                                                                                                                                                                                           <  1:252707/251‑1 (MQ=255)              cATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATg                                                                                                                                                                                                                                     <  2:65539/251‑1 (MQ=255)                 cAAAAACAGATTGGAGTTTGAGCGTCAGCCAGCTTTTCAACAAGCCGTGATGTGTGCAATAAAGCGACTCGAACGTTAAGGACGCTGTGGAAGTGGCGCGGCCAGACATCCGGAGGGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATCAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATggg                                                                                                                                                                                                                                   >  2:235518/1‑250 (MQ=255)                  aaaaGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTATAAGTGTTGTAAGGTAGGTATAAAATGGGAt                                                                                                                                                                                                                                 >  1:404221/1‑251 (MQ=255)                   aaaGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGACAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGTAAATAAAATTGTTTCCAATAGAAATTTTTAACATGTTGTTTTTCTAAGTGTTGTACGGTAGGTATAAAATGGGATg                                                                                                                                                                                                                                >  2:827000/1‑251 (MQ=255)                      gCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATgg                                                                                                                                                                                                                               >  1:325417/1‑249 (MQ=255)                              gAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGCTGTGTGCCATCAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGAGGCGCGGTCAGACATGCGGAGTGAATCAAAAGATAATTATCACGTAGTCATATTAATATGCGAATGGTTATCATTACAATTGGTAATACAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTATAAGTGTTGTAGGGTAGGTATAAAATGGGATGGAGCCTCTGCt                                                                                                                                                                                                                     >  2:614646/1‑251 (MQ=255)                                 tttGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGCACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCt                                                                                                                                                                                                                  >  1:629844/1‑251 (MQ=255)                                     cgcgTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGc                                                                                                                                                                                                               >  1:828165/1‑250 (MQ=255)                                     cgcgTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCa                                                                                                                                                                                                              <  2:104721/251‑1 (MQ=255)                                           agccagcTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTc                                                                                                                                                                                                        <  1:786937/251‑1 (MQ=255)                                                 cTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTc                                                                                                                                                                                                    >  1:753704/1‑249 (MQ=255)                                                      cAACCAGCCGTGATGTGTGCCCTAAAGCGAGTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTGCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATg                                                                                                                                                                                              <  2:455899/250‑1 (MQ=255)                                                        aCCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACt                                                                                                                                                                                           <  1:438566/251‑1 (MQ=255)                                                         ccAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTACGGAAGCTGGGGGAGAGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTACTTATCAAGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGCATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTc                                                                                                                                                                                          >  2:278086/1‑251 (MQ=255)                                                              cGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCatgat                                                                                                                                                                                     <  1:40599/251‑1 (MQ=255)                                                                          cATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTAc                                                                                                                                                                         >  1:772935‑M1/1‑250 (MQ=255)                                                                                     tCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtgga                                                                                                                                                               >  2:411339‑M1/1‑239 (MQ=255)                                                                                       gAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaatt                                                                                                                                                            >  2:634235‑M1/1‑237 (MQ=255)                                                                                             tAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtc                                                                                                                                                       <  1:581816‑M1/250‑20 (MQ=255)                                                                                               atgaaGCTGTGGTAGTGGAGCGGTCAGACATGAGGAGAGCATCAAAAGTTAATTATCACGAAGTCATATTAATCTGAGAATGGTTATCATTAAAATTGGAAATAAAATTGATTCCCATAGACAGTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATACAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTAcgcaagggggtattgggtcgct                                                                                                                                                    >  2:507479‑M1/3‑229 (MQ=255)                                                                                                    gCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtat                                                                                                                                                >  1:59809‑M1/1‑224 (MQ=255)                                                                                                       gtgGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtattccg                                                                                                                                            <  2:635974‑M1/251‑31 (MQ=255)                                                                                                        tgGTAGTGGCGCGGTCAGCCATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgc                                                                                                                                           >  2:227354‑M1/1‑220 (MQ=255)                                                                                                            aGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccag                                                                                                                                        <  1:411339‑M1/250‑35 (MQ=255)                                                                                                            aGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTAcgcaaggtggaattgtgtcgcagtatacggccag                                                                                                                                        <  1:601715‑M1/250‑35 (MQ=255)                                                                                                               ggCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaac                                                                                                                                    <  1:459641‑M1/251‑39 (MQ=255)                                                                                                                     gTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTGATCATTACAATTGGAAATAAAATGGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggag                                                                                                                              <  2:504224‑M1/251‑45 (MQ=255)                                                                                                                       cAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtgtaattgtgtcgcagtataccgccagaaacggggatt                                                                                                                             <  1:614646‑M1/250‑46 (MQ=255)                                                                                                                                 gAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTACGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGAGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtatcccgccagaaacggggcttaatactcttt                                                                                                                   >  2:819578‑M1/1‑195 (MQ=255)                                                                                                                                 gAGTGCATCAAAAGATAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttc                                                                                                                  >  2:120860‑M1/1‑195 (MQ=255)                                                                                                                                  aGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttac                                                                                                                 <  2:649262‑M1/251‑58 (MQ=255)                                                                                                                                  aGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttac                                                                                                                 >  2:882274‑M1/1‑194 (MQ=255)                                                                                                                                      cATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgt                                                                                                             <  1:356985‑M1/251‑62 (MQ=255)                                                                                                                                       aTCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgtt                                                                                                            <  1:303275‑M1/251‑63 (MQ=255)                                                                                                                                         cAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttac                                                                                                          <  1:120385‑M1/251‑65 (MQ=255)                                                                                                                                         cAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttac                                                                                                          >  1:307806‑M1/1‑187 (MQ=255)                                                                                                                                          aaaaGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacc                                                                                                         <  2:591747‑M1/251‑66 (MQ=255)                                                                                                                                                   ttATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacatgttacccggctacgg                                                                                                >  1:723125‑M1/1‑177 (MQ=255)                                                                                                                                                       cACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagc                                                                                            >  2:523515‑M1/1‑173 (MQ=255)                                                                                                                                                        aCGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagct                                                                                           <  1:535713‑M1/251‑80 (MQ=255)                                                                                                                                                             gTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatccc                                                                                      <  1:201702‑M1/251‑85 (MQ=255)                                                                                                                                                                  atTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagtt                                                                                  <  1:654375‑M1/250‑89 (MQ=255)                                                                                                                                                                    tAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtatctatcccagttga                                                                                >  2:863550‑M1/1‑160 (MQ=255)                                                                                                                                                                           agaATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgggtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtaactatcccagttgaagatcacc                                                                        >  2:530356‑M1/1‑153 (MQ=255)                                                                                                                                                                            gaATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcaccc                                                                       >  1:508906‑M1/1‑152 (MQ=255)                                                                                                                                                                                 gTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaa                                                                   >  1:309357‑M1/1‑147 (MQ=255)                                                                                                                                                                                      cATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagtt                                                                                  <  1:580691‑M1/230‑89 (MQ=255)                                                                                                                                                                                      cATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagtt                                                                                  >  2:580691‑M1/1‑142 (MQ=255)                                                                                                                                                                                      cATTACAATTGGAAATAAAATTGTTTCCAATAGACATTGTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacgaaccgt                                                              <  2:641596‑M1/250‑109 (MQ=255)                                                                                                                                                                                          aCAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttaccc                                                         <  2:772935‑M1/251‑114 (MQ=255)                                                                                                                                                                                                      aaaaTTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggcttaatactctttacctgttacacggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacgg                                             >  2:457238‑M1/1‑126 (MQ=255)                                                                                                                                                                                                         aTTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtcc                                          <  1:634235‑M1/251‑129 (MQ=255)                                                                                                                                                                                                               tCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttatccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcac                                     <  1:792224‑M1/250‑134 (MQ=255)                                                                                                                                                                                                               tCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacg                                    >  1:540894‑M1/1‑117 (MQ=255)                                                                                                                                                                                                                  aaTAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgt                                 >  2:508818‑M1/1‑114 (MQ=255)                                                                                                                                                                                                                   aTAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgtt                                >  2:384604‑M1/1‑113 (MQ=255)                                                                                                                                                                                                                    tAGACATTTTTAACATGTTGCTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttc                               >  2:805443‑M1/1‑112 (MQ=255)                                                                                                                                                                                                                      gACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgc                             <  1:497700‑M1/251‑142 (MQ=255)                                                                                                                                                                                                                        cATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgccc                           >  1:589522‑M1/1‑108 (MQ=255)                                                                                                                                                                                                                          tttttAACCTGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgccc                           >  1:92614‑M1/1‑106 (MQ=255)                                                                                                                                                                                                                            tttAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggaccatcacgcgttcgcccagc                        >  2:392924‑M1/1‑104 (MQ=255)                                                                                                                                                                                                                                  aTGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttgg                 <  2:218871‑M1/251‑154 (MQ=255)                                                                                                                                                                                                                                    gttgttTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtc               <  1:771682‑M1/251‑156 (MQ=255)                                                                                                                                                                                                                                     ttgttTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtc               >  1:483612‑M1/1‑95 (MQ=255)                                                                                                                                                                                                                                       gttTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatccaagttgaagatcacccacaaactgttacccagtcgcatacggtccatcacgcgttcgcccagcaacttggtcatc            >  2:552663‑M1/1‑93 (MQ=255)                                                                                                                                                                                                                                          tttCTACGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgccccgcaacttggtcatctct         >  2:710147‑M1/1‑90 (MQ=255)                                                                                                                                                                                                                                          tttCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctct         >  1:312771‑M1/1‑90 (MQ=255)                                                                                                                                                                                                                                                 gtgtTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttccatat  >  2:697279‑M1/1‑83 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      GCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCC  >  NC_000913/4518111‑4518601

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