Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,587	C→T	T137I (ACC→ATC)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,587	0	C	T	100.0%	137.6 / NA	45	T137I (ACC→ATC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/23);  total (22/23)

CTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGA  >  NC_000913/4638340‑4638833                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        cTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCtt                                                                                                                                                                                                                                                      >  1:177790/1‑251 (MQ=255)             cttaTCGCTATTCCTGTGACCGAACTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGCAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGccctt                                                                                                                                                                                                                                          >  2:227906/4‑248 (MQ=255)               taTCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCttt                                                                                                                                                                                                                                         <  2:161447/249‑1 (MQ=255)                  cGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTAtt                                                                                                                                                                                                                                      <  2:177790/250‑1 (MQ=255)                  cGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCACTTTAtt                                                                                                                                                                                                                                      <  2:488695/250‑1 (MQ=255)                   gCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGCGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATtgt                                                                                                                                                                                                                                    >  2:370931/1‑251 (MQ=255)                        tCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGAtt                                                                                                                                                                                                                                <  2:18549/250‑1 (MQ=255)                                      tttATACGGTGCAGGACGAGGATAAAACCGTGTAGCGGAACCGCTGTTTTATCCAATTTTGGTTAACTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGACATTTATTGTGATTCGCGTCGGGGATgcg                                                                                                                                                                                                                 >  2:213149/1‑251 (MQ=255)                                          taCGGTGCAGGCAGAGGATAAAACCGTGGAGAGGAAACGAAGTTTTCTCCATTTTTGGTTACCTGAGTCATTGCTGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGATAAGGACTGGCAAAGTGATTTAACGCTCAAAGCCGATTTCACTGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCCAGCCCTTTCTTGAGATTAGCGTCGGGAATGCGCGTg                                                                                                                                                                                                             >  2:572509/1‑251 (MQ=255)                                                       gagtaTAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑tggtgtggt                                                                                                                                                                                               <  1:126207/247‑1 (MQ=255)                                                                   gTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGa                                                                                                                                                                                   <  1:263800/251‑1 (MQ=255)                                                                      gAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCTATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCTTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGt                                                                                                                                                                                <  2:496634/251‑1 (MQ=255)                                                                              aCGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAgc                                                                                                                                                                                         <  2:506901/234‑1 (MQ=255)                                                                              aCGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAgc                                                                                                                                                                                         >  1:506901/1‑234 (MQ=255)                                                                                gAAGTTTTATCCATTTTTCGTTACCTGAGCCATTGATGCTTGATGACAATCAGCACGTGGACGAACGCAAGCTAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCT‑TGGTGTGGGGAAAGCGCCTGAAGTTaacggaccg                                                                                                                                                                      >  2:93339/1‑248 (MQ=255)                                                                                       tATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTg                                                                                                                                                                <  1:468220/250‑1 (MQ=255)                                                                                         tCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGAcc                                                                                                                                                             >  1:130359/1‑251 (MQ=255)                                                                                            aTTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAGTCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCAtt                                                                                                                                                          <  2:113067/251‑1 (MQ=255)                                                                                              ttttGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGAGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGAGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTg                                                                                                                                                         >  1:362060/1‑250 (MQ=255)                                                                                                 tGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAAcc                                                                                                                                                     >  1:337832/1‑251 (MQ=255)                                                                                                 tGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAAc                                                                                                                                                      <  1:107097/250‑1 (MQ=255)                                                                                                   gTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCg                                                                                                                                                   <  1:308262/251‑1 (MQ=255)                                                                                                      aCCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTATGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAt‑t                                                                                                                                                                                                       <  1:383953/196‑1 (MQ=255)                                                                                                      aCCTGAGTAATTGATGGTTGATGGCAATCAGCACGTGGCAGAACGCAAGATAGGGCTTTATACCGGTCAGGTCTGGCACAGAGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCACATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGGAt‑t                                                                                                                                                                                                       >  2:383953/1‑196 (MQ=255)                                                                                                                   atggttgatggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTTATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGGGGTGAAAGCGACTGAAGTTAACGGAACGGCGCt                                                                                                                                                                 <  1:203069/221‑1 (MQ=255)                                                                                                                                aTCAGAACATGGAAGAACGCAAGATAGGGCTTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAACGCCGCTTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGAAAGGCGGGCAGGGc                                                                                                                      >  1:229075/1‑251 (MQ=255)                                                                                                                                       cGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCata                                                                                                               <  2:448762/251‑1 (MQ=255)                                                                                                                                          ggAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTCGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGCTAACGGAACGGCGCGGGCACTTGAACCCGGCACCGGGTTAGAGCACGGCGGACCGGGAGTGCatat                                                                                                              >  2:598920/1‑249 (MQ=255)                                                                                                                                          ggAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTGTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGGCCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCatat                                                                                                              <  1:598920/249‑1 (MQ=255)                                                                                                                                             agaaCGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCtt                                                                                                         <  1:511508/251‑1 (MQ=255)                                                                                                                                                actcAAGATAGTGACTTGTACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATGTCGATGGTTCGCGTCTTAGCTAGCGCAACGCGCCAAAGATCATCTTAGGCAAGCCCTTTAGTGTGATTGTCGTCGGGGATTCGCGTTGTAC‑TGGTGTTGTGTAAGCGCCTGATGTGAACGGAACGTCGCTGACCATTGAACCCCGCACCTGGTTAGAGCAAGGCTGGCTTGTCGTGCATAGCCCTTTAc                                                                                                      <  2:50620/248‑1 (MQ=255)                                                                                                                                                       aTAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGGACCGGGTTAGAGCAAGTCGGGCAGGGCGTGCATATCCCGTTACCCGAAgg                                                                                               <  2:82772/251‑1 (MQ=255)                                                                                                                                                          gggATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGa                                                                                            <  1:302400/251‑1 (MQ=255)                                                                                                                                                                tataCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGGAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCGGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGc                                                                                       <  1:227906/250‑1 (MQ=255)                                                                                                                                                                      ggTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTAGATGTTTCACGTCTTAGCGAGCTAAACGCGCCAAATATCATCTTAGGCCAGCCCTTTATTGTGATTAGCGACGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCAGAAAGGGACTGGAGGaacc                                                                                 >  1:290543/1‑248 (MQ=255)                                                                                                                                                                                   acagtgagttcacgctatacgCCGATTTCGATGTTTCGCGTCTTACCGAGCTCAACGCCCCAAATATCATATTAGTCACGCCCTTTATTGTGATTAGCGTCGGGGCTGCGCGTGGTAT‑TGGCGTGGTGATAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAACAGACGCTGAAGCTg                                                                   >  2:95029/21‑251 (MQ=255)                                                                                                                                                                                             aaCGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGccc                                                         >  1:431291/1‑251 (MQ=255)                                                                                                                                                                                                     aaGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGCACGGCGCTGACCCTTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGAGTGCATATACCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTAAAGCTTAATATGGCCCTGAAtttt                                                 >  2:557990/1‑250 (MQ=255)                                                                                                                                                                                                      aGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGCACGGCGCTGACCATTGAACCCGGCACCGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAGGGGGACTGGCGGAAGCAGAACCTGGAGccggataatgccctgtatata                                                 >  2:110076/1‑230 (MQ=255)                                                                                                                                                                                                               cGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTAcc                                        <  2:345456/250‑1 (MQ=255)                                                                                                                                                                                                                gATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCACGCCCTTTATTGTGATTAGCGTCAGGGATGCGCGTGGCAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCGGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCCGGGCGTGCATATCCCTTTACCCGAAGGGGCCTGGCGGAAGCAGAAACTGAAGCTGGATATGGCCCTGAATTTAAGCGGTGCCgg                                      >  1:234160/1‑251 (MQ=255)                                                                                                                                                                                                                          ggctgaGCGAGCTCAACGCGGCAAATATCATCTTAGGCAAGCCCTTTAGTGTGATTAGCGTCGGGGATGCGCGTGTTATGTGG‑GTGGTGAAAGCGCCAGAAGTTAACGGAACGGCGCTGACCATTGAGCCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGGATATCCCGTGATCCGAAGGGGTCTGTCGGAGGCAGAGCCTGAAGCTCAATATGGGCCTGAATGTGAGCGGTCCGGGCGCTCTTTCt                            <  2:538537/246‑1 (MQ=255)                                                                                                                                                                                                                            cTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCCTTGAACCCGGCACCGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGGAGCAGAACCTGAAACTGAATAATGCCCTGAATTTAAAGGGTACACGCGATCCTTCtgg                          >  2:354458/1‑250 (MQ=255)                                                                                                                                                                                                                                          cgcgCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCCTATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGGATATGGCCCTGAATTTTAGCGGTACCGGCGATCTTTCTTTGGGGCCTgccgggg            >  2:155259/1‑250 (MQ=255)                                                                                                                                                                                                                                                  atatCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGc    >  1:236351/1‑251 (MQ=255)                                                                                                                                                                                                                                                    atcatcTTAGGCAAGCCCTTTATTGTGAGTAGCGTCGGAGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGTTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGa  <  2:170958/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        CTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGA  >  NC_000913/4638340‑4638833

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