Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,313,027	T→G	intergenic (‑36/‑323)	alsR ← / → rpiB	DNA‑binding transcriptional repressor AlsR/allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,313,027	0	T	G	100.0%	154.1 / NA	42	intergenic (‑36/‑323)	alsR/rpiB	DNA‑binding transcriptional repressor AlsR/allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (15/27);  total (15/27)

AGCTTTGATACCTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGAAGGC  >  NC_000913/4312783‑4313263                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             agcTTTGATACCTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATcac                                                                                                                                                                                                                                         <  2:451751/250‑1 (MQ=255) agcTTTGATACCTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATcac                                                                                                                                                                                                                                         >  1:451751/1‑250 (MQ=255)     ttGATACCTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATcaccac                                                                                                                                                                                                                                      <  1:257519/249‑1 (MQ=255)      tGATACCTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAg                                                                                                                                                                                                                                    <  1:64971/250‑1 (MQ=255)      tGATACCTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCNTGAATTGAAGACGCCACTAGATCACCACAGc                                                                                                                                                                                                                                   >  2:223228/1‑251 (MQ=255)                   atcatcGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGaaa                                                                                                                                                                                                                        >  1:51262/1‑249 (MQ=255)                   atcatcGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGaaa                                                                                                                                                                                                                        >  1:42010/1‑249 (MQ=255)                       ccgcctcataTACCGCCAGTGCTTCTGCGACATCGTCAATTGCGGTTGCACAACTCAGGTTACCGGGTTTGAGTATCCTCTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTATCCCTAGACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATATCCACAGCAAGCACAGAAAATATGc                                                                                                                                                                                                                  <  2:62368/242‑1 (MQ=255)                            tCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGc                                                                                                                                                                                                              <  2:422447/250‑1 (MQ=255)                             cAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCg                                                                                                                                                                                                             <  1:232840/250‑1 (MQ=255)                               gATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACt                                                                                                                                                                                                          <  1:333740/251‑1 (MQ=255)                                   ccGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTcac                                                                                                                                                                                                      <  1:316247/251‑1 (MQ=255)                                    cGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTcaca                                                                                                                                                                                                     <  1:349124/251‑1 (MQ=255)                                     cccAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACaa                                                                                                                                                                                                    <  1:159185/250‑1 (MQ=255)                                       cAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACaaa                                                                                                                                                                                                   >  1:10742/1‑250 (MQ=255)                                                 gCGACATCTTTAATTGCGGGTGCTCAACTCAGGTTACCGGTTTTGAGTCACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGAACACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAaaaataat                                                                                                                                                                                           <  2:354315/248‑2 (MQ=255)                                                          ttAAGTGCTGGTGCACAACTCAGGTTACNGGGTTTGAGTAACCACTCCACGATGCTGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGTTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATaaa                                                                                                                                                                                  <  2:230882/248‑1 (MQ=255)                                                            aaTTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGttt                                                                                                                                                                              <  1:433240/250‑1 (MQ=255)                                                                 cGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGACTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCAGCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGattatttta                                                                                                                                                                        <  1:396480/251‑1 (MQ=255)                                                                  gggttcaaaaCTAAGGTTCCCGGGTTTGCGTGACCACTCCAAGATGCGGCTTTCA‑TGTTGTGTCATTCCTTCCTGGTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTGCGGAAGCGCTGAGTCAAACTCTTACTGGCTCATCACGGTCTCTTGTCCTTGAAGTGAAGACGCCACTAGATCACCACATCAAGAACATAAAAGATGCTGTCGACTTCACAAAAAAAGAAAAGTTATAAAGtatattttat                                                                                                                                                                       <  2:383496/243‑1 (MQ=255)                                                                          aaCTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAAc                                                                                                                                                                >  2:424634/1‑250 (MQ=255)                                                                           acccAGGTTACCGGGTTTGAGTACGCACTCCCAGATGCCACTTTCAGTTTTCT‑TTATTCCTTCCTGCTTCATTCGACGTTTAGGCTCTAACCCTATACCTGTCTGAAGCGCTGGATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTTAATTGAAGACGCCTCTAGATCACCACTACAAGAACAGAAAAGATGCTGGCGACTTTTCAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGTACaa                                                                                                                                                              <  1:435223/248‑1 (MQ=255)                                                                                                    cACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTAttt                                                                                                                                     <  2:42010/251‑1 (MQ=255)                                                                                                     aCTCCACGATTCGGCTTTCA‑TTTTCTGTCATTCCTGCCTGCTTCATTCGCAGGTAAGGCGCGAACCCTATACCGTTCGGAAGCGCTGAATCAAGCTCTGACTGGCTCATCTCGGTCTCTTGTCCGTGAAGTGAAGACGCCACTTGATCACCGGAGCTGGAACAGAAAAGATGCTGGCTACGTCGCAAAAAAATTAAAGTTATAAAGTTTAATTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTc                                                                                                                                    <  1:440510/251‑1 (MQ=255)                                                                                                        ccccGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCAtt                                                                                                                                 <  2:23599/248‑1 (MQ=255)                                                                                                             aTGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTattt                                                                                                                            >  2:341575/1‑248 (MQ=255)                                                                                                              ggctgCTTTCA‑GTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCCGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTa                                                                                                                           <  2:345214/247‑1 (MQ=255)                                                                                                                cGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTACTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTGCATTTAATTAGt                                                                                                                         <  1:255/251‑1 (MQ=255)                                                                                                                                   aTTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTAt                                                                                                       >  1:370284/1‑251 (MQ=255)                                                                                                                                            tGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTc                                                                                              >  1:177324/1‑251 (MQ=255)                                                                                                                                                tCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTg                                                                                            <  2:152459/249‑1 (MQ=255)                                                                                                                                                 cATTCGCATGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATATTATTTGTTTTTCTGGtt                                                                                         <  2:407894/251‑1 (MQ=255)                                                                                                                                                  aTTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGttt                                                                                        <  1:99314/251‑1 (MQ=255)                                                                                                                                                      gCAGGCAAGGCGCTAACCCTATACCGTTAGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGAGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTAATTTACTTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGtttttt                                                                                     >  2:272270/1‑250 (MQ=255)                                                                                                                                                              ggCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAg                                                                            <  1:223228/251‑1 (MQ=255)                                                                                                                                                                cgcTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTNATTTGTTTTTCTGGTTTTTTAgatgat                                                                              >  2:218487/1‑247 (MQ=255)                                                                                                                                                                                              aTCAAACTCTGACTGGCTCAACACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAACATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCt                                             >  1:327612/1‑250 (MQ=255)                                                                                                                                                                                                cAAACTCTGACTGGCTCATCACGTTCTCGTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATTGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTGTTTCTGGTTGTTTAGATGATAGATGGTGCTGGATGGTGATTGTCTTTAACGCTAt                                           <  2:449244/250‑1 (MQ=255)                                                                                                                                                                                                                             tgtacTTGAATTGAAGTCGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGAGAGTGGCTTAAAATTTTATTTCTTTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGATGGATGGTGAGTGTCTTTAACGCTNTTGCGCTGTCAGATGTAACTATGTAAAAc               >  2:220039/5‑249 (MQ=255)                                                                                                                                                                                                                                   tgaattgaaGACGCCACTAGATCACCACAGCAAGAACCGAAAAGATGCTGGCGACTTCACAAAAAAATAAACTTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTTTTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAgg       >  1:329942/1‑251 (MQ=255)                                                                                                                                                                                                                                        tgaaGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGAAgg   <  2:346295/250‑1 (MQ=255)                                                                                                                                                                                                                                        tgaaGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGAAGGc  >  1:445465/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             AGCTTTGATACCTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCA‑TTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGAAGGC  >  NC_000913/4312783‑4313263

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