Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	796,803	G→A	D84N (GAT→AAT)	modC →	molybdate ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	796,803	0	G	A	100.0%	141.3 / NA	40	D84N (GAT→AAT)	modC	molybdate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/19);  total (21/19)

GGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGACAGCACCGGAATTGCTGTTGCTGGATGAACCGCTGGCGTCACTGGATATTCCGCGTAAACG  >  NC_000913/796559‑797047                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     ggAACTAAATTTTTCCCAGACGTGGGGCAACCATTGCCTGACTATAAAAGAACCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATgcgc                                                                                                                                                                                                                                                <  1:12969/251‑1 (MQ=255)  aatcTGAATTTTTCCCAGCCGTTGGCCACCAATTGCCTGACTATAAATGAACCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGTTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGcgcg                                                                                                                                                                                                                                               <  1:209269/248‑1 (MQ=255)        aaTTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGtt                                                                                                                                                                                                                                         >  1:247084/1‑251 (MQ=255)          gttttACCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGTCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCTTCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTACATGATGCCGAAACAGGTAGCAGCCTGAGGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGAGCGGCTGTTcc                                                                                                                                                                                                                                       >  2:245502/2‑251 (MQ=255)                                    gCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTg                                                                                                                                                                                                             >  1:57800/1‑251 (MQ=255)                                      cTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTgcg                                                                                                                                                                                                           >  1:360266/1‑251 (MQ=255)                                      cTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTgcg                                                                                                                                                                                                           >  1:31723/1‑251 (MQ=255)                                                                cAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAACCTTCGCTGATTAACGCCATCAGTGAACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAACGTGCGTGGCAATCTGCGCTACGGCATGTCGAACAGTATGGTCg                                                                                                                                                                                 >  1:232279/1‑251 (MQ=255)                                                                 aaTGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGa                                                                                                                                                                                <  2:93349/251‑1 (MQ=255)                                                                  aTGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGAt                                                                                                                                                                               >  1:175484/1‑251 (MQ=255)                                                                    ggCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCa                                                                                                                                                                             <  1:342643/251‑1 (MQ=255)                                                                    ggCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCa                                                                                                                                                                             <  1:342669/251‑1 (MQ=255)                                                                         cACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTc                                                                                                                                                                         <  1:174614/250‑1 (MQ=255)                                                                             gCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGGAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATa                                                                                                                                                                     >  1:322217/1‑250 (MQ=255)                                                                               tATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGc                                                                                                                                                                  >  2:327733/1‑251 (MQ=255)                                                                                         gTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGCAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGGCGAAAAGGATGGTCGATCAGTTCGATAAGCTGGTGGCGCt                                                                                                                                                        >  1:46550/1‑251 (MQ=255)                                                                                           ctcCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGGCGAAAAGCATGGTCGATCAGTTCGATAAGCTGGTGGCGCttt                                                                                                                                                      >  1:303318/1‑251 (MQ=255)                                                                                                   ccGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCAtt                                                                                                                                              <  1:290920/251‑1 (MQ=255)                                                                                                    cGGAAAAACTTCGCTGCTTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGCTTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTg                                                                                                                                             >  2:177425/1‑251 (MQ=255)                                                                                                         aaaCTTCGCTGATTAACGCCATCAGTGGACTGCCGCGACCGCAACAAGGGCGGATTGCCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGGGTTGGCTATGTTTTTCAGACTGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGACGCTCAGTTCGATAAGCTGGTGGCGATTTTAGGCATTGAACCg                                                                                                                                        >  1:221663/1‑251 (MQ=255)                                                                                                              tCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTg                                                                                                                                     <  2:264846/249‑1 (MQ=255)                                                                                                                gCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTg                                                                                                                                 <  1:126193/251‑1 (MQ=255)                                                                                                                   gATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCGGGTGGCGCTTTTAGGCAGTGAACCGTTGCTTGAcc                                                                                                                              <  2:247084/251‑1 (MQ=255)                                                                                                                                     aCTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGt                                                                                                            <  2:262387/251‑1 (MQ=255)                                                                                                                                         aCGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCgg                                                                                                        <  1:153033/251‑1 (MQ=255)                                                                                                                                                 gCAAAAAGGGCGGCTTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCCGCCTGTCCGGAGGCGaa                                                                                                 >  2:210741/1‑250 (MQ=255)                                                                                                                                                  cAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGaaaa                                                                                               <  2:345281/251‑1 (MQ=255)                                                                                                                                                           cGGATTGTCCTCAATGGGCGGGTCCTAAATGATGCCGCAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGAGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCCTGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGCCCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGt                                                                                      >  2:24057/1‑251 (MQ=255)                                                                                                                                                            ggATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGGTGATTGACCGTTTACCAGGCATCCTGTCCGGAGGCGAAAAACAGCGCGTg                                                                                     <  1:147183/251‑1 (MQ=255)                                                                                                                                                             gATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTgg                                                                                    <  2:196807/251‑1 (MQ=255)                                                                                                                                                                 gTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACCGCGCGTGGCGAt                                                                                >  1:72982/1‑251 (MQ=255)                                                                                                                                                                    cTCAATGGGCGGGTCCTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTgg                                                                             >  1:93926/1‑251 (MQ=255)                                                                                                                                                                         ggggCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCggg                                                                        <  2:88595/250‑1 (MQ=255)                                                                                                                                                                              gggTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTg                                                                   >  1:277489/1‑251 (MQ=255)                                                                                                                                                                                  aCTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTg                                                                >  2:60745/1‑250 (MQ=255)                                                                                                                                                                                      aaTGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCTTTTTCCAGGCAGCCTGTCCGGAGGCGAAAAACAGCTCGTGGCGATTGGTCGGGCTTTGCTGACAGc                                                           <  1:162943/251‑1 (MQ=255)                                                                                                                                                                                       atgatgCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGACAGc                                                           >  1:22182/1‑250 (MQ=255)                                                                                                                                                                                                                             cGTCGCGTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGACAGCACCGGAATTGCTGTTGCTGGATGAACCGCTGGCGTCACt                    <  2:92464/251‑1 (MQ=255)                                                                                                                                                                                                                                cgcgTTGGCTATGTTTTTCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGACAGCACCGGAATTGCTGTTGCTGGATGAACCGCTGGCGTCACTgg                  <  1:84581/250‑1 (MQ=255)                                                                                                                                                                                                                                               tttCAGAATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGACAGCACCGGAATTGCTGTTGCTGGATGAACCGCTGGCGTCACTGGATATTCCGCGTaaatg  >  2:119247/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     GGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGACAGCACCGGAATTGCTGTTGCTGGATGAACCGCTGGCGTCACTGGATATTCCGCGTAAACG  >  NC_000913/796559‑797047

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