Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	819,146	G→A	Q33Q (CAG→CAA)	moaD →	molybdopterin synthase sulfur carrier subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	819,146	0	G	A	100.0%	125.7 / NA	40	Q33Q (CAG→CAA)	moaD	molybdopterin synthase sulfur carrier subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/19);  total (21/19)

AATGGAAGCATTAACCGCGGCCTCCGTGGCGGCGCTGACCATTTATGACATGTGCAAAGCGGTGCAAAAAGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGGAGAAGAGTACCCGTGGCTGGCGGAGCGTGACGAAGACGGTGCGGTAGTC  >  NC_000913/818905‑819390                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      aaTGGAAGCATTAACCGCGGCCTCCGTGGCGGCGCTGACCATTTATGACATGTGCAAAGCGGTGCAAAAAGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCt                                                                                                                                                                                                                                              <  2:352023/251‑1 (MQ=255)  aTGGAAGCATTAACCGCGGCCTCCGTGGCGGCGCTGACCATTTATGACATGTGCAAAGCGGTGCAAAAAGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTg                                                                                                                                                                                                                                             >  1:197524/1‑251 (MQ=255)         cATTAACCGCGGCCTCCGTGGCGGCGCTGACCATTTATGACATGTGCAAAGCGGTGCAAAAAGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAATTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCag                                                                                                                                                                                                                                        >  1:200174/1‑249 (MQ=255)               ccGCGGCCTCCGTGGCGGCGCTGACCATTTATGACATGTGCAAAGCGGTGCAAAAAGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCg                                                                                                                                                                                                                                <  1:241951/251‑1 (MQ=255)                 gcgGGCTCCGTGGCGGCGCTGACCATTTATGACATGTGCAAAGCGGTGCAAAAAGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCt                                                                                                                                                                                                                              <  2:120249/251‑1 (MQ=255)                                               gACATGTGCAAAGCGGTGCAAAAAGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACGTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGGTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTa                                                                                                                                                                                                 <  2:301978/250‑1 (MQ=255)                                                        aaaGCGGTGCAAAAAGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCg                                                                                                                                                                                      >  2:308874/1‑251 (MQ=255)                                                           gCGGTGCAAAAAGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCa                                                                                                                                                                                   <  2:184555/251‑1 (MQ=255)                                                                  aaaaaGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGAc                                                                                                                                                                            >  2:114471/1‑251 (MQ=255)                                                                   aaaaGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACg                                                                                                                                                                           >  2:144299/1‑251 (MQ=255)                                                                    aaaGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGc                                                                                                                                                                          >  2:82778/1‑251 (MQ=255)                                                                                      ccGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATc                                                                                                                                                        >  1:65092/1‑251 (MQ=255)                                                                                           aCGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGctg                                                                                                                                                   >  2:113065/1‑251 (MQ=255)                                                                                           aCGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGct                                                                                                                                                    <  1:197581/250‑1 (MQ=255)                                                                                             gTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGACGCGTTACGCCAACACATGGCTGCGCAGAGCGATAGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGGGAGTTTTGACCATCCGctgac                                                                                                                                                 >  2:288501/1‑251 (MQ=255)                                                                                                 gCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCtgactgac                                                                                                                                             <  1:220082/251‑1 (MQ=255)                                                                                                     gCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGcg                                                                                                                                         <  2:330132/251‑1 (MQ=255)                                                                                                     gCGAAGAGCGGCGGCAAGTCGGGAGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCAAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCCACCAGACGCTGGGGAGTTTTGACCATCCGCTGACTGACGGcg                                                                                                                                         >  2:228048/1‑251 (MQ=255)                                                                                                        aaGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGcgacg                                                                                                                                      <  1:250146/251‑1 (MQ=255)                                                                                                                 ggcAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCtt                                                                                                                             <  2:26987/251‑1 (MQ=255)                                                                                                                           gTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCccgcc                                                                                                                   <  2:102499/251‑1 (MQ=255)                                                                                                                                               gattatTAATGTTTTTTTTTGCGCCCGGGGGCGCGTGTTGGGGGGAACAGATGAAACCGAAGGGGTGGCGGATTTCCCAAG‑TGTTGAGGCGTTACGCAAACCAAGGGCGGCGCGGAGGGATCGCTGGGCGCTGGCGCTGGAGGATGGCAAATTAG‑TGGGTGGCGTCACCCGGACTCTGGTGAGTTTTGGCCTCCCGGTGACTGAGGGCGACGTAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATg                                                                                               <  2:357125/247‑1 (MQ=255)                                                                                                                                                atgatTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATgg                                                                                              >  1:18226/1‑251 (MQ=255)                                                                                                                                                  gatTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGcc                                                                                            >  2:106139/1‑250 (MQ=255)                                                                                                                                                        aGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAAcc                                                                                      >  1:336606/1‑251 (MQ=255)                                                                                                                                                         gtttttttttgCGCCGAGGGGCGCGAGTTGGGGGGAACAGATGCAACCGAGGGGGCGGCGGTTTCCCCAAC‑TGTGGAAGCGTTCCGCCAACACATGGCTGCGCGGAGCGATCGCGGGGCGCTGGCACTGAAAGAGGGAAAATTAC‑TGGCTGCCGTCAACCAGCCGCTGGTGATTTTGGACCACCCGGTGCCTGACGGGGACGAAGTAGTTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAAcc                                                                                      <  2:255935/246‑1 (MQ=255)                                                                                                                                                               ttttcCGCCCAGGTGCGGGAGTTGGTGGGAACAGAGGCAACCGAAGTGGCGGGGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCGGGGCGCGGGCGCTGGAAGATGGCAAATTAC‑TGGCGCCCGTCACCCAGACGCTGGTGAGTTTTGACCACCCGCTGATGGACGGCGACGAAGTAGTTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAAtt                                                                                <  2:282976/250‑1 (MQ=255)                                                                                                                                                                   tCGCCCGGGTGCGCGAGTGGGGGGAAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTGGAGGCGTTACGCCAACACAGGGCGGCGCAGAGCGATCGCTGGGCGCGGGCGCGGGAAGATGGAAAATTAC‑GGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGCCCATCCGCTGACTGACGGCGACAAAGAAGTTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAAttgttg                                                                            <  2:287502/250‑1 (MQ=255)                                                                                                                                                                   tCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGCCCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCACAAttgttgt                                                                           >  2:241578/1‑251 (MQ=255)                                                                                                                                                                    cGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAAttgttgtt                                                                          >  1:183474/1‑251 (MQ=255)                                                                                                                                                                      cccAGGTGCGCGAGTTGGGGGGAACAGATGCAACCGAAGTGGCTGCGGTTTTCCCAAC‑GGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCGGGGCGCGGGCGCGGGAAGATGGAAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGAGGGCGACGAAGTAGTTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTgg                                                                        <  1:364231/251‑1 (MQ=255)                                                                                                                                                                           gcgcgcGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGc                                                                   >  1:45965/3‑251 (MQ=255)                                                                                                                                                                                 gAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGTTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGt                                                             >  2:264128/1‑251 (MQ=255)                                                                                                                                                                                      ggggggAACAGAGGCAACGGAAGGGGCGGCGGATTCCCCAAC‑TGTGGAGGCGTTCCGCCAACAAAGGGCGGCGCAGAGCGACCGTGGGGCGCTGGCGCGGGAAGAGGGAAAATTAC‑TGGCTGCCGTCAACCAGGCGCTGGTGAGTGTGGACCATCCGCGGACGGACGGCGAGGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGc                                                        <  2:319488/248‑1 (MQ=255)                                                                                                                                                                                               aGATGCAACCGCAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGgaga                                                >  2:378250/1‑250 (MQ=255)                                                                                                                                                                                                       aacgaagtgggtgagaattCCCAACCTGGTGAAGGGTTACGACAACACAGGGTTGCGCGGACCGTTCGCTGGGCGCTGGCGTGGGAAGAGGGCAA‑TTACCTGGCTACCGTCAACCAGACGCGGGTGATTTTTGCCCATCCGCTGACTGACGGCGACGAAGAAGCTTTCTCCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGTAGAAGAGTAcc                                        <  2:49465/234‑1 (MQ=255)                                                                                                                                                                                                               gCTGCGGATTTCCCAAT‑TGTTGAGGCGTTACCCAAACACATGGCGGCGCAGAGCGTTCGCGGGGCGCGGGCGCTGAAAGATGGCAATTTCC‑TGGCTGCCGTCAACCAGACGCGGGTGATTTTTGACCTCCCGCTGACTGACGGCGAAGAAGTAGTTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGGAGAAGAGTACCCGtggctgg                               <  2:197524/251‑1 (MQ=255)                                                                                                                                                                                                                 tGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGGAGAAGAGTACCCGTGGCTGGCg                             >  1:23650/1‑251 (MQ=255)                                                                                                                                                                                                                   cGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGGAGAAGAGTACCCGTGGCTGGCggg                           >  2:118212/1‑250 (MQ=255)                                                                                                                                                                                                                                              cGCCAACACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGGAGAAGAGTACCCGTGGCTGGCGGAGCGTGACGAAGACGGTGCGGTAGTc  >  1:261176/1‑250 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      AATGGAAGCATTAACCGCGGCCTCCGTGGCGGCGCTGACCATTTATGACATGTGCAAAGCGGTGCAAAAAGATATGGTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAAC‑TGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTAC‑TGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGGAGAAGAGTACCCGTGGCTGGCGGAGCGTGACGAAGACGGTGCGGTAGTC  >  NC_000913/818905‑819390

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