Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,146,847	C→T	D232N (GAT→AAT)	alkA ←	DNA‑3‑methyladenine glycosylase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,146,847	0	C	T	100.0%	153.1 / NA	43	D232N (GAT→AAT)	alkA	DNA‑3‑methyladenine glycosylase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/22);  total (21/22)

GATGCGACGCTAACGCGTCTTATCATGCCTACAAATCGCTCATTCCCCAGGCCGGATAAGGCGCTCGCACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATCTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCG  >  NC_000913/2146604‑2147092                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      gATGCGACGCTAACGCGTCTTATCATGCCTACAAATCGCTCATTCCCCAGGCCGGATAAGGCGCTCGCACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCg                                                                                                                                                                                                                                                  <  1:447790/249‑1 (MQ=255)        cGCTAACGCGTCTTATCATGCCTACAAATCGCTCATTCCCCAGGCCGGATAAGGCGCTCGCACCGCATCCGGCGACCAACGTCACGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATCATCATCCGGCAGGAAAACATTTTTCGCCTgccagc                                                                                                                                                                                                                                         >  1:215776/1‑251 (MQ=255)                  tCTTATCATGCCTACAAATCGCTCATTCCCCAGGCCGGATAAGGCGCTCGCACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAAAATTTTTCGCCTGCCAGCCACgcaaagc                                                                                                                                                                                                                               >  1:222582/1‑251 (MQ=255)                            ccTACAAATCGCTCATTCCCCAGGCCGGATAAGGCGCTCGCACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTc                                                                                                                                                                                                                     <  1:454451/251‑1 (MQ=255)                                  aaTCGCTCATCCCCCTGACCGGATAAGGCACTCCCCCCGCATCCGGCGACCAACGTCATGCTTCGTCTTGTTGCCACCCTGCCGTGTTCCAGATATGCCACAGCGCATAAGAACGCCAGGTCTTCCAGCGCTCGGCATAACGGCGGATTTGCTCCGGTGTCATTCCCGCAAATCGCTGTTTCATCCTATAATCACCCGGCAGGCAAACATTTTTCGCCTGCCATCCGCGCAAAGCAAAATAATTCGCCGTc                                                                                                                                                                                                               <  2:101670/251‑1 (MQ=255)                                   aTCGCTCATTCCCCAGGCCGGATAAGGCGCTCGCACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTcc                                                                                                                                                                                                              <  1:427033/251‑1 (MQ=255)                                      gCTCATTCCCCAGGCCGGATAAGGCGCTCGCACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAgc                                                                                                                                                                                                           <  2:43604/251‑1 (MQ=255)                                                 aGGCCGGATAAGGCGCTCGCACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATAccc                                                                                                                                                                                                <  1:405592/251‑1 (MQ=255)                                                 aGGCCGGATAAGGCGCTCGCACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATAccc                                                                                                                                                                                                <  1:353007/251‑1 (MQ=255)                                                  ggCCGGATAAGGCGCTCGCACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCCGCCACGCAAAACCAAATAATTCGCCGTCCAGCGCCCGATACCCg                                                                                                                                                                                               >  1:39282/1‑251 (MQ=255)                                                          aaGGCGCTCGCACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGtt                                                                                                                                                                                       >  2:55930/1‑251 (MQ=255)                                                           aGGCGCTCGCACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGttt                                                                                                                                                                                      <  2:54796/251‑1 (MQ=255)                                                                    cACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGttt                                                                                                                                                                             <  1:390940/251‑1 (MQ=255)                                                                        gCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCAt                                                                                                                                                                         >  2:241931/1‑251 (MQ=255)                                                                          aTCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCACCAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCg                                                                                                                                                                       <  1:127732/251‑1 (MQ=255)                                                                                 gACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGc                                                                                                                                                                                                                                                 <  2:331346/170‑1 (MQ=255)                                                                                 gACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGc                                                                                                                                                                                                                                                 >  1:331346/1‑170 (MQ=255)                                                                                           aTGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGGGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGccc                                                                                                                                                      >  1:266391/1‑251 (MQ=255)                                                                                             gCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCg                                                                                                                                                     <  1:266175/250‑1 (MQ=255)                                                                                                                       ataCCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCt                                                                                                                                                                    >  1:331627/1‑209 (MQ=255)                                                                                                                       ataCCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCt                                                                                                                                                                    <  2:331627/209‑1 (MQ=255)                                                                                                                                                       aGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATTATCCGGCAGGAAAACATTTTTTGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGCTACCCGGAAAAGGTTGCAGCGGTTTCATCGCCTGCTCCACATCGCCCGGTATTTTCATTGGTAAGGGGCCCTCCGGCGCCGCACTTGCCAGATGAATCAGCGCctct                                                                                          >  2:230013/1‑251 (MQ=255)                                                                                                                                                       aGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGCGTCATTCCCGGCAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCCGCTCCACATCGCCCGGCATTGTCATTGGTGAGGTGCCCTCCAGCGCCGCATTTGCCAGCAGAATCAGAGCctcc                                                                                          >  1:192436/1‑250 (MQ=255)                                                                                                                                                       aGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGAGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACAACGCCCGGGATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTTCCAGATGAATCAGCGCctct                                                                                          >  1:391380/1‑251 (MQ=255)                                                                                                                                                        gggCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTg                                                                                         <  2:428964/251‑1 (MQ=255)                                                                                                                                                        gggCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAAACGATGTTTAATCAGATAATCATCCGGCAGGCAAACATTTTTCGCCTGACAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGAAACCCGGAAAAGTTTGCAGCGTTTTCACCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGAAGAATCAGCGCCTCTg                                                                                         >  1:64195/1‑251 (MQ=255)                                                                                                                                                                   ccTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAAc                                                                              <  2:213506/250‑1 (MQ=255)                                                                                                                                                                      cGGCATAACGGCGGATTTGCGCCGGGGTCATTCCCGGAAATCGCTGTTTAATCCGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGAAACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCCTTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGc                                                                           >  1:232031/1‑251 (MQ=255)                                                                                                                                                                         cATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCAt                                                                         >  1:173345/1‑250 (MQ=255)                                                                                                                                                                                gcggATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACg                                                                 <  1:386757/251‑1 (MQ=255)                                                                                                                                                                                gcggATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACg                                                                 >  1:339595/1‑251 (MQ=255)                                                                                                                                                                                           ccGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGcc                                                      <  1:230013/251‑1 (MQ=255)                                                                                                                                                                                                  aaTTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATCCCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCTGGt                                               <  2:401349/250‑1 (MQ=255)                                                                                                                                                                                                        cGGAAATCGCTGTTTAATCAGATAATCATCCGGAAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGAATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTctgctg                                         >  1:416830/1‑251 (MQ=255)                                                                                                                                                                                                          gAAATCGCTGTATAATCAGATAATCATCCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTctgctgct                                       >  2:175155/1‑251 (MQ=255)                                                                                                                                                                                                                   gtttagtaaggtAATCATCCGACAGGAAAACATTTTTAGCCTGCCAGACACGCAAAGCAAAATAAGTCGCCGTGAAGTGCCCGATAACCGGAAACATTTGCAGCGTTTTCATCGACTGCTCCACAGCGCCCGGTATTGTCATTGGTAAGGTGCCATCCAGCGCCGCATTTGCAAAATGAATCAGCGCCTCTGCCCGTTTCAGCGGCATATCTTACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGc                              >  2:296048/12‑251 (MQ=255)                                                                                                                                                                                                                        aTCAGATAATCATCCGGCAGGAAAACATTTCTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACAGCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCGTTTGCCAGATGAATCAGCGCCTCTGCCCTTTTCAACGGCATTCCTAACGCTTTTAATGCCTGCGGTTCTGCTGCTGCCATCCGCTGAgg                         <  2:154511/251‑1 (MQ=255)                                                                                                                                                                                                                               aaTCATCCGACAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAAATCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCAGTGGTAAGGTGCCTTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTCTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGATGCGTCggg                  <  2:215776/251‑1 (MQ=255)                                                                                                                                                                                                                                   atcCGGCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGTTGCTGCCAGCCGCTGAGGCGTCGGGAAGc              >  2:233147/1‑251 (MQ=255)                                                                                                                                                                                                                                    cccGGCAGGAAAACATTGTTCGCCTGCCACCCACGCAAAGCAAAATAATTCGCCGTCCTGCGCCCGGTACCCGGAAAAGTTTGCATCGTTTTCATCGCCTGCTCCGCATCGCCCGGTAGTGTCACTGGTAAGGTGCCCTCCATCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCGATCGCTTGTAATGCCTGCGGGTCTGCTGCTGCCAGGCGCGGTGGCGTCGGGAAGCa             <  2:192436/250‑1 (MQ=255)                                                                                                                                                                                                                                       ggCAGGAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGc              >  1:106645/1‑247 (MQ=255)                                                                                                                                                                                                                                           ggAAAACATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCTATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCATCCGCTGAGGCGTCGGGAAGCAGATATAc      <  2:222582/251‑1 (MQ=255)                                                                                                                                                                                                                                               aaCATTTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCATATATACTCCg  >  1:330616/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      GATGCGACGCTAACGCGTCTTATCATGCCTACAAATCGCTCATTCCCCAGGCCGGATAAGGCGCTCGCACCGCATCCGGCGACCAACGTCATGCTTCGTCTGGTTGCCAGCCTTCCGTATACCAGATATGCAACAGCGCATAAGAACGCCAGGGCTTCCAGCGCTCGGCATAACGGCGGATTTGCGCCGGTGTCATTCCCGGAAATCGCTGTTTAATCAGATAATCATCCGGCAGGAAAACATCTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCG  >  NC_000913/2146604‑2147092

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