Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,012,448	T→C	K123K (AAA→AAG)	metR ←	DNA‑binding transcriptional dual regulator MetR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,012,448	0	T	C	100.0%	210.4 / NA	59	K123K (AAA→AAG)	metR	DNA‑binding transcriptional dual regulator MetR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/31);  total (28/31)

TGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGATTTAAAATCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAAC  >  NC_000913/4012156‑4012709                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                      tgtTATCGACGCTTTTCAGTGACGTGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTaaatc                                                                                                                                                                                                                                                                <  2:644431/300‑3 (MQ=255) tgtTATCGACACTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTaaatc                                                                                                                                                                                                                                                                <  2:104093/300‑3 (MQ=255) tgtTAGCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTaaatc                                                                                                                                                                                                                                                                <  1:312785/300‑3 (MQ=255)           gCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTAcc                                                                                                                                                                                                                                                     <  1:558891/300‑1 (MQ=255)           gCTTTTCAGGGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTAcc                                                                                                                                                                                                                                                     <  2:58003/300‑1 (MQ=255)            cTTTTCAGTGACGGGCTGCCGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCCCTGCGCGGCAGAATATCGGACGTCATTACCAGATTCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTATAcc                                                                                                                                                                                                                                                     >  1:13587/1‑299 (MQ=255)              tttCAGTGACTGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTCAA‑TCCATCTCTACCTGc                                                                                                                                                                                                                                                  >  1:558423/1‑300 (MQ=255)              tttCAGTGACGGGCTGACGCCTGCAGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTAGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAANATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGATGCGGGTCAAATGTCACGCCCGACTTAAA‑TCAATCTATACCt                                                                                                                                                                                                                                                    >  2:510862/1‑298 (MQ=255)                        gggCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTc                                                                                                                                                                                                                                         <  1:601422/299‑1 (MQ=255)                                        gCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTaa                                                                                                                                                                                                                        >  1:80418/1‑300 (MQ=255)                                        gCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTAGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCAACTTAAA‑TCCATCTATACCTGCGGCCAGTTCTTATGGAAGTTTTCTa                                                                                                                                                                                                                         >  1:80857/1‑299 (MQ=255)                                                          aGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGc                                                                                                                                                                                                           <  2:401994/295‑1 (MQ=255)                                                               tCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCNCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATAc                                                                                                                                                                                                 <  1:135228/300‑1 (MQ=255)                                                                  attCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGc                                                                                                                                                                                              <  1:570014/298‑1 (MQ=255)                                                                  aGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGc                                                                                                                                                                                              <  1:569906/300‑1 (MQ=255)                                                                       aCTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCNAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATgg                                                                                                                                                                                         <  2:504417/300‑1 (MQ=255)                                                                       aCTACGCTGCACCGGAGAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATgg                                                                                                                                                                                         <  2:623711/300‑1 (MQ=255)                                                                               gCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAACATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTAACGNCCGACTTAAA‑TCCATCTATACCTGAGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATggagctc                                                                                                                                                                                    >  2:494677/1‑292 (MQ=255)                                                                               gCACCGGATAAATTAATAGCGTCTAGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCAAGACCGACTTAAA‑TCCATCTATACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGATATGGCACTCAAt                                                                                                                                                                                 >  2:167152/1‑300 (MQ=255)                                                                                    ggATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCaa                                                                                                                                                                            <  2:296249/300‑1 (MQ=255)                                                                                    ggATAAAGTAAGAGCGTCTCGCTGGCGAGATCGGCAGGTGGAATTCGCGTGTGCGCAGCAAGTGGATGGTCAGGTGCGAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCTGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCaa                                                                                                                                                                            <  2:296242/300‑1 (MQ=255)                                                                                         aaTTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACNTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATgcgc                                                                                                                                                                       >  1:124092/1‑300 (MQ=255)                                                                                                aGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGAcgcg                                                                                                                                                                <  2:413429/300‑1 (MQ=255)                                                                                                      tcGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTctgc                                                                                                                                                           <  2:529031/299‑1 (MQ=255)                                                                                                                        ggcgtAATTCTCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACGTCAGAGTCGAACATCGGCGAATAATGCAGGCCACGGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCCGGTCAAATTTCACGCCCGACTTAAA‑TCCATCTCTACCTg                                                                                                                                                                                                                                                   <  2:185700/190‑1 (MQ=255)                                                                                                                                                   ggATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCAAGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTa                                                                                                             >  2:19182/1‑300 (MQ=255)                                                                                                                                                         tCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTAACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGt                                                                                                       <  1:42671/300‑1 (MQ=255)                                                                                                                                                         tCAGGGGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCATGCCACTGCGCTTCAGAATATCGGACTTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACT‑GAAATCCAGCTCTGTCGGTGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATGCAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCCTGCAGGGCCCGGCTAATTTGCGGCAGTACCTGGt                                                                                                       <  2:163824/300‑1 (MQ=255)                                                                                                                                                              tGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGcc                                                                                                  >  1:375240/1‑300 (MQ=255)                                                                                                                                                               gCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCa                                                                                                 >  1:203066/1‑300 (MQ=255)                                                                                                                                                               gCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTANGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCa                                                                                                 >  2:497330/1‑300 (MQ=255)                                                                                                                                                                             cACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGAtt                                                                                   >  1:635399/1‑300 (MQ=255)                                                                                                                                                                               cTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTtc                                                                                 >  1:103429/1‑300 (MQ=255)                                                                                                                                                                                           aCATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTTTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCATCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGgtgt                                                                     <  2:412174/300‑1 (MQ=255)                                                                                                                                                                                           aCATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGgtgt                                                                     <  2:412402/300‑1 (MQ=255)                                                                                                                                                                                                tgCGAATAATGCAGGCTACTGCGCGGCAGAAGATCGGGCATCATTACCAGAGCCAACCCTCAATGTGGCAAGGCGGGCTGCGGGGCAAATGTCACGCCCGACTTAGA‑TCCATCCCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTTTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACACGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCGACAGGATTTCTCCCTGCGGTGTAAAGc                                                                <  2:133348/299‑1 (MQ=255)                                                                                                                                                                                                ggCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGc                                                                <  1:432667/300‑1 (MQ=255)                                                                                                                                                                                                ggCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCNGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATTGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGc                                                                <  1:453977/300‑1 (MQ=255)                                                                                                                                                                                                ggCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCNGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGc                                                                <  1:453764/300‑1 (MQ=255)                                                                                                                                                                                                   gAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCGTTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGt                                                              <  2:615277/299‑1 (MQ=255)                                                                                                                                                                                                   gAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGNCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTa                                                             >  1:129326/1‑300 (MQ=255)                                                                                                                                                                                                   gAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTa                                                             >  1:646577/1‑300 (MQ=255)                                                                                                                                                                                                             ggCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCggctggct                                                   >  1:111366/1‑300 (MQ=255)                                                                                                                                                                                                               ccACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCNGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGctct                                                 >  1:454211/1‑300 (MQ=255)                                                                                                                                                                                                               ccACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTAACTGCGGCAAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGctct                                                 >  2:620669/1‑300 (MQ=255)                                                                                                                                                                                                                      gcgGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGAGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGAAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCAGTGTAAAGCGTAGCGGCTGGCTCCTACGcac                                          >  2:170643/1‑300 (MQ=255)                                                                                                                                                                                                                         gCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGTTCTTACGCACaaa                                       >  2:437326/1‑300 (MQ=255)                                                                                                                                                                                                                            gAATATCGGACGTCATTACCAGATCCAGCTCTCACTGTTGCAAGGCGGGCTGCGGGTCAAATGTCAAGCCCGACTTAAA‑TCCATCTATAACTGCGGCCAGTTCTTATGGAAGTTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCATGGCTAATTTGCGGCAGTACCTGGTTTGCAGGTTGCAACAGGATTTCTAACTGCGGTGTAAAGCGTAGCGGCTGGCCCTCACGCACAAATa                                     >  2:72508/1‑299 (MQ=255)                                                                                                                                                                                                                                                 gATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTc               <  2:587178/300‑1 (MQ=255)                                                                                                                                                                                                                                                 gATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCGATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTc               <  2:587467/300‑1 (MQ=255)                                                                                                                                                                                                                                                  aTCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCAAAGGCGTTGTTc               >  2:563502/1‑299 (MQ=255)                                                                                                                                                                                                                                                  aTCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCAAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGNAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCgttgtttc              >  1:141664/1‑298 (MQ=255)                                                                                                                                                                                                                                                      aGCTCTCCCTGTTGCAAGGCGGGCTGCTGGTCAAATTTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGAt          <  1:51150/300‑1 (MQ=255)                                                                                                                                                                                                                                                          ctcCCTGTTGCAAGGCGAGCTGCGGGTCAAATATCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATCACCGGAAGCCAAGGCGTTGTTCCAGATCgat      >  2:173581/1‑298 (MQ=255)                                                                                                                                                                                                                                                            cccTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCACTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTa     >  1:181157/1‑299 (MQ=255)                                                                                                                                                                                                                                                              cTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAAc  >  1:86491/1‑300 (MQ=255)                                                                                                                                                                                                                                                              cTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGCGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAACGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATc         >  2:343277/1‑293 (MQ=255)                                                                                                                                                                                                                                                                gTTGCAATGCGGGCTGCGGTTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATTGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAAc  <  2:103429/298‑1 (MQ=255)                                                                                                                                                                                                                                                                 ttGCAATGCGGGCTGCGGGTCAAATGTCACGCCCGACTTAAA‑TCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCNTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAAc  <  2:468783/297‑1 (MQ=255)                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                      TGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGATTTAAAATCCATCTCTACCTGCGGCCAGTTCTTATGGAAATTTTCTAACGCGGGTGTCAGCCACTGAATACAGCTATGGCACTCAATGGCAATGCGCAGACGCGTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAAC  >  NC_000913/4012156‑4012709

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