Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,754	T→C	Y41H (TAT→CAT)	mraZ →	DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,754	0	T	C	100.0%	234.6 / NA	73	Y41H (TAT→CAT)	mraZ	DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/33);  total (40/33)

GTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTATCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAGTTTGAGCTGTGGGATG  >  NC_000913/89511‑89997                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    gTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGGGGGGCAAAGTGGGCATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGAAAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCACCAccc                                                                                                                                                                                                                                              >  2:259898/1‑251 (MQ=255)        tCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCctg                                                                                                                                                                                                                                       >  1:63195/1‑251 (MQ=255)                tttACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCttt                                                                                                                                                                                                                                <  1:362749/250‑1 (MQ=255)                 ttACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAAGTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTa                                                                                                                                                                                                                               <  1:665496/250‑1 (MQ=255)                  tACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTAc                                                                                                                                                                                                                              <  1:169423/250‑1 (MQ=255)                   aCCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTAccc                                                                                                                                                                                                                            >  1:258255/1‑251 (MQ=255)                    ccTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCCCCGTAAACTCCTTTCAGTGGGGAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTAcccc                                                                                                                                                                                                                           >  2:173483/1‑251 (MQ=255)                        ggTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCcctgc                                                                                                                                                                                                                       <  1:519408/251‑1 (MQ=255)                         gTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCCCCGTAAACTCCTTTCAGTGGGAGATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCCAAGGGCGCTTATCAATGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACCTTCATCACCCGGGCCTCATGCTTTACCCcccgcc                                                                                                                                                                                                                      >  2:510429/1‑251 (MQ=255)                          gaTGCCTTGTCGCTGGCCTGAAAAGCTTTTTCTACGCGCCCTAAACTCCTTTCAGTGGGAAATTTTGGGGCAAAGGGGGGATAAGGGGTGAAGCTGGCATGTTCCGGGGAGCAACGTTAGGCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCTGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATTGTCTGCACCATTGACATTCATCACCCGTTCCTGCTGTTTTACCCcctgcct                                                                                                                                                                                                                     <  2:336939/250‑1 (MQ=255)                                  gtgACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATggg                                                                                                                                                                                                              <  2:63195/250‑1 (MQ=255)                                         gCTTGACAAGCTTGTCCCCCGCGCCTGAAGCTCGTTTCAGTGGGAAAGTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCGCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGGCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATc                                                                                                                                                                                                      <  2:106116/251‑1 (MQ=255)                                          cTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCg                                                                                                                                                                                                     >  1:44104/1‑251 (MQ=255)                                             gACAAGCTTTTCCTCCTCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAACGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCCAATGGTTTGCACCATTGCAATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGATTGGGAAATTATCGAGc                                                                                                                                                                                                  >  2:219076/1‑251 (MQ=255)                                             gACAAGCTTTTCCTCCGCTCCGTAAACTCCTTGCAGTGGGAAATTGGGGGGCAACGTGGAAATAAGGGGTGAGGCTTGCATGTTCCTGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTTCCTACCCTTTCCCGGGAACAGCTGCTTGAGAACTCTGCCGTTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTTCTTTACCCCCTGCCTGAATGGGAATTTATCGAGc                                                                                                                                                                                                  <  2:282088/251‑1 (MQ=255)                                                    ttttCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGATGCTTTACCCCCTGCGTGAATGGGAAATTATCGGGCAAAAAtt                                                                                                                                                                                           >  2:474455/1‑251 (MQ=255)                                                            cgatCCGTAAACGTCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGATTGCATGTTCCGGGGAGCAACGTTAGGCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCGGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTc                                                                                                                                                                                   <  2:258255/248‑1 (MQ=255)                                                                   tAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCttt                                                                                                                                                                                                                                >  1:147374/1‑199 (MQ=255)                                                                   tAAACTCCTTTCAGTGGGAAATTGGGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCTGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCTCCATTGACATTCATCACCCGTGCCTGCTGCttt                                                                                                                                                                                                                                <  2:147374/199‑1 (MQ=255)                                                                      aCTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCAt                                                                                                                                                                         <  1:359670/251‑1 (MQ=255)                                                                       cTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGCGCAAAAATTATCGCGTCTGTCGAGCATg                                                                                                                                                                        >  1:599684/1‑251 (MQ=255)                                                                                         tgtgGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgt                                                                                                                                                      >  1:426032/1‑251 (MQ=255)                                                                                            ggggCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGtgtg                                                                                                                                                   >  2:195823/1‑251 (MQ=255)                                                                                                           tAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTAAATATCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATAGAGAAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGttt                                                                                                                                    >  2:320711/1‑250 (MQ=255)                                                                                                            aaGGGGTGAGGCTGGCATGTTCCGGGGAGCAAAGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAACTGGTTTGCACCATTGCCATTCATCACCCGCGCCTGCTGCTTTACCCCCTGCATGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCTTGAACCCGGATAAGCGCCGTGTGCAGCGACTACTGttgg                                                                                                                                   >  2:473746/1‑249 (MQ=255)                                                                                                            aaGGGGTGAGGCTGGCATGTCACGGGGAGCAAAGTTAGTCACTCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTACTTGAGAACGCTGCCGGTAAAATGGTTTGCACCATTGCAATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTAACGAGCAAAGATTATCGCGTCTGTACAACATGAACCCGGTTGCGCGCCGTATGCAGCGACTACTGta                                                                                                                                     >  2:666400/1‑248 (MQ=255)                                                                                                                ggTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATATCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTc                                                                                                                                >  1:619015/1‑250 (MQ=255)                                                                                                                  tGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATg                                                                                                                             >  1:641612/1‑251 (MQ=255)                                                                                                                   gAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGc                                                                                                                            <  1:586146/251‑1 (MQ=255)                                                                                                                       caggCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAACTGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAACTTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGACCCCGGTTGCGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGc                                                                                                                        >  2:250847/3‑251 (MQ=255)                                                                                                                              gTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGt                                                                                                                  <  1:259898/250‑1 (MQ=255)                                                                                                                                   ggggAACAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCATACTGTTAGGTCATGCCAGCGAATGTCAgatg                                                                                                            >  1:154941/1‑251 (MQ=255)                                                                                                                                    gggAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatgg                                                                                                           <  2:286949/251‑1 (MQ=255)                                                                                                                                    gggAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatgg                                                                                                           >  1:232885/1‑251 (MQ=255)                                                                                                                                         cAACGTTAGTCAATCTCGACAGCAAAGGGCTCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGgc                                                                                                      <  1:636127/251‑1 (MQ=255)                                                                                                                                          aaCGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGgcg                                                                                                     >  1:124696/1‑251 (MQ=255)                                                                                                                                                    aaTCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCcctg                                                                                                                                                                                                                        >  2:131316/1‑126 (MQ=255)                                                                                                                                                    aaTCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCcctg                                                                                                                                                                                                                        <  1:131316/126‑1 (MQ=255)                                                                                                                                                         cGACAGCAAAGGGTGCTTATCAGTGCCTTTCTGTTCCCGGGACCACCTGTTTTAGAGCGCTGCTGGTCAAAGGGTTTGCATCCCTGATATTCATCCTCAGTGCCTGCTGCTTTACCCCCTGCCTGACTGGGAAATTATCGAGCATAAATCATCGCGTCTGTCGAGCATGTACCCGGTTGAGCGCCGTGTGCAGTGCCTCCTGTTTGGTCAGGCTAGCGAACTTCAGATGGATGGCGCAGGTCGATTGTTaa                                                                                      <  1:517289/251‑1 (MQ=255)                                                                                                                                                            cAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCCCCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTATGCAGCGCCTACTGTTAGGTCATGCCCGCGAATGTCAGATGGATGGCGCCGGTCGATTGTTaaacg                                                                                   >  2:80632/1‑248 (MQ=255)                                                                                                                                                             aGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATcgc                                                                                  <  1:219076/251‑1 (MQ=255)                                                                                                                                                               cAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCgcgc                                                                                >  1:658178/1‑251 (MQ=255)                                                                                                                                                                 aaGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGcc                                                                               >  1:381773/1‑250 (MQ=255)                                                                                                                                                                   gggCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGt                                                                            >  1:197328/1‑251 (MQ=255)                                                                                                                                                                   gggCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGt                                                                            <  1:350930/251‑1 (MQ=255)                                                                                                                                                                    ggCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatggat                                                                                                         >  1:318081/1‑221 (MQ=255)                                                                                                                                                                    ggCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatggat                                                                                                         <  2:318081/221‑1 (MQ=255)                                                                                                                                                                              agggCCTACCCGTCACCGGGAACAGTTGCTTTAGACCTCTGCCGGTCAAATGGTGTGCTCCATTGACATGCATCACCCGTGCCTGCTGCTTTACCACCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGGCATGCCAGCGAATGTCAGATGGATGGCGCAGGTAGATTGTTAATCTCGCCAGTACTGCTGCAAc                                                                 <  2:160753/248‑1 (MQ=255)                                                                                                                                                                                   cTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGTATGTc                                                                                                                 >  2:364915/1‑198 (MQ=255)                                                                                                                                                                                   cTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTc                                                                                                                 <  1:364915/198‑1 (MQ=255)                                                                                                                                                                                   cTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGcc                                                            >  1:533071/1‑251 (MQ=255)                                                                                                                                                                                     aCCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTATTGCGGCAACATGCCg                                                           <  1:378222/250‑1 (MQ=255)                                                                                                                                                                                             ccGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACa                                                   >  1:96164/2‑250 (MQ=255)                                                                                                                                                                                               gggAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACaaa                                                 >  1:655108/1‑250 (MQ=255)                                                                                                                                                                                                    cAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCAAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGt                                           <  1:177797/251‑1 (MQ=255)                                                                                                                                                                                                      gctgctTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGTCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGa                                         <  1:409041/251‑1 (MQ=255)                                                                                                                                                                                                         gctTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATg                                       <  1:674625/250‑1 (MQ=255)                                                                                                                                                                                                                    cTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGACGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTATTCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGGGATGCTGGTTggcc                             >  1:123834/1‑247 (MQ=255)                                                                                                                                                                                                                     tGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTCTACCCCCTGCCTGAATGGGAAATTATCGCGCAAAAATTATCGCGTCTGTCGAGCATGAAGCCGGTTGAGCGCCGTGTGCCGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTCATCGCGCCAGTACTGCGGCAACATGCCGGTCTGCAAAAAGAAGTGATGCTGGTTGGACAGt                          >  2:365191/1‑251 (MQ=255)                                                                                                                                                                                                                      gCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAAGGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGtt                         <  2:44104/251‑1 (MQ=255)                                                                                                                                                                                                                       ccGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTc                        >  1:578879/1‑251 (MQ=255)                                                                                                                                                                                                                          gtaaaaTGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCGTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTAGCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGAGTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTcaac                     <  2:589785/248‑1 (MQ=255)                                                                                                                                                                                                                          gTCAAATGGTTTGCACCATTGACATTCATCACCAGTGCCTGCTGCTTTACACCCTGCCTGAATGGGAAATTATCGAGCACAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGCTGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTcaac                     >  2:445109/1‑251 (MQ=255)                                                                                                                                                                                                                             aaaTGGTTTGCTCCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATTGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGt                          <  1:103737/243‑1 (MQ=255)                                                                                                                                                                                                                             aaaTGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTcaaca                    >  1:59355/1‑249 (MQ=255)                                                                                                                                                                                                                             aaaTGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAg                  <  1:320711/251‑1 (MQ=255)                                                                                                                                                                                                                             aaaTGGATTACACCCTTGACATTCATCACCCGTGCCTGATGCTTTACACCCTGCCTGAATGGGAAATTATAGAGCAAAAATTATCGCGTCTGTCGAGCATGCACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAga                          >  2:103737/1‑242 (MQ=255)                                                                                                                                                                                                                                       cACCATTGACATTCATCTCCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCCAAAATTATCGAGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGAGCAGCGCCTACTGTTACGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCCACAGTACTGCGGCCACACGCCCGGCTGACAACAGAAGTGATGCTGGTTGGCCAGTTCAACAAGGTTGAGCtgt        >  1:647218/1‑251 (MQ=255)                                                                                                                                                                                                                                          cATTGACATTCATCACCCTTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAGTTTGAGCTGTggg     <  1:436400/251‑1 (MQ=255)                                                                                                                                                                                                                                              gACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTCTCGAGCAAAAATTATCGCGTCTGTAGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGAGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGAAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGGCAGTTCAACAAGTATGAGCTGGGGGATg  >  2:434037/1‑250 (MQ=255)                                                                                                                                                                                                                                              gACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAGTTTGAGCtgtg       >  2:552652/1‑245 (MQ=255)                                                                                                                                                                                                                                              gACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAGTTTGAGCtgtg       <  1:552652/245‑1 (MQ=255)                                                                                                                                                                                                                                              gACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAGTTTGAGCTGTGGGATg  <  1:445109/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    GTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTATCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAGTTTGAGCTGTGGGATG  >  NC_000913/89511‑89997

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