Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,628,160	T→C	S282S (TCT→TCC)	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,628,160	0	T	C	100.0%	219.8 / NA	70	S282S (TCT→TCC)	nadR	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/30);  total (40/30)

TACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGG  >  NC_000913/4627917‑4628404                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     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)                                                        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)                                                                    acggAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACTCACAATACATTTATTTTGCAGTGAAATATGCCAATAAAGTGGcgtt                                                                                                                                                                            <  1:230320/248‑4 (MQ=255)                                                                    acggAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAAACGCGCTGGGCCACGCACAATACATTTGTTTTGCAGGGAAATATGCCACTAAAggtgcggt                                                                                                                                                                            >  2:14936/4‑244 (MQ=255)                                                                         aGTGACGCCGTTCTTTGTACGAACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCACTATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACAACTCGGCGGGGATGGAATCGCAATGCAGTATTCCGATTCCGATAAAATCGCGACGGGCCAAGACCAATACCTTGATTTTGCAATGAAATAAGCCacaacaggggcgtttattt                                                                                                                                                                       >  2:603970/1‑233 (MQ=255)                                                                         aGTGAAGCCGTTCTTTGTACGTACCGTGGCGAGCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGGAAACAAACGTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTAttg                                                                                                                                                                       <  1:22767/251‑3 (MQ=255)                                                                                ccGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGCTTACGGTAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTAttgaaactt                                                                                                                                                                 >  1:504760/1‑242 (MQ=255)                                                                                  gttctttGTCCGTACCGAGGCGATACTTGGTGGCGAGTCGCGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATCCCACCAGTGCGTGGGAATATGGTCGCGCTTATGTCTTTTCACACCTCGGCGGTGATGAACTCGCCATGCAGTATTTCGATTACGATCAAAACGCGCTGGGCCACGAAACATACATTGATTTTGAAGTGAAAAATTCCAATAAAGTGGCGTTTAttgtcaatgatt                                                                                                                                                              >  2:308118/5‑240 (MQ=255)                                                                                    tctttGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGAttt                                                                                                                                                             <  2:85747/248‑1 (MQ=255)                                                                                         gtacgttccgaggCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACCATACATTGAGTTTGCAGTGAAATATTCCAATAAAGTGCCGTTTATTGACCCTGATTTTGTct                                                                                                                                                        <  2:368552/239‑2 (MQ=255)                                                                                            cGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTcactac                                                                                                                                                    >  1:198565/1‑248 (MQ=255)                                                                                               aCCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTcactaccc                                                                                                                                                  <  1:475905/250‑6 (MQ=255)                                                                                                 cgtcgCGAACCTTGGTGGCGAGTCGAACGGTCAATCCACCCTGGTAAACAAACTTTTCAATATCTTCAATACCAACATTTCGTGGAAATATGGTCGCTATTATGTCTTTTCACACCTCGGCGGTGATGATATCGCATTGAAGTATTCCGACTACGATAAAATCGCGCTGGGCCACGCACAATACACTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTctctaccctgg                                                                                                                                               <  2:485653/247‑11 (MQ=255)                                                                                                          ccTTGGTGGCGAGTCGAGCGGTAGATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTGCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGACTACGATAAAATCTCGCTGGGCCCCGCCCAATTCATTTATTTTGCAGTGAAATATGCCATTAAAGTGGCTTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCa                                                                                                                                      <  2:381853/251‑1 (MQ=255)                                                                                                              ggtggcgAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCaaaa                                                                                                                                   <  1:134908/247‑1 (MQ=255)                                                                                                               gtggcgAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGGGGCGTTTATTTACACTGATTTTGTCACTACCCAGGCGTTCTGCaaaaa                                                                                                                                  >  2:189226/3‑250 (MQ=255)                                                                                                                      gTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAg                                                                                                                          <  2:161208/251‑1 (MQ=255)                                                                                                                        cGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAggg                                                                                                                        <  2:429332/251‑1 (MQ=255)                                                                                                                         gAGCGGTACATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGGGATGAGATCGCATTGCAGTATTTCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGAAAATGATTTTTTCACTACCCAGGCGCTCTGCAAAAAGTACGAAgggg                                                                                                                       >  2:157322/1‑250 (MQ=255)                                                                                                                         gAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCCTTGCAGTATTTCGTTTTCGATAAAATCGCTCTGGGCCACGCCCAATACATTTATTTTGCAGTGAAATATGCCAATTAAGTGGCGTTTATTTACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAATTACGAAGGGc                                                                                                                       >  2:416604/1‑251 (MQ=255)                                                                                                                          aGCGGTAAAGCCGCCCTGGTAAGCAAACTTGCCAATATCTTCAATCCCACCAGTGCGTTGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTTCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGAGTTTGTCACTACCCAGGCGTTCTGCAAAAAGGACTAAgg                                                                                                                         <  2:251242/248‑1 (MQ=255)                                                                                                                             ggTAAATCCACCCTGGTAAACAACCTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTAAGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGGGAAATATGCCAATAAAGTGGCGTTTATTGACACCGATTTTTTCACTAACCAGGCGTTCTGCAAAAAGAACGAAGGGCgtg                                                                                                                    >  1:171481/1‑248 (MQ=255)                                                                                                                             ggTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCgtg                                                                                                                    >  1:266671/1‑248 (MQ=255)                                                                                                                              gTAAATCCTCCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTTCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACTATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATGGACACTGATTTTGGCACTACCCAGGCGTTCTGCAAGAAGTACGAAGGTCGGGag                                                                                                                  <  2:603462/251‑2 (MQ=255)                                                                                                                              gTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGAAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGAAAATGATTTTGTCACTACCCAGGGGTTCTGCAAAAAGTACGAagggaggg                                                                                                                    >  1:72499/1‑249 (MQ=255)                                                                                                                                  aTCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATc                                                                                                              >  1:244427/1‑251 (MQ=255)                                                                                                                                           ggTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGt                                                                                                       <  2:69856/249‑1 (MQ=255)                                                                                                                                                    acTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCCGTATTCCGGTTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGCTTTTGTCACTACCCAGGCGTTCTGCCAAAAGTACGAAGGGCGTGAGCATCCGTTCGtcaaggcgttg                                                                                             >  2:195227/1‑240 (MQ=255)                                                                                                                                                                   caataccaccAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGGTGAGATCGCATTGCAGTATTCCGAATACGATAAAATCACGCTTGGCCACGCACCATACCTTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGGTTTTTGACACTAATTTTTTAACAACCCAGGCGTTCTGCAAAAAGTACAAAGGGAGCGAGCATCCGTTCGTACAGGCGTTGATTGGTGAATACCGtt                                                                             >  1:112635/5‑251 (MQ=255)                                                                                                                                                                    aataccaccAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGttt                                                                            <  1:379554/248‑1 (MQ=255)                                                                                                                                                                          accaGTGCGTGGGAATCTGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCCGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGCTTATTGACACTGATTTTGTCACTAACCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCt                                                                      >  1:531537/1‑251 (MQ=255)                                                                                                                                                                               tGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGGGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATCAAGTGGCGTTTATTGACACTGATTTTGTCACAAACCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAAAAACGTTTCGATCTGGTg                                                                  >  1:424677/1‑250 (MQ=255)                                                                                                                                                                                  gTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATcc                                                              >  2:133185/1‑251 (MQ=255)                                                                                                                                                                                     ggAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAATATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGTTTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGTTTTTGGCACTACCCAGGCGTTCTGCAAAAAGTACCAAGGGCGGGAGCATCCGTTCGGACAGGCGTTGATTGATGAATACCGTTTCGGTCTGGGGATCctcg                                                           >  2:186740/1‑249 (MQ=255)                                                                                                                                                                                     ggAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGACTTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGAGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACCGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCctgc                                                           >  1:394516/1‑251 (MQ=255)                                                                                                                                                                                        atatGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGATCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTgg                                                        >  1:46073/1‑251 (MQ=255)                                                                                                                                                                                                  gATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCAGGGCCACGCACAATACATTGATTTTGCAGGGAAAAATGCCAATAACGGGGCGTTTAATGACACTGATTTTGTCACTACCCAGGCGTTCTGAAAAAAGTACGGAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGCGATCCTGCTGGAGAACAacac                                              >  1:54388/1‑251 (MQ=255)                                                                                                                                                                                                    ttATGTCTTTTCACACCTCGGCGGTGATGCGATCGCATTTCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGAAGTGAAATATGCCAATAAAGTGGGGTTTATTGAAAATGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGAACGAAGGGCGGGAGCATCCGTACGGAAAGGCGGTGTTTGGTGAATACAGTTTCGATCTGGTGATCCTGCTGGAGAACCACACGc                                            >  2:225840/1‑251 (MQ=255)                                                                                                                                                                                                            tttCACACCTCGGCGGTTATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCGGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCgtgg                                       <  1:514479/248‑1 (MQ=255)                                                                                                                                                                                                                     tCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCCCAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGGGGGGGGCGGATggg                           >  1:373615/1‑250 (MQ=255)                                                                                                                                                                                                                      cggcggTGAGGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACTCCGTGGGTGGCGGATGGtt                          <  2:500478/251‑1 (MQ=255)                                                                                                                                                                                                                       ggcggTGATGAGATCGCATTGCATTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGttt                         <  2:266671/251‑1 (MQ=255)                                                                                                                                                                                                                       ggcggTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACCTTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGttt                         >  1:589436/1‑251 (MQ=255)                                                                                                                                                                                                                         cggTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATTCATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTGATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGGATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTTGGTGGCGTATGGTTTAc                       <  2:105205/251‑1 (MQ=255)                                                                                                                                                                                                                           gTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGc                     >  1:474966/1‑251 (MQ=255)                                                                                                                                                                                                                            tgatgaGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTtgatt                                                                                          >  1:6614/1‑181 (MQ=255)                                                                                                                                                                                                                            tgatgaGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTtgatt                                                                                          <  2:6614/181‑1 (MQ=255)                                                                                                                                                                                                                                gagaTCGCATTGCAGTATGCCGATTACGATAAAATCGCGCTGGGCCACGCACAATTCATTTATTTTGCAGTGAAATATGCCAATAAAGTTGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCTTGAGCATCCGTTCGTACAGGCGTTGATTGAGGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGTCGGATGGTTTACTCAGCCt                <  2:388023/251‑1 (MQ=255)                                                                                                                                                                                                                                     cGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCa           <  2:184955/251‑1 (MQ=255)                                                                                                                                                                                                                                      gCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGc            >  1:164844/1‑249 (MQ=255)                                                                                                                                                                                                                                       cATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTAGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGCGGGTGGCGGCTGGTTTACGCCGCCTCGGCAGt         >  2:464843/1‑251 (MQ=255)                                                                                                                                                                                                                                        aTTGCAGTATTCCGATGACGATAAAATCGCGCTGGGCCACGCACAAGACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGtt        <  2:404129/251‑1 (MQ=255)                                                                                                                                                                                                                                        aTTGCAGTAATCCGATTACGATAAACTCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGtt        >  2:328956/1‑251 (MQ=255)                                                                                                                                                                                                                                         ttGCAGTATTCCGATTACGATACAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTc       >  1:160325/1‑251 (MQ=255)                                                                                                                                                                                                                                          tGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCAAAATACATTGATTTTGCAGTGAAATATGCCAATAAAGAGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTATGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGAGGGAGGCGGATGGTTTACGCAGCCTCGGCAGTTCg      >  2:163207/1‑251 (MQ=255)                                                                                                                                                                                                                                              gTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTgg  <  1:133185/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     TACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGG  >  NC_000913/4627917‑4628404

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