Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,098,150	G→A	D139N (GAT→AAT)	rhaS →	DNA‑binding transcriptional activator RhaS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,098,150	0	G	A	100.0%	164.2 / NA	46	D139N (GAT→AAT)	rhaS	DNA‑binding transcriptional activator RhaS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/26);  total (20/26)

CATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATCTGTGTCTGACCAATGTGCTGTATCGCTCGCCGGATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACATCGGCAGCTTAAGCAGCAAACGGGACTGACGCCTCAGCGATACCTGA  >  NC_000913/4097903‑4098396                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       cATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATCTGTGTCTGACCAATGTGCTGTATCGCTCGCCGGATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAAtt                                                                                                                                                                                                                                                     <  1:55301/251‑1 (MQ=255)               cGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATCTGTGTCTGACCAATGTGCTGTATCGCTCGCCGGATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGGGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGcc                                                                                                                                                                                                                                       >  1:145690/1‑251 (MQ=255)                                      gATCGGCATCTGTATGAACATACCGATAATCTGTGTCTGACCAATGTGCTGTATCGCTCGCCGGATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGAAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGc                                                                                                                                                                                                                 >  2:311065/1‑250 (MQ=255)                                         cGGCATCTGTATGAACATACCAATAATCGGTGTCTGACAATTGTGCTGTATCGCTCGCCGGATCGATTTAAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAAAAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAAt                                                                                                                                                                                                              <  1:269779/250‑1 (MQ=255)                                                  tATGACCAGACCGATAATCTGTGTCTGACCAAGGTGCTGTATCGCTCGCCGGATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGAGGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGctct                                                                                                                                                                                                     <  1:222112/250‑1 (MQ=255)                                                   aTGAACATACCGATAATCTGTGTCTGACCAATGTGCTGTATCGCTCGCCGGATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGAAAAATAATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtg                                                                                                                                                                                                   <  1:125808/251‑1 (MQ=255)                                                                tAATCTGTGTCTGACCAATGTGCTGTATCGCTCGCCGGATCGATTTCAGTTTCTCGCCGGACTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTaaacgc                                                                                                                                                                                       >  2:188494/1‑247 (MQ=255)                                                                 aatttgtgTCTGACCAATGTCCTGTATCGCTCGCCGGATCGATTTCAGTTCCCCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGAGGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCATATTGCAGCAGGTGCGACAGCTGGTTGCACAGAGGGAACAGCAGGAAGGGAAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAg                                                                                                                                                                                     <  2:256150/247‑1 (MQ=255)                                                                    cggtgtCTGACCAATGTGCTGTATCGCTCGCCGTATCTATTTCATTTTCTCGCCGGGCTGATTCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCAAAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTTCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGtt                                                                                                                                                                                   <  1:23964/248‑1 (MQ=255)                                                                     ggtgtCTGACCAATGTGCTGTATCGCTCGCCGGATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCAACAAGAGCGGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTg                                                                                                                                                                                 <  2:129564/250‑1 (MQ=255)                                                                                aaTGTGCTGTTCCGCTCGCCGGATCGTTTTCATTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGAGGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCGGGGGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCt                                                                                                                                                                      <  2:347894/251‑1 (MQ=255)                                                                                   gtgCTGTATCGCTCGCCGGATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTgg                                                                                                                                                                    >  1:120384/1‑250 (MQ=255)                                                                                        gTATCGCTCGCCGGATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTCACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACa                                                                                                                                                              <  1:282526/251‑1 (MQ=255)                                                                                         tATCGCTCGCCGGATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAAc                                                                                                                                                               >  1:195789/1‑249 (MQ=255)                                                                                                 gcCGGATTGATTTCAGTTTCCCCCCGGGCTGAATCAGTGGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCAc                                                                                                                                                     <  2:91672/251‑1 (MQ=255)                                                                                                 gcCGGATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCAc                                                                                                                                                     >  1:128420/1‑251 (MQ=255)                                                                                                    ggATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCCTCACGtc                                                                                                                                                  >  1:235947/1‑251 (MQ=255)                                                                                                    ggATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGtc                                                                                                                                                  >  1:321442/1‑251 (MQ=255)                                                                                                    ggATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGtc                                                                                                                                                  >  1:59199/1‑251 (MQ=255)                                                                                                    ggATCGATTTCAGTTTCTCGCAGGGCTGAATCAGTTGCTGCCACAAGAGCTGGCTGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTACAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCCTCACGtc                                                                                                                                                  >  2:306773/1‑251 (MQ=255)                                                                                                      aTCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGAAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGtctc                                                                                                                                                >  2:28771/1‑251 (MQ=255)                                                                                                                          ggtcTGAATCCGTTGCTGCCACAAGAGCTGGATGGGCATTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACATCAGGAAGGGTAAAATAATTTTCCCTCTCCCGCCAGTCGCGAGATCTTTTTTATGCAATTACTGCTCTTGCTGCGTAAAATCATTTTGCAGGATAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggct                                                                                                                            <  1:175215/248‑1 (MQ=255)                                                                                                                          gggCGGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTAAACCACAGCGTATTGCAGCGGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTTCTTCTGGCctggct                                                                                                                            <  1:296431/251‑1 (MQ=255)                                                                                                                            gcggAATCAGTTGCGGCCACAAGAGCTGGAGGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctggcctgtcagg                                                                                                                          <  2:321442/248‑6 (MQ=255)                                                                                                                            gCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAg                                                                                                                                                                                     >  2:71090/1‑192 (MQ=255)                                                                                                                            gCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAg                                                                                                                                                                                     <  1:71090/192‑1 (MQ=255)                                                                                                                              tGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggctgg                                                                                                                          >  1:169739/1‑249 (MQ=255)                                                                                                                               gAATCAGTGGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCGGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTggagg                                                                                                                       <  1:282952/251‑1 (MQ=255)                                                                                                                                    ggtggctgcCAAAAGAGCTGGATGGGCAGTACCCGTCCCCCTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCAt                                                                                                                  <  2:107323/247‑1 (MQ=255)                                                                                                                                        gctgcaaaaaGAGCGGGATGGGAAGTATCCGTCTCACGGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCAtttt                                                                                                               <  2:275420/242‑1 (MQ=255)                                                                                                                                                    gcgggATGGGCAGTTTCCGTCTCACTGGCGCGTTAACCACAGCGTTTGGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTCCCCCCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGa                                                                                                  <  2:111801/248‑1 (MQ=255)                                                                                                                                                        gATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTg                                                                                              >  2:195782/1‑251 (MQ=255)                                                                                                                                                           gggCAGTATCCGTCCCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGa                                                                                           <  2:127770/251‑1 (MQ=255)                                                                                                                                                                  aTCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTg                                                                                    <  1:195782/251‑1 (MQ=255)                                                                                                                                                                      gTCTCCCTGGCGCGTTCACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCCGGAAGGGGCAAATACTTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGCACCTGGAAAACAGCGCATCCCGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGCTGAGGCGAATTGGGATGCCGTGGCgg                                                                                >  2:220128/1‑251 (MQ=255)                                                                                                                                                                              ggggcgTTAACCACAGCGTATTACAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAGGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAAt                                                                          <  2:30931/246‑1 (MQ=255)                                                                                                                                                                                     tAACCACAGCGTATTGCAGCAGTTGCGACAGCTGGTGGCACAGATGGAACAGAAGGAAGGTGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCtt                                                                  <  2:25665/250‑1 (MQ=255)                                                                                                                                                                                                  ttGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCt                                                    >  1:38308/1‑251 (MQ=255)                                                                                                                                                                                                      agcagGTGCGAAATCTGGTTGCACAGATGGAACAGCAGGAAGGGAAAAATAATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAt                                                <  1:10582/251‑1 (MQ=255)                                                                                                                                                                                                         agGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACCGTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGCACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACATCgg                                             >  1:24159/1‑251 (MQ=255)                                                                                                                                                                                                                acatcgggTTGCCCAGAGGGAACAGCAGGAAGGGGAAAAGAATTTGCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTGGCGTCTGGCCTGGCTGGAGGACCAGTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTTCGTACGCTACATCGGCAGCTTa                                      <  1:357119/245‑1 (MQ=255)                                                                                                                                                                                                                    cTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATAATTTACACTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACATCGGCAGCTTAagc                                   >  1:221751/1‑250 (MQ=255)                                                                                                                                                                                                                       gtggCACAGTTGGAACATCAGGAAGGATAAAATAATTTCCCCCCGACCGCCAGTCGCGAGATCTTGTGTATGCAATTACTGCGCTTGCTGCGTAAAAGCAGTTTGCATGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGTCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTCCAGCGGCTGCTTATgcagca                               <  2:209302/248‑1 (MQ=255)                                                                                                                                                                                                                                          agGAAGGGGAAACTAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAACACCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACATCGGCAGCTTAAGCAGCAACCGGGACTGACGCCTccg            >  2:350264/1‑249 (MQ=255)                                                                                                                                                                                                                                                  gAAAATAATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACATCGGCAGCTTAAGCAGCAAACGGGACTGACGCCTCAGCGATACCt    >  1:271951/1‑251 (MQ=255)                                                                                                                                                                                                                                                    aaaTAATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATTCCGTGGCGGATCAATTTTCGCTTTCACTGCGTACGCTACATCGGCAGCTTAAGCAGCAAACGGGACGGACGCCTCAGCCAGACCTGa  <  2:230807/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       CATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATCTGTGTCTGACCAATGTGCTGTATCGCTCGCCGGATCGATTTCAGTTTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACATCGGCAGCTTAAGCAGCAAACGGGACTGACGCCTCAGCGATACCTGA  >  NC_000913/4097903‑4098396

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