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	95.7%	139.9 / ‑0.7	47	E193D (GAG→GAT)	nadR	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (25/20);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.82e-01

GCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGA‑TCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGT  >  NC_000913/4627659‑4628136                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    gcagcaGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAACCTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGCGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCgctc                                                                                                                                                                                                                                     >  2:141205/1‑248 (MQ=255) acgacgattttccgatctggATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATc                                                                                                                                                                                                                                     <  2:222744/233‑1 (MQ=255)  cagcagCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGCAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATcc                                                                                                                                                                                                                                    >  1:397446/1‑251 (MQ=255)      agcCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAGGGATCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGaaa                                                                                                                                                                                                                                <  1:479447/251‑1 (MQ=255)         ccaCCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTACATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGCAGCCGATGCGCCACCGTATATTGACCATATGGGGATCGCTACGGTGATGGTCGATCCGacaagt                                                                                                                                                                                                                             >  2:427413/3‑246 (MQ=255)          caccttGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTa                                                                                                                                                                                                                            <  1:135070/246‑1 (MQ=255)                   ggATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATc                                                                                                                                                                                                                                     >  1:222744/1‑233 (MQ=255)                          cTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATAGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGTTCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATc                                                                                                                                                                                                              >  2:83939/1‑249 (MQ=255)                            gcgctGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAg                                                                                                                                                                                                            <  2:388537/245‑1 (MQ=255)                                  gTTATTACAAACTTTTAAATATCAGAAAAATCTTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTgcgc                                                                                                                                                                                                    >  2:116136/1‑251 (MQ=255)                                        acAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TGCAGCCCGCCCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATc                                                                                                                                                                                               <  2:477519/249‑1 (MQ=255)                                                        cAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCAGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAACCGGACCTTTATGTGTATCAGCGGTGCGCAGCTCCGCGAAAACCCGATCCg                                                                                                                                                                              >  2:122377/1‑251 (MQ=255)                                                             aaaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACt                                                                                                                                                                         >  1:80947/1‑251 (MQ=255)                                                               atatTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAGAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTgg                                                                                                                                                                       >  1:343130/1‑251 (MQ=255)                                                                         acATGATTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGCAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCTATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGATAACCCGTTCCGCTACTGGGATGATATTc                                                                                                                                                             >  2:95368/2‑251 (MQ=255)                                                                                tgcAACGAAGAGGGTATGGAGCCGGATCCGCACTGCGGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATAGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAGCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCTGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTAcgg                                                                                                                                                       <  1:222220/248‑3 (MQ=255)                                                                                                   aGCCGTTTCCGCCCGGCTGGGATGGGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTtgta                                                                                                                                   <  2:397446/251‑2 (MQ=255)                                                                                                           ccGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGt                                                                                                                           >  1:400409/1‑251 (MQ=255)                                                                                                                  gCTGGGATGTGTGGAGCAACGGCATAAAAAAGTTTATGGCTGACAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCCGGGGATCGCTACGGAGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAAACCGTTCCGCTAATGGGAATATATTCCTACGGAAGTGAAGACGTTCTTTGTACGTACCGTgg                                                                                                                         >  2:73849/1‑246 (MQ=255)                                                                                                                                      ggCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCttt                                                                                                                                                                                                                       <  2:115674/132‑1 (MQ=255)                                                                                                                                      ggCATCAAAAAGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCttt                                                                                                                                                                                                                       >  1:115674/1‑132 (MQ=255)                                                                                                                                        cATCAAAAAGGTTGTGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATCTGGAACATCTGGGGATCGATACGGTGCTGGTCGCTCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGAGCAGCCGTTCTTTGGACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCggg                                                                                               >  2:208931/1‑249 (MQ=255)                                                                                                                                              aaaGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATcc                                                                                                                                                                                                                                    >  1:185700/1‑111 (MQ=255)                                                                                                                                              aaaGTTTATGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATcc                                                                                                                                                                                                                                    <  2:185700/111‑1 (MQ=255)                                                                                                                                               aaGTTTATGGCTGAAAAAGGGA‑TCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATacac                                                                                       >  2:58279/1‑251 (MQ=255)                                                                                                                                               aaGTTTATGGCTGAAAAAGGGA‑TCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCAc                                                                                       <  1:153531/251‑1 (MQ=255)                                                                                                                                                   ttATGGCTGAAAAGGGGA‑TTCAGCCGGACATGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTCCCTTTATGAGTATCCGCGGTGCGCAGATCCGCGAAAACCCGATCCGCTACTGGGAATATATTCCTACGGAAGGGAAGCCGTTCTTTGTACGTAACGTGGCGGTCCTTGGTGGCGAGTGGAGAGGAAGATCCACCCTg                                                                                   >  2:472683/1‑251 (MQ=255)                                                                                                                                                     aTGGCTGAAAAGGGGA‑TTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGt                                                                                 <  1:83939/251‑1 (MQ=255)                                                                                                                                                            gaaaaGGGAATTCAGCCGGACCTGCTCTACACCTCGGAAGAAGCCGATGCGCAACAGTATATGGAACATCTGGGGATCGCTACGGTGCTGGTCGATCCGCAAAGTCCCTTTATGAGTATCAGCGGTGCGCAGATCCGCGACAACCCGTTACGCTACAGGGCATATATTCCTACGGACGTGAAGCCGTTCTTTGTACGTACCGTGGCGCTCATTGGTGGCGAGTCGATCGGTAAATCAACCCTGGAaaaca                                                                            >  2:412514/2‑250 (MQ=255)                                                                                                                                                               aggggA‑TTCAGCCGGACCGGATCTACTCCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACtt                                                                       <  2:398325/249‑1 (MQ=255)                                                                                                                                                                    gagtcAGCCTGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCa                                                                   <  2:379211/247‑1 (MQ=255)                                                                                                                                                                                   gatctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCATTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCaatac                                                      <  1:141205/249‑4 (MQ=255)                                                                                                                                                                                        acacCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGttc                                                                                                                                         <  2:142099/161‑2 (MQ=255)                                                                                                                                                                                        acacCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGttc                                                                                                                                         >  1:142099/1‑160 (MQ=255)                                                                                                                                                                                          acCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAATCCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGc                                             >  1:375063/1‑251 (MQ=255)                                                                                                                                                                                               ggAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTg                                                                      >  1:409997/1‑221 (MQ=255)                                                                                                                                                                                               ggAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTg                                                                      <  2:409997/221‑1 (MQ=255)                                                                                                                                                                                                    aagCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCt                                                             >  1:174273/1‑225 (MQ=255)                                                                                                                                                                                                    aagCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCt                                                             <  2:174273/225‑1 (MQ=255)                                                                                                                                                                                                                   aGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAAACCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGAGCGTGGGAATATGGTCGCGATTATgac                    >  2:322976/1‑249 (MQ=255)                                                                                                                                                                                                                    gTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGAGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGTTTATGTCt                   >  2:309822/1‑251 (MQ=255)                                                                                                                                                                                                                        atGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGAGTATGTCttgt                <  2:80947/250‑3 (MQ=255)                                                                                                                                                                                                                              cATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGCAAACCCGCTCCGCTACTGGGAATATATTCCTACGGAAGTGCAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATCCCACCAGAGCGTCGGAATATGGTCGCGATTATTTCTTTTCACAccc         >  2:127109/1‑250 (MQ=255)                                                                                                                                                                                                                              cATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCt         <  1:378398/251‑1 (MQ=255)                                                                                                                                                                                                                                 cTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTcgg      >  1:128399/1‑251 (MQ=255)                                                                                                                                                                                                                                  tGGGGACCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTcgg      <  2:19079/250‑1 (MQ=255)                                                                                                                                                                                                                                     ggATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGt  <  1:448256/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    GCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGA‑TCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGT  >  NC_000913/4627659‑4628136

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: ‑