Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,627,986	C→G	T224T (ACC→ACG)	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,986	0	C	G	100.0%	183.4 / NA	57	T224T (ACC→ACG)	nadR	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (27/30);  total (27/30)

AAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC  >  NC_000913/4627745‑4628231                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      aaGAGGGCATGGAGCCGTATCCGCACGGCTGTGATGTGGGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAg                                                                                                                                                                                                                                              <  2:563007/251‑1 (MQ=255)        cATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGACGTGAAGCCGttct                                                                                                                                                                                                                                       >  2:520433/1‑249 (MQ=255)               ccGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTtgtacg                                                                                                                                                                                                                                <  1:23708/251‑4 (MQ=255)                  tATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAACGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACCCCTCGGAAGGAGCCGCTGCGCCACAGAATATGGAACATCTGGGGATCGAAACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATACGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGttctttttacgtcc                                                                                                                                                                                                                             >  2:611019/1‑243 (MQ=255)                      cGCACGGCTGGGATGCGTGGAGAAACGGCATCAAAAAGTTTATGGCTGAACAGGGGATTCAGCCGGCCCTGCTCTACACCTCGGAAGAAGCCGATGCGCGACAGTATCTGGAACATCTGGGGCTCGAAAAGGTGCTGGTAGATCCGAAACGTACGTTTCTGAGTATCAGCGGTGCGCAGATCCGCGACAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGCTCTTTGTACGTACCGTg                                                                                                                                                                                                                         >  1:406243/1‑251 (MQ=255)                               gggATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAATCCGTTCTTTGTACGTACCGTGGCGATCCtt                                                                                                                                                                                                                <  1:119209/251‑1 (MQ=255)                               gggATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCtt                                                                                                                                                                                                                <  1:246801/251‑1 (MQ=255)                               gggATGTGTGGAGCAACGGCATAAAACAGTTTATAGCTGAAAAGGGGATTGAGCCGGACCTGATCTACACCTAAGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATAGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGGAAACCCGTTCCGCTACTGGGAATATATTCCTACGGCAGTGACGCCGTTCTTTGTACGTACCGTGACGATCCtt                                                                                                                                                                                                                >  2:423587/1‑251 (MQ=255)                                 gATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAAAGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTCCTGGGCATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTgg                                                                                                                                                                                                              >  1:381853/1‑251 (MQ=255)                                    gtggcgagtaaCGGAATCAGAAAGTTTATGGTTGATAAGGGGAATCAGCCGGACCTGATCTACACCTGGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTTAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTggtg                                                                                                                                                                                                            <  2:597657/241‑3 (MQ=255)                                             aaCGGCATCAAAAAGTTTATGGCTGAAAAGGGGGTTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCg                                                                                                                                                                                                    >  2:339640/1‑249 (MQ=255)                                                          aGTTGATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTAGCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCTAGCGGTAAATCCAcc                                                                                                                                                                                     <  1:520383/251‑1 (MQ=255)                                                           gTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGCAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACAGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCAccc                                                                                                                                                                                    >  1:359585/1‑251 (MQ=255)                                                           gTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCAccc                                                                                                                                                                                    >  1:404129/1‑251 (MQ=255)                                                                tGGCTGCAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACATTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGt                                                                                                                                                                                >  2:597992/1‑250 (MQ=255)                                                                          ggggATTCAGCCGGACCTGATCTACACCTCGGAAGCAGCCGATGCGCCACAGTATATGGCACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGAGGCGATCCTTGGTGGCGAGACGCGCGGGAAATCAACCCTGGTAAACAAACTTg                                                                                                                                                                     >  2:379554/2‑251 (MQ=255)                                                                                 cAGCCGGACCTTATCTACACCGCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCTGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCGACTGGGAATATATTCCTACGGAAGTTAAGCCGTGCTTTGTACGTACCGTGGCGATCCTGGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAatat                                                                                                                                                              <  1:565385/251‑1 (MQ=255)                                                                                   gCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCt                                                                                                                                                            <  1:188088/251‑1 (MQ=255)                                                                                   gCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACAGCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGAAGATCCGCGAAAACCCGTGCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACTTTCCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCt                                                                                                                                                            <  1:598677/251‑1 (MQ=255)                                                                                    ccGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCtt                                                                                                                                                           >  1:388023/1‑251 (MQ=255)                                                                                       gacctgatctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGGAAACAAACTTGCCAATATCTTCa                                                                                                                                                         >  2:168333/4‑250 (MQ=255)                                                                                                   acCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGc                                                                                                                                            <  1:339640/251‑1 (MQ=255)                                                                                                       cGGAAGAAGCCTATGCGCCACAGTATATGGAACATCTGGGGAGCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGAGCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTgg                                                                                                                                        <  2:225938/251‑1 (MQ=255)                                                                                                                ccGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATggt                                                                                                                               <  1:129937/251‑2 (MQ=255)                                                                                                                      cgcCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGAt                                                                                                                         <  1:168333/251‑1 (MQ=255)                                                                                                                        ccACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGAtt                                                                                                                        >  1:161208/1‑250 (MQ=255)                                                                                                                         cacaGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTAt                                                                                                                      <  1:619238/251‑1 (MQ=255)                                                                                                                            aGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCGTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTc                                                                                                                   <  2:551242/251‑1 (MQ=255)                                                                                                                               atatGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGAGCCGCGAAAACCCATTCCGCTACGGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCGCCCTGGTAAACAAACTTGCCAATATCTTCAAGACCGCCAGTGCGTGGGAATATGGGCGCGATTATGTCtt                                                                                                                 <  2:571692/250‑1 (MQ=255)                                                                                                                               atatGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTAGGAGTATCAGCGGTGCTCAGATCCGCGAAATCCCGTTCCGCTACGGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCTCCAGTGCGTGGGAATATGGTCGCGATTATGTCtt                                                                                                                 <  2:555860/250‑1 (MQ=255)                                                                                                                                 atGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTc                                                                                                              >  1:26752/1‑251 (MQ=255)                                                                                                                                 atGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTc                                                                                                              >  1:601712/1‑251 (MQ=255)                                                                                                                                    gAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGATTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGTGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTcaca                                                                                                           <  1:579423/251‑1 (MQ=255)                                                                                                                                            ggggATCGATACGGTGCTGGTCGATCCGAAACGTTCCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTc                                                                                                       <  1:623519/247‑1 (MQ=255)                                                                                                                                              ggATCGATACGGTGCTGGTCGATCCGAAACGTTCCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACGGGGAATTTATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGAGCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGTTCGCGATTATGTCTTTTCACTCCTCGGCGGt                                                                                                 <  2:397725/251‑1 (MQ=255)                                                                                                                                                           gCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCa                                                                                     <  1:41476/250‑1 (MQ=255)                                                                                                                                                              ggTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGc                                                                                 <  1:520433/251‑1 (MQ=255)                                                                                                                                                                           cGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTccgatt                                                                     >  1:85747/1‑245 (MQ=255)                                                                                                                                                                             tACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGCGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTTTTccgatta                                                                    >  1:538505/1‑243 (MQ=255)                                                                                                                                                                             tACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACg                                                                  >  1:136359/1‑251 (MQ=255)                                                                                                                                                                                  ttATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGGGATCGCATTGCAGTATTCCGATTACGATaa                                                              >  1:555636/1‑250 (MQ=255)                                                                                                                                                                                         tATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTTCGATTACGATAAAATCgcgc                                                      >  2:475905/1‑251 (MQ=255)                                                                                                                                                                                            cAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATcaaatcat                                                   <  1:328719/251‑8 (MQ=255)                                                                                                                                                                                                 gTGCGCAGATCCGCGACAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTAACGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGACGTATTCCGTTTCAGATAAAATCGCGCTGGGCCAc                                              >  2:53901/1‑251 (MQ=255)                                                                                                                                                                                                     gcAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTTCGATTACGATAAAATTGCGCTGGGCCAcccac                                          >  2:96748/1‑247 (MQ=255)                                                                                                                                                                                                       aGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAAGCGCGCTGGTCCACGCACaa                                        <  1:543910/251‑1 (MQ=255)                                                                                                                                                                                                            cgcgAAAACCCGTTCCGCTACGGGGAATATATTCCTACGGAAGGGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAAGACCACCAGTGCGTGGGAATATGGGCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACAt                                   <  1:598537/251‑1 (MQ=255)                                                                                                                                                                                                              cgAAAACCCGTTTCGCTACTGGGAATATATTCCTGCGGAAGTGAAGCCGTTCTTTGTACGGACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGAGTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCAGTGCAGTAGTCCGATTTCGATAAAAGCGCGCTGGGCCACGCACAATACATTg                                 <  2:61642/251‑1 (MQ=255)                                                                                                                                                                                                               gAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTgt                                >  1:375723/1‑250 (MQ=255)                                                                                                                                                                                                                aaaaCCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTTCGATTACGATAAAATCGCGCTGGGCCCCGCACAATACAtttat                               >  1:429332/1‑248 (MQ=255)                                                                                                                                                                                                                  aaCCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGAttt                             <  1:96748/251‑1 (MQ=255)                                                                                                                                                                                                                        ttCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACAcc                                                                                                         <  2:468340/169‑1 (MQ=255)                                                                                                                                                                                                                        ttCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACAcc                                                                                                         >  1:468340/1‑169 (MQ=255)                                                                                                                                                                                                                                    gAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCa            >  1:184955/1‑250 (MQ=255)                                                                                                                                                                                                                                        atatTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAAGACATTGATTTTGCAGTGAAATATGCCAATaaa       <  2:601712/251‑1 (MQ=255)                                                                                                                                                                                                                                            tCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTgg   >  1:569248/1‑251 (MQ=255)                                                                                                                                                                                                                                             ccTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATTCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGAGTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGc  <  1:53132/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      AAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC  >  NC_000913/4627745‑4628231

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