Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	431,516	A→C	Y166* (TAT→TAG)	tsx ←	nucleoside‑specific channel‑forming protein Tsx

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	431,516	0	A	C	100.0%	145.8 / NA	39	Y166* (TAT→TAG)	tsx	nucleoside‑specific channel‑forming protein Tsx
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (19/20);  total (19/20)

AGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCATAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTCCATAAACAGCGGAGAACCGTGGTTCCAGATACCTTTAGCATCGG  >  NC_000913/431269‑431757                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  aGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAg                                                                                                                                                                                                                                                <  1:122517/251‑1 (MQ=255)              cAGAGCCAGAATATGGCTGGACGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGCTCGGTAATCGGCACAAAGGATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGAc                                                                                                                                                                                                                                    >  1:403106/1‑250 (MQ=255)              aagagCCAGAATATGGCTGGACGCGATAGAGTTATGAGTACGGGTCGTAATACCGTTGATTGCGTTACCGCTTGCATCTCCGAAATCGGAACCCCAGTCGAAGTTGGTGAAGCAGATGTAGCTCAGATTACCGCCCCACAGATCGGTAATCGGCACAAAGTATTGAATTTTGAAACGGTAACCGTCCCACTCGTTGTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACa                                                                                                                                                                                                                                   <  1:245785/250‑1 (MQ=255)                gagCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCTCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAgg                                                                                                                                                                                                                                 <  2:369256/251‑1 (MQ=255)                    cAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCACTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGccc                                                                                                                                                                                                                              >  2:345247/1‑248 (MQ=255)                                      aaTAGAGTTATTAGTCCGGGTCTTAATACCGTTGATTGCGTTACCGCTGCCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGCAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAAACGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCCGGCCAGGGTCg                                                                                                                                                                                                            >  2:404460/2‑250 (MQ=255)                                          gagtgaataGTATGGGTATGAATACCGGTGATTGCGTTCCCTCGGTCATCCCCTAAATCGGAACCCCAGTCGTAGATGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCTGCACAAAGTATTTAATTTTGAAACGGTGATCGTCCCACTCTTTTTCGTTCGCTTCGTCATAGTTCTGCCTCTGGTAGTTTGCCTAGACGTTCATGGACAGGCTCATCTGTAGTCCATTGTAGATATc                                                                                                                                                                                                       <  2:376396/244‑1 (MQ=255)                                           agTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCg                                                                                                                                                                                                      <  1:243404/251‑1 (MQ=255)                                                 taGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTAccc                                                                                                                                                                                                >  1:209058/1‑251 (MQ=255)                                                         gTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCAt                                                                                                                                                                                        <  2:143958/251‑1 (MQ=255)                                                             tAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATg                                                                                                                                                                                       <  2:221597/248‑1 (MQ=255)                                                                                          tCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACg                                                                                                                                                       >  1:119000/1‑251 (MQ=255)                                                                                                aaaTCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCAt                                                                                                                                                 <  2:14676/251‑1 (MQ=255)                                                                                                    cGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCg                                                                                                                                             <  1:215160/251‑1 (MQ=255)                                                                                                     ggAACCCCAGTCGAAGTTGGTGAAGCCGCTGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCTACCATCTTTATTACGACCCATGTCGt                                                                                                                                            >  2:205663/1‑251 (MQ=255)                                                                                                              gTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAgt                                                                                                                                   <  1:221010/251‑1 (MQ=255)                                                                                                              gTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAgt                                                                                                                                   >  2:347873/1‑251 (MQ=255)                                                                                                                cGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAgttg                                                                                                                                 >  1:286949/1‑251 (MQ=255)                                                                                                                           aaGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCCGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGt                                                                                                                       >  2:248472/1‑250 (MQ=255)                                                                                                                                gATGTAGCTCAGCTGACCGCCCCCCAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCCTCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTAcccct                                                                                                                 >  2:266529/1‑248 (MQ=255)                                                                                                                                 aTGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCActc                                                                                                                <  1:83432/251‑1 (MQ=255)                                                                                                                                    tAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCttt                                                                                                             >  2:71590/1‑251 (MQ=255)                                                                                                                                              gACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGAc                                                                                                     <  1:105788/249‑1 (MQ=255)                                                                                                                                                    cccACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCg                                                                                                                                             >  1:153861/1‑203 (MQ=255)                                                                                                                                                    cccACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCg                                                                                                                                             <  2:153861/203‑1 (MQ=255)                                                                                                                                                                    cGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTc                                                                              >  1:190941/1‑250 (MQ=255)                                                                                                                                                                            aGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCg                                                                     <  2:247198/251‑1 (MQ=255)                                                                                                                                                                               aTTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACAAGGTGCTCTGGCGACCACCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGAt                                                                   <  1:178083/250‑1 (MQ=255)                                                                                                                                                                                 ttAATTTTGAAACGATAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATgg                                                                 >  1:312160/1‑250 (MQ=255)                                                                                                                                                                                  tAATTTTGAACCGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCgcagg                                                                 >  2:268995/1‑245 (MQ=255)                                                                                                                                                                                   aaTTTTGAAACGGTAACCGGCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCGCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGGTAAGGTCAGTAGTGGTCAGCTTGTCGATGgaga                                                              <  2:197009/251‑1 (MQ=255)                                                                                                                                                                                    aTTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGaa                                                             <  1:226796/251‑1 (MQ=255)                                                                                                                                                                                         gAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCGTTGAACGGACCGAAGCGAAGGTCAGTATTGGTCAGCGTGTCGATGGAGAAACGTg                                                        <  2:45412/251‑1 (MQ=255)                                                                                                                                                                                            aCGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATCTCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGtt                                                     >  2:374750/1‑251 (MQ=255)                                                                                                                                                                                                    cGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCTGTATTTGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTcc                                             >  2:76209/1‑251 (MQ=255)                                                                                                                                                                                                    cGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTGCGATTTcc                                             <  2:23916/251‑1 (MQ=255)                                                                                                                                                                                                                                            cTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGAcc                                                                                                    <  1:69340/156‑1 (MQ=255)                                                                                                                                                                                                                                            cTGGTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGAcc                                                                                                    >  2:69340/1‑156 (MQ=255)                                                                                                                                                                                                                                               gTATTTCGCCTAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTCCATAAACAGCGGAGAACCGTGGTTCCAGATACCTTTAGCATCgg  >  1:185377/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  AGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAACCGTCCCACTCGTTTTCGTTCGCTGCGCCATAGTTCTGCCACTGGTATTTCGCATAGACGTTCATGGACAGGCTCATCGGCAGGCCAGTGTCGATATCGGTACCCAGACCCATGTACCAGGTGCTCTGGCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTCCATAAACAGCGGAGAACCGTGGTTCCAGATACCTTTAGCATCGG  >  NC_000913/431269‑431757

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