Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,996,049	T→A	L156* (TTG→TAG)	yigA →	DUF484 domain‑containing protein YigA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,996,049	0	T	A	100.0%	136.3 / NA	38	L156* (TTG→TAG)	yigA	DUF484 domain‑containing protein YigA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (18/20);  total (18/20)

AAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTTGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGCGT  >  NC_000913/3995830‑3996286                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              aaGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCtt                                                                                                                                                                                                                <  2:6834/251‑1 (MQ=255)        gTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGAc                                                                                                                                                                                                         <  2:270908/251‑1 (MQ=255)           ttATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGAcc                                                                                                                                                                                                        >  1:12075/1‑249 (MQ=255)               cGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCAGTAAGCCGGCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAgc                                                                                                                                                                                                   <  2:152756/250‑1 (MQ=255)                 ccTACTCTACCTGCAGCGCAGTCTCACCGCCGCCTGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAgctgc                                                                                                                                                                                                >  1:155165/1‑251 (MQ=255)                    aCTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCGGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCATTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGATCAGCACTATCTTGGGCCGCTTAACGGACCAGAgctgct                                                                                                                                                                                               <  2:64218/249‑1 (MQ=255)                         aCCTGCAGCGGAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCGGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTAGTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCtggtgg                                                                                                                                                                                          <  2:123849/249‑1 (MQ=255)                           cTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCt                                                                                                                                                                                      >  1:177599/1‑251 (MQ=255)                                        tcACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGaa                                                                                                                                                                          >  1:131693/1‑250 (MQ=255)                                          aCCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGc                                                                                                                                                                       <  1:80390/251‑1 (MQ=255)                                                             acgaTATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATg                                                                                                                                                    <  2:12075/251‑1 (MQ=255)                                                                   tGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGGGGCGGTGTCGATg                                                                                                                                              >  1:170475/1‑251 (MQ=255)                                                                               ttCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCg                                                                                                                                    <  1:195443/249‑1 (MQ=255)                                                                                   ccGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTg                                                                                                                              <  1:229316/251‑1 (MQ=255)                                                                                     gCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGAt                                                                                                                            <  2:51472/251‑1 (MQ=255)                                                                                       tGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGAttt                                                                                                                          >  1:92591/1‑251 (MQ=255)                                                                                        gggCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGAttt                                                                                                                          <  1:141054/250‑1 (MQ=255)                                                                                         ggCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGGGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTgg                                                                                                                        >  2:141048/1‑251 (MQ=255)                                                                                                 aTCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGGCCAGAGCTGCTGGGGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGGGATGTCGATGCTGGGAAGGGATGCTGGTTTgggggtgggg                                                                                                                >  2:62322/1‑244 (MQ=255)                                                                                                     cGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGTGTGTCGTGc                                                                                                               <  2:92591/248‑1 (MQ=255)                                                                                                        ccTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGttt                                                                                                          >  1:310753/1‑250 (MQ=255)                                                                                                           ggCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCa                                                                                                      <  1:162087/251‑1 (MQ=255)                                                                                                                       tCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGccc                                                                                          >  2:88740/1‑250 (MQ=255)                                                                                                                              cTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACt                                                                                    >  1:8593/1‑250 (MQ=255)                                                                                                                                 tttCCGGATCGATGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGAACAGAGATGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCCGTCACTATc                                                                                 >  2:16764/1‑250 (MQ=255)                                                                                                                                     cGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATcaacaa                                                                            <  1:309143/251‑1 (MQ=255)                                                                                                                                                cgcTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAAc                                                                 >  2:91582/1‑251 (MQ=255)                                                                                                                                                 gcTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACg                                                                <  1:29064/251‑1 (MQ=255)                                                                                                                                                              gTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATg                                                   <  1:195896/251‑1 (MQ=255)                                                                                                                                                                cGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGaa                                                 <  1:203109/251‑1 (MQ=255)                                                                                                                                                                         cTCATCTGGCATTCAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGACCTATCTTGGGCAGCTTAACGGACCAGAGCTGGTGGGGGTGCGACCGGAAGCGAAAGAGGTGGGCTCGGCGGCGATGCAGATGCTGGGAAGCGATGCTGATTTGGGTGTGGTGCTGTTAACCAGTAGCGATGCCCGTCCATATGAACAAGGGTAAGGAAGGGAGTTACTTCCTGAACTTgagtgg                                        >  2:175274/1‑246 (MQ=255)                                                                                                                                                                            aTCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGGGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATg                                     >  2:199535/1‑251 (MQ=255)                                                                                                                                                                                  cATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCg                               <  1:78888/251‑1 (MQ=255)                                                                                                                                                                                             cAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGCACAGCACTATCTTGGGCCGATTAACGGCCCAGAGATGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGCGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTgg                     >  1:23836/1‑250 (MQ=255)                                                                                                                                                                                                         ccGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTg         >  1:114450/1‑250 (MQ=255)                                                                                                                                                                                                           gCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAAc      <  2:78866/251‑1 (MQ=255)                                                                                                                                                                                                           gCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAAc      <  2:197735/251‑1 (MQ=255)                                                                                                                                                                                                               cGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGCGt  >  1:304707/1‑251 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              AAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTTGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGTTACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGCGT  >  NC_000913/3995830‑3996286

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