Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	22,105	T→G	G233G (GGT→GGG)	ribF →	bifunctional riboflavin kinase/FAD synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	22,105	0	T	G	100.0%	219.1 / NA	61	G233G (GGT→GGG)	ribF	bifunctional riboflavin kinase/FAD synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (34/27);  total (34/27)

ACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGC  >  NC_000913/21958‑22249                                                                                                                                                    |                                                                                                                                                 aCTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGa                                                                                                                                                <  1:122978/150‑1 (MQ=255) aCTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGa                                                                                                                                                <  1:31932/150‑1 (MQ=255)  cTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGaa                                                                                                                                               >  2:439363/1‑150 (MQ=255)   tGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGaaa                                                                                                                                              >  1:64035/1‑150 (MQ=255)    ggggCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAg                                                                                                                                            <  1:846788/151‑1 (MQ=255)         aCCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGtt                                                                                                                                       >  2:227003/1‑151 (MQ=255)               ttGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCg                                                                                                                                  <  1:174598/150‑1 (MQ=255)                        ccGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCaa                                                                                                                         >  2:623258/1‑150 (MQ=255)                               tgtAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCgg                                                                                                                  >  2:308664/1‑150 (MQ=255)                                gtAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGaa                                                                                                                <  1:682486/151‑1 (MQ=255)                                            gatgaATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACg                                                                                                     <  1:215356/150‑1 (MQ=255)                                                aaTTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGc                                                                                                >  2:582829/1‑151 (MQ=255)                                                     gggCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGt                                                                                            >  1:293149/1‑150 (MQ=255)                                                      ggCGCCCTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGt                                                                                            <  2:158315/149‑1 (MQ=255)                                                       gcgcACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTAtt                                                                                         <  1:675608/151‑1 (MQ=255)                                                       gcgcACTACAGGTTTCCCGACGGCGAATGTACAGCTGCGCCGTCAAGTTTCCCCGGAGAAAGGGGTTTATGCGGTAGAAGTGCTGGGGCTCGAGGAAAAGCCGTTACCCGGCGTGGCAAAAATCGGAACACGCCCAACGGTTGCCGGTAtt                                                                                         >  1:246911/1‑151 (MQ=255)                                                           aCTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGcc                                                                                     >  1:932366/1‑151 (MQ=255)                                                            cTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGcc                                                                                     >  1:739687/1‑150 (MQ=255)                                                              ataGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCcagc                                                                                  >  1:741313/1‑151 (MQ=255)                                                                    ttCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTgg                                                                            >  2:98542/1‑151 (MQ=255)                                                                     tCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTgg                                                                            >  2:537560/1‑150 (MQ=255)                                                                       ccGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGa                                                                           >  1:721952/1‑149 (MQ=255)                                                                       ccGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAg                                                                         <  1:367685/151‑1 (MQ=255)                                                                        cGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAg                                                                         >  2:572742/1‑150 (MQ=255)                                                                        cGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGt                                                                        >  2:361508/1‑151 (MQ=255)                                                                            ggCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCAt                                                                    >  2:931323/1‑151 (MQ=255)                                                                              cGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCAttt                                                                  <  1:62866/151‑1 (MQ=255)                                                                               gAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTg                                                                 >  2:751085/1‑151 (MQ=255)                                                                                     acagCTGCGCAGTTAGGTTTCCGAGGTGAAAGGGGATTATGCGGCAGAAGTGCTGGGACTCGGGGAAAAGCCGCTACCCGGCGTGGCAAACATCGGAACACGCCCAAGGGTTGCGGGTAATCGCCAGCAGCTGGAAGTGCATTGGtt                                                               >  1:698434/4‑147 (MQ=255)                                                                                     aCCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGAt                                                           <  2:690574/151‑1 (MQ=255)                                                                                       cGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGt                                                         >  1:762899/1‑151 (MQ=255)                                                                                       cGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGt                                                         >  2:347666/1‑151 (MQ=255)                                                                                           gcgcCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGc                                                      <  1:226998/150‑1 (MQ=255)                                                                                            cgcCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCaa                                                    <  1:623238/151‑1 (MQ=255)                                                                                              ccGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATg                                                  <  2:192997/151‑1 (MQ=255)                                                                                              ccGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATg                                                  <  1:236214/151‑1 (MQ=255)                                                                                                  cAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGAcc                                              >  2:896887/1‑151 (MQ=255)                                                                                                  cAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGAcc                                              <  2:802267/151‑1 (MQ=255)                                                                                                   aGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGAcc                                              <  1:81637/150‑1 (MQ=255)                                                                                                   aGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCt                                             >  1:767379/1‑151 (MQ=255)                                                                                                   aGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCt                                             >  1:309551/1‑151 (MQ=255)                                                                                                    ggTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCtt                                            >  1:167888/1‑151 (MQ=255)                                                                                                    ggTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCtt                                            <  1:577126/151‑1 (MQ=255)                                                                                                         ccccGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACgg                                       >  2:4963/1‑151 (MQ=255)                                                                                                           ccGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTc                                     <  2:439492/151‑1 (MQ=255)                                                                                                            cGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCg                                    <  1:837990/151‑1 (MQ=255)                                                                                                             ggTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGc                                   <  1:98754/151‑1 (MQ=255)                                                                                                                 aaaGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCata                               <  2:164512/151‑1 (MQ=255)                                                                                                                 aaaGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCat                                >  1:362984/1‑150 (MQ=255)                                                                                                                      ggTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACa                           >  2:153491/1‑150 (MQ=255)                                                                                                                           aTGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAagtagt                     >  1:216851/1‑151 (MQ=255)                                                                                                                           aTGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAagtagt                     <  1:760616/151‑1 (MQ=255)                                                                                                                            tGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGtg                    >  1:851089/1‑151 (MQ=255)                                                                                                                            tGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGtg                    >  2:910567/1‑151 (MQ=255)                                                                                                                               ggTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGtgctg                 >  1:503384/1‑151 (MQ=255)                                                                                                                                    aaGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTa             <  2:514999/150‑1 (MQ=255)                                                                                                                                      gTGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTaaaa          <  2:277412/151‑1 (MQ=255)                                                                                                                                       tGCTGGGCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTaaaa          >  2:25638/1‑150 (MQ=255)                                                                                                                                           gggCCTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTaaaaaaa       <  1:447642/149‑1 (MQ=255)                                                                                                                                              ccTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGc  <  2:23473/151‑1 (MQ=255)                                                                                                                                               cTCGGGGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGc  >  2:665959/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                 ACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGC  >  NC_000913/21958‑22249

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