Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,467,674	T→G	L17R (CTG→CGG)	mgtA →	Mg(2(+)) importing P‑type ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,467,674	0	T	G	100.0%	128.5 / NA	39	L17R (CTG→CGG)	mgtA	Mg(2(+)) importing P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (21/18);  total (21/18)

GAACCTGATCCCACGCCTCTCCCTCGACGGAGATTAAAACTTTTCCGGTAAGCCCGTCTTTTCACGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCA  >  NC_000913/4467434‑4467919                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      gAACCTGATCCCACGCCTCTCCCTCGACGGAGATTAAAACTTTTCCGGTAAGCCCGTCTTTTCACGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGt                                                                                                                                                                                                                                             >  1:422725/1‑251 (MQ=255)  aaCCTGATCCCACGACTATCCCTCGACGGAGATTAACACTTTTCCGGTAAGCCCGTCTTTTAACGGCGTTACCGCATGCGTAAGGACGTGCCGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGACATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTg                                                                                                                                                                                                                                            >  2:208475/1‑251 (MQ=255)                                       aCTTTTCCGGTAATCCAGTCTTTTCACGGCGTTACCGGATGCGGAAGGCCGTGACGTTTTAACGTCCCTGCTCAGATTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGAAAGTGCGGAGGGACTCCTTATGTTTTAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATACATTGCCTGCCGCGCAGCAGACCGTGAATACggt                                                                                                                                                                                                       >  2:331695/1‑251 (MQ=255)                                          tttCCGGTAAGCCCGTCTTTTCACGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACggtggt                                                                                                                                                                                                    >  1:198356/1‑251 (MQ=255)                                                   aGCCCGTCTTTTCCCGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTc                                                                                                                                                                                           <  1:338191/251‑1 (MQ=255)                                                   aGCCCGTCTTTTCAGGGCGTTACCGGATGCGGAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTTCTGTTCCAGACAGTGCGGAGGGACTCCTTATTTTGAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTc                                                                                                                                                                                           <  2:285641/251‑1 (MQ=255)                                                                                        cGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGGGATGccc                                                                                                                                                      >  1:240542/1‑251 (MQ=255)                                                                                           gACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCCGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGGGATTCCCga                                                                                                                                                    >  2:325829/1‑250 (MQ=255)                                                                                                tttAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaag                                                                                                                                               <  1:346024/250‑1 (MQ=255)                                                                                                   aaCGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAtt                                                                                                                                           <  1:374229/251‑1 (MQ=255)                                                                                                         ccTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCTGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAGGAAGAATTGTGGaa                                                                                                                                     <  1:82273/251‑1 (MQ=255)                                                                                                                           cTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGCAGAAGAATTGTGGAAAACGTTCTACACCCATcc                                                                                                                   >  2:56845/1‑251 (MQ=255)                                                                                                                               aGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACGCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGaa                                                                                                               <  1:263659/251‑1 (MQ=255)                                                                                                                                  tAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGCGGGAATCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGATCGACACCCATCCGGGAggg                                                                                                            >  1:237157/1‑251 (MQ=255)                                                                                                                                    aGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGt                                                                                                                                                                                                                <  1:376270/149‑1 (MQ=255)                                                                                                                                    aGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGt                                                                                                                                                                                                                >  2:376270/1‑149 (MQ=255)                                                                                                                                    aGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGtt                                                                                                          <  1:61151/251‑1 (MQ=255)                                                                                                                                                            ttCTGTAATGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGGTCATCGTGATCCCTTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGACGCAGAATTGTGGAAAAAGTTCGATCCCCCTCCGGAAGGGTTAAATCAgatcgaaacgcacactgc                                                                                  >  2:355838/1‑234 (MQ=255)                                                                                                                                                             tCTGAACTGTGTCAGACAGTGCGGAGTGACTCCTTATGTTTAAAGAAATTTTTACCCGTCTCTTGCGCCATTGACCGTCCCGTCGGGTTCAGCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGGCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGTAAGGGTTAAATCAGGCGGAATTGGAATCTGc                                                                                  <  2:200491/250‑1 (MQ=255)                                                                                                                                                                gTACTGTGTCAGACAGTGCGGAGGGACTCCTTATGTTGAAAGAAATTGTTACCCGGCTCATTCGCCATTTACCTGCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGccc                                                                                <  2:319630/249‑1 (MQ=255)                                                                                                                                                                   cTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGa                                                                            >  1:161244/1‑250 (MQ=255)                                                                                                                                                                     gTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTTAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGaaca                                                                         <  1:297033/251‑1 (MQ=255)                                                                                                                                                                          aGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGGGGAAAACGTTCGACACCCCTCCGGAAGGGTTAAATCAGGCGGGAGTGGAATCTGCCCGCGAACAACATg                                                                    >  1:351206/1‑251 (MQ=255)                                                                                                                                                                          aGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGGGCGCATGGCTTGAAAAGGGCGGGGAGGCCCGAAGAAGAATTGTGGAAAACGTTCGAAACCCCTCCGGAGGGGTTAAATCAGGGGGAAGGGGAATCTGCCCCCGaacaaacgg                                                                    >  1:100622/1‑247 (MQ=255)                                                                                                                                                                             cAGTGCGGAGGGACTCCTTATGTTTAAATAAATTTTTACCCGTCTCATTCGCCATTTCCCTTCCCGTCGGGTTCATCCTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATTCCCAAAGAAGAATTGTGGAAAACGTTCGACGCCCATCCGGAAGAGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAGCAACATGGTg                                                                 <  1:376459/251‑1 (MQ=255)                                                                                                                                                                                   ggAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGGCCGAAGAAGAATTGGGGAAAACGTTCGACACCCCTCAGGAAGGGGTAAATCAGGCGGAAGTGGAATCTGCCCCCGCACAACATGGTGAaaat                                                            >  1:199177/1‑250 (MQ=255)                                                                                                                                                                                       ggACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGCAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGGAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAAtt                                                       >  1:428074/1‑251 (MQ=255)                                                                                                                                                                                                     tAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGGAGAAGAATTGTGGAAAACGTTCGACAACCCTCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCCCGGACAACATGGGGAAAATAAATTACCCGCACaacaac                                         >  2:48765/1‑251 (MQ=255)                                                                                                                                                                                                      aaaGAAATTTTTACCCGGCTCATTCGCCATTGACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAAcc                                        <  1:428494/251‑1 (MQ=255)                                                                                                                                                                                                       aaGAAATGTTTCCCCGGCTCATTCGCCATGTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCg                                       <  2:422725/251‑1 (MQ=255)                                                                                                                                                                                                                 tACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGGAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGGAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGGGAAAATAAATTACCCGCACAACAACCGTCGCCggggg                             >  1:345423/1‑247 (MQ=255)                                                                                                                                                                                                                         tCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGGACAt                     >  1:342857/1‑251 (MQ=255)                                                                                                                                                                                                                          cATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCGGTGGTGGGTACAtt                    <  1:343001/251‑1 (MQ=255)                                                                                                                                                                                                                            ttCGCCATTTACCTTCACGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGACAATGGCGGTGATGCCCGAAGAAGGATTGTGGAAAACGTTCGGCACCCATCCGGAAGGGTTAAATCCGGCGGAAAAGGAATCTGACCGCGGACAACAAGGGGAAAATAAATATCCCGCACAAAAACCGTCGCCGTGGGGGGTACAtttt                  >  1:185715/1‑250 (MQ=255)                                                                                                                                                                                                                            gtcggcATTTACCGTCCCGTCGGGTTCATCGTGCGCCATTGGCGGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGGGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGTGTTATATCTGGCGTAAGTGGAATCGGCCCGCGAACGACATGGTGAAAATAAATGACCCGCGCAACAACCTTCGGCGGAGGGGGTACAttt                   <  2:444774/245‑1 (MQ=255)                                                                                                                                                                                                                              cGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGGAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCCCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGGGGGGACATTTAtt                >  1:146433/1‑250 (MQ=255)                                                                                                                                                                                                                                    ttACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCt           <  2:342857/250‑1 (MQ=255)                                                                                                                                                                                                                                           cccGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGGAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGGAGGGGTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGGGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGGACATTTATGGGTCTGCTATCGc   >  1:154928/1‑251 (MQ=255)                                                                                                                                                                                                                                            ccGTCGGGTTCATCGTGTTCCATTGCCTGTCTCGCAGCAGACAGTGAATACGGTGGTCCCGCCGGCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGTGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTTTGGGTCTGCTATCGCa  <  2:23408/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      GAACCTGATCCCACGCCTCTCCCTCGACGGAGATTAAAACTTTTCCGGTAAGCCCGTCTTTTCACGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCA  >  NC_000913/4467434‑4467919

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