Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,011,802	C→T	P243S (CCG→TCG)	bioP →	biotin transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,011,802	0	C	T	100.0%	208.8 / NA	60	P243S (CCG→TCG)	bioP	biotin transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/26);  total (34/26)

CCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTC  >  NC_000913/4011555‑4012045                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     ccTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCgg                                                                                                                                                                                                                                                   <  1:5436/251‑1 (MQ=255)     aTGGAATCTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTTGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAgg                                                                                                                                                                                                                                               <  2:245008/251‑1 (MQ=255)     aTGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAgg                                                                                                                                                                                                                                               >  1:242557/1‑251 (MQ=255)      tGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAggg                                                                                                                                                                                                                                              >  1:8211/1‑251 (MQ=255)        gAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGgct                                                                                                                                                                                                                                            <  1:167349/251‑1 (MQ=255)         aaaCTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGgctg                                                                                                                                                                                                                                           <  1:265916/251‑1 (MQ=255)                       ccACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGc                                                                                                                                                                                                                             >  1:5606/1‑251 (MQ=255)                         acaGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGGAAACCTGGCt                                                                                                                                                                                                                            >  2:76758/1‑250 (MQ=255)                                   gCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCAcc                                                                                                                                                                                                                  >  1:113397/1‑250 (MQ=255)                                    ccTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGAGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGCGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTAGGCAGGGCTGATGGTAAACCTGGCTATCTGGCACCa                                                                                                                                                                                                                 >  1:240131/1‑250 (MQ=255)                                     cTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAAc                                                                                                                                                                                                               >  1:145714/1‑251 (MQ=255)                                          cGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCg                                                                                                                                                                                                          >  1:302861/1‑251 (MQ=255)                                          cGTGGTTATATTTTGGCGTGTTTATGGTCGCAGTGATTGAATGGTTCCTGTTGGGACATGCGCAGAAAATGACGCAAACCACGATGCAATGGGGCATTATGGTGATTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGAGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATG‑ACAGTTTCGGCAGGGCTGCTGGGAAACCTTGCTATACGGCACCACCAGCCg                                                                                                                                                                                                          >  2:41987/1‑251 (MQ=255)                                              gTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACt                                                                                                                                                                                                      <  2:145714/251‑1 (MQ=255)                                              gTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACt                                                                                                                                                                                                      >  1:55559/1‑251 (MQ=255)                                                   aTCTTGGCGCGTTTATGGGGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACAACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCCTTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGcc                                                                                                                                                                                                  >  1:236308/1‑250 (MQ=255)                                                       tGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCa                                                                                                                                                                                                 <  2:111327/247‑1 (MQ=255)                                                               ttCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTAc                                                                                                                                                                                      <  1:316235/250‑1 (MQ=255)                                                                    gtggCAGTGATTGCATGGTTCTTGCTGGGAAATGCGAAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTGGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGgcgc                                                                                                                                                                                <  2:11434/251‑1 (MQ=255)                                                                       gCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCt                                                                                                                                                                             >  1:184261/1‑251 (MQ=255)                                                                           tGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTg                                                                                                                                                                         >  1:207285/1‑251 (MQ=255)                                                                           tGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATAATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTg                                                                                                                                                                         >  1:8222/1‑251 (MQ=255)                                                                           tGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCGGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACCGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTg                                                                                                                                                                         >  2:267962/1‑251 (MQ=255)                                                                            gATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGa                                                                                                                                                                        <  1:6207/251‑1 (MQ=255)                                                                                             tGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGACGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACTGCGCGACGCAGGTGGACGCCGGACCGCTGGGCATTATGAATTATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTGATTCCAGGGTCGCGGGTGAGCCTGGACTCACTGTgg                                                                                                                                                       >  2:303749/1‑251 (MQ=255)                                                                                                        cAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTa                                                                                                                                             >  1:242214/1‑250 (MQ=255)                                                                                                           aaaaTGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTg                                                                                                                                         >  1:67849/1‑251 (MQ=255)                                                                                                                   gCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTc                                                                                                                                 <  1:292226/251‑1 (MQ=255)                                                                                                                    cAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTcc                                                                                                                                <  1:183461/251‑1 (MQ=255)                                                                                                                     aaaCCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGAACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCcg                                                                                                                               >  2:28133/1‑251 (MQ=255)                                                                                                                       aCCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCcgcg                                                                                                                             >  2:276701/1‑251 (MQ=255)                                                                                                                                 aaTGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACa‑a                                                                                                                  <  2:146645/251‑1 (MQ=255)                                                                                                                                         aTTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGc                                                                                                             >  1:201899/1‑248 (MQ=255)                                                                                                                                         aTTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCCGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGc                                                                                                             >  2:33506/1‑248 (MQ=255)                                                                                                                                           tCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAAATACGGCGCGAAGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGGTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGGGGGAGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAgatg                                                                                                        >  1:324949/1‑251 (MQ=255)                                                                                                                                            cTGGTGTTTCTTGGCGTCGCGGCTTCCGGGATTCGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGACCGCTGGGCATTATGAATACTATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCCCTGCCCCACGATTATTACAGGCGCGCTGGTGATCCTGGCCCCACTGTGGGGGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AAAGGCAgatg                                                                                                        >  2:250178/1‑250 (MQ=255)                                                                                                                                               gtgtTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATct                                                                                                    <  2:207285/251‑2 (MQ=255)                                                                                                                                                    tCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGTAGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGCTGATCGCAGGc                                                                                               >  1:228073/1‑251 (MQ=255)                                                                                                                                                    tCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAgg                                                                                                <  1:76758/250‑1 (MQ=255)                                                                                                                                                       tGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTg                                                                                            >  1:101291/1‑251 (MQ=255)                                                                                                                                                                                 aCTTTATGTGGGACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCt                                                                   <  2:284991/250‑1 (MQ=255)                                                                                                                                                                                    ttATGTTGAGCTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAAGAGGCAC‑GTTTCGGCAGGGCTGCGGGTGAACCTGTCTATCTGGCACCAACAGCCGCTCTGGCCAACGTTTATTACAGGCGCGCGGGTGATCCTGGCCTCGCTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCGTCACA‑AACGGCAGATTGTCGGAGGTGTGATTGCGCGCTGGGCGAATAAACGCtt                                                                  <  2:332071/248‑1 (MQ=255)                                                                                                                                                                                    ttATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTcc                                                                <  1:318379/250‑1 (MQ=255)                                                                                                                                                                                      aTGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGt                                                              <  1:147516/250‑1 (MQ=255)                                                                                                                                                                                      aTGTGGAACTACGACCCGACGCAGGTGGCAGCCGGAATGCTGGGCATTATGTATACTCTGCAC‑GTTTCGACAGGGCAGCTGGTACACCTGGCTCTCTGACACCCACAGCCGCACTGGCCAACGTTTATTACAGGCGAGCTGGTGCACCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCCCTCTTC‑CATCAAGGCCCATGAACGCCGGCTTGTTTGCGCGCTGAGCGACTAAACGCTCCCga                                                              >  2:368596/1‑249 (MQ=255)                                                                                                                                                                                                aCGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCCC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTg                                                     <  1:36597/249‑1 (MQ=255)                                                                                                                                                                                                   gcgcGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACgct                                                <  2:55559/251‑1 (MQ=255)                                                                                                                                                                                                     gcgACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGCTGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGCCgctgc                                              >  2:36724/1‑251 (MQ=255)                                                                                                                                                                                                     gcgACGCAGGTGGACGCCGGAACGCTCGGCCTTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCACCAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑CACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACgctgc                                              >  2:169539/1‑251 (MQ=255)                                                                                                                                                                                                      cgACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACgctgct                                             <  2:5606/251‑1 (MQ=255)                                                                                                                                                                                                        aCGCAGGTGGACGCCGGACCGCTGGGCATTAGGGATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCg                                           <  2:83196/251‑1 (MQ=255)                                                                                                                                                                                                                    gCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCAcg                               <  1:28133/251‑1 (MQ=255)                                                                                                                                                                                                                       ggAACGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCAGCGCGCAcggcg                            <  1:292679/251‑1 (MQ=255)                                                                                                                                                                                                                         aaCGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACggcggc                          >  2:314973/1‑251 (MQ=255)                                                                                                                                                                                                                          aCGCTGGGCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACATGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGgcggcg                         >  2:134990/1‑251 (MQ=255)                                                                                                                                                                                                                                ggCATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGAACAGt                   >  1:201031/1‑251 (MQ=255)                                                                                                                                                                                                                                  cATTATGAATAATATGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCg                 >  1:285715/1‑251 (MQ=255)                                                                                                                                                                                                                                               atGCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACCGTCGGCTCCACaaaccc    >  2:274619/1‑247 (MQ=255)                                                                                                                                                                                                                                                 gCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCCCCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTc  <  2:133450/251‑1 (MQ=255)                                                                                                                                                                                                                                                 gCAC‑GTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTc  <  2:8211/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     CCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCAC‑GTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACA‑AACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTC  >  NC_000913/4011555‑4012045

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