Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,627,893	G→T	E193D (GAG→GAT)	nadR →	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,627,893	0	G	T	100.0%	177.1 / NA	58	E193D (GAG→GAT)	nadR	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (32/26);  total (32/26)

CAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGA  >  NC_000913/4627660‑4628138                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      cagcagCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCGAATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATcc                                                                                                                                                                                                                                      >  2:362649/1‑251 (MQ=255)                  ggATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGa                                                                                                                                                                                                                     <  1:188639/251‑1 (MQ=255)                     tcgtctgctctcgttagtacaaccTTTTAAATTTCAGACAAATATTCGCTATCATGCTTTAAACGAAGAGGGAATGGAGCATTATTCGGACGGCTGGGATGGGTGGAGCAACTGCATCAAAAAGTTTATGTCTGTAAAGGGTATTGAGCCGGACGTT‑ATCTACACCTCGGAAGAAGCCGATGCGCCCCAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGt                                                                                                                                                                                                                   <  2:45637/227‑1 (MQ=255)                           gcgctGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCg                                                                                                                                                                                                            >  1:512062/5‑251 (MQ=255)                                ggTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTgcg                                                                                                                                                                                                       <  1:480301/251‑1 (MQ=255)                                   tATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTgcgc                                                                                                                                                                                                      >  1:548739/1‑249 (MQ=255)                                    attacAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAg                                                                                                                                                                                                    <  1:572529/246‑1 (MQ=255)                                    attacAAACTTTTAAAGATCAGAAAAATAGTCGCTTTCATGATTTCAACGAAGAGGGCATGGAGCCGTTTCCGGTCTGCGGGGATGTGGGGAGCAACGGCAGCAAAAAGTTTATGGTTGAAAAGGATATGCAGCCCGCCCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCGGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAg                                                                                                                                                                                                    <  2:672892/246‑1 (MQ=255)                                        cAAACTTTTAAATATCAGAAAAATACTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGAGGTGCGCAGATCcg                                                                                                                                                                                               >  2:667311/1‑251 (MQ=255)                                        cAAACTTTTAAAAATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGCTTATGGATGAAAAGGGGATTCAGACGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTAGCTCCGAAACGTACCTTTATGAGTATAAGCGGTGCGCAGATccc                                                                                                                                                                                               >  2:312347/1‑250 (MQ=255)                                                 aaaTATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGAACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGGCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCCGATCCGCGAAAAccc                                                                                                                                                                                      >  1:76345/1‑251 (MQ=255)                                                   atatCAGAAAAATATTCGCATTCAGGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGACTGGGATGTGTGGCGCAACGGCATCAAAAAGTCTATGGCTGACAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGGGACACAGTATATGGAACATATGGGGATCGCTACGGTGCTGGTCGAACGTAAACGTACCTTTATGAGTCTAAGGGGTGCGCAGATCCGCGCAAACCCGt                                                                                                                                                                                    >  2:527331/1‑251 (MQ=255)                                                       cAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCg                                                                                                                                                                                <  1:295456/251‑1 (MQ=255)                                                          aaaaaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCtc                                                                                                                                                                             >  1:561247/1‑250 (MQ=255)                                                           aaaaTAATCGCAGTCATGCTTTCAACGAAGAGGGCATGGATCCGTAGCCGCACGGCTGGGATGTGTGGAGCAACGGCACCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTAc                                                                                                                                                                            <  1:618578/251‑1 (MQ=255)                                                             aaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTg                                                                                                                                                                          >  1:475808/1‑251 (MQ=255)                                                              atatTCGCATTCATGCTTTCAACGAAGAGGGGATGGAGCCGTATCCGCACGGCTGGGAGGTGTGGAGCAACGGCATCAGAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACTCCTCGGAAGAATCCGATGCGCCACAGTATATTGAACATCTGGGGATCGTTACGGTGCTGGTCGATCCGAAACGTTCCTTTATGAGTAGCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTgg                                                                                                                                                                         <  1:384599/251‑1 (MQ=255)                                                               tatTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTcc                                                                                                                                                                                 <  1:100543/242‑1 (MQ=255)                                                               tatTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTggg                                                                                                                                                                        >  1:175525/1‑251 (MQ=255)                                                               tatTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAAAATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCCCAGTTCCGCGAAAACCCGTTcc                                                                                                                                                                                 >  2:100543/1‑242 (MQ=255)                                                                    gCATTCATGCTTACAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCTTTCCGCTACTGGGAatat                                                                                                                                                                   <  1:426006/251‑1 (MQ=255)                                                                           tGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCAGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGAACCTTTATGATTATCAGCCGTGCACAGATCCGCGAAAACCCGGTcctattatggggatataatcct                                                                                                                                                             >  2:355115/1‑230 (MQ=255)                                                                                tCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTAcgga                                                                                                                                                        <  1:544492/250‑4 (MQ=255)                                                                                 caaagaggagGGCATGTAGCTGTTTCCGCACGGGTGGGATGTGGTGTGCATCGTCCTCAAAAATGTTATGGCTGAAAATGGGATTCAGACGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGAGCCGCGAAAACCCGTTCCGCGACTGGGAATGTATTccgacggaa                                                                                                                                                       <  2:92499/243‑8 (MQ=255)                                                                                                 gAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTtgt                                                                                                                                      >  1:205443/1‑251 (MQ=255)                                                                                                                    gggATGTCTGGATCAACGGCATCAAAAAGGTTATCACAGTAATGGGGATTCATCCGGA‑CTGGATCTACGCCTCGGAATAGTCCGATTCGCCACAGGATATTGAACATGTGGGTACCGATACGGTGCTGGTCTATCCGAAACGGATCTCTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGGTCCGCTACTGGGAATATTATCCTACGGAAGTGAAGCCGTGCCTTGTACGTACCGCGTCGATccgt                                                                                                                   <  2:99347/251‑3 (MQ=255)                                                                                                                           gtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTggtgtcg                                                                                                            >  1:111290/1‑246 (MQ=255)                                                                                                                                     ggCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCgg                                                                                                  <  2:205443/251‑1 (MQ=255)                                                                                                                                      gCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGt                                                                                                 >  1:140523/1‑251 (MQ=255)                                                                                                                                           aaaaaGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATc                                                                                            >  2:210179/1‑251 (MQ=255)                                                                                                                                                               ggggATTCAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTg                                                                        <  1:362649/250‑1 (MQ=255)                                                                                                                                                                    ttcAGCCGTACCTG‑ATCTACACCTCGGCAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGAAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATCTTCCTACGGAAGTGAAGCCGTTCTTTGTACGGACCGTg                                                                                                                            >  2:664094/3‑194 (MQ=255)                                                                                                                                                                    ttcAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTg                                                                                                                            <  1:664094/192‑1 (MQ=255)                                                                                                                                                                    ttcAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAat                                                                   >  1:371286/3‑251 (MQ=255)                                                                                                                                                                      cAGCCGGACCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGACGAGCGGTACATCCACCCTGGTAAACAAACTTGCCAatat                                                                 >  2:371793/1‑251 (MQ=255)                                                                                                                                                                       catccggacctgaactccactgtggAAGAAGCCGATGCGCCACAGTATATGGACCATCTTGGGATCGATACGGTGCTGGTCGATCCGAAACGTTCCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGGGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCTCCCTGGTAAACAAACTTGCTAata                                                                  <  2:142257/227‑1 (MQ=255)                                                                                                                                                                           ggacctg‑atctACACCGTGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCa                                                            <  1:383580/247‑1 (MQ=255)                                                                                                                                                                            gacctg‑atctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCaa                                                           <  2:140523/248‑1 (MQ=255)                                                                                                                                                                            gacctg‑atctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCa                                                            >  1:626983/4‑250 (MQ=255)                                                                                                                                                                               ctg‑atctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTAc                                                                                                                                                           >  2:424428/1‑152 (MQ=255)                                                                                                                                                                               ctg‑atctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTAc                                                                                                                                                           <  1:424428/152‑1 (MQ=255)                                                                                                                                                                                      tACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGCTCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATaccacc                                                    >  1:502867/1‑249 (MQ=255)                                                                                                                                                                                          ccTCGGAAGAAGCCGATGCGCCACAGTATATGGACCATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTTCCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATAGCTTCCATACCACCAGTGCg                                              <  2:111290/251‑1 (MQ=255)                                                                                                                                                                                           cTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGt                                             >  1:600848/1‑251 (MQ=255)                                                                                                                                                                                             cGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGc                                                                                                                          >  1:521074/1‑172 (MQ=255)                                                                                                                                                                                             cGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTgg                                           <  1:676904/251‑1 (MQ=255)                                                                                                                                                                                             cGGAAGAAGCCGAGGAGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGc                                                                                                                          <  2:521074/172‑1 (MQ=255)                                                                                                                                                                                                            cgcCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGACTATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGGGCGTGGGAATATGGTcgcgct                            >  2:336285/1‑249 (MQ=255)                                                                                                                                                                                                              ccACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGAtt                           >  1:90310/1‑250 (MQ=255)                                                                                                                                                                                                                acaGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATg                        >  1:19668/1‑251 (MQ=255)                                                                                                                                                                                                                   gTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCt                     >  1:630423/1‑251 (MQ=255)                                                                                                                                                                                                                         ggAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTca                <  1:202391/250‑1 (MQ=255)                                                                                                                                                                                                                          gAACACCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTcac               >  1:181362/1‑250 (MQ=255)                                                                                                                                                                                                                           aaCATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCacac             >  1:588537/1‑251 (MQ=255)                                                                                                                                                                                                                            aCATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGAAAACAAACTTGCCAATATCTTCAATACAACCAGTGCGTGGGAATATGGTCGCGATTATGTATTTTAACAcc            >  2:12817/1‑251 (MQ=255)                                                                                                                                                                                                                                  ggggATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTc          <  1:667311/247‑1 (MQ=255)                                                                                                                                                                                                                                     gATCGATACGGGGCTGGTCGATCCGAATCGTCCCTTTAGGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATAGATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTTGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGt    <  2:76345/250‑1 (MQ=255)                                                                                                                                                                                                                                      aTCGATACGGTGCTGGTCGATCCGAAACGTCCCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGCATATATTCCTCCGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTACATCCACCCTGGAAAACCAACTTGCCAATATCTTCACTACCACCATTGCGTGGGAATCTGGTCGCGCTTATGTCTTTTCACCCCCCGGCGGttg  >  2:247766/1‑249 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      CAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTG‑ATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGA  >  NC_000913/4627660‑4628138

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