Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,191,997	G→C	P213P (CCC→CCG)	thiG ←	thiamine biosynthesis ThiGH complex subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,191,997	0	G	C	100.0%	281.3 / NA	77	P213P (CCC→CCG)	thiG	thiamine biosynthesis ThiGH complex subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (30/47);  total (30/47)

CGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGG  >  NC_000913/4191849‑4192142                                                                                                                                                     |                                                                                                                                                  cGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACcgg                                                                                                                                                 <  2:1446461/151‑1 (MQ=255)  gATCGCTGAGGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACcgg                                                                                                                                                 <  2:1428368/150‑1 (MQ=255)   aTCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGAt                                                                                                                                               <  2:1306293/151‑1 (MQ=255)   aTCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGAt                                                                                                                                               <  1:978683/151‑1 (MQ=255)     cGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCg                                                                                                                                             >  2:626550/1‑151 (MQ=255)     cGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAACGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCg                                                                                                                                             >  2:576028/1‑151 (MQ=255)         gAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCcg                                                                                                                                         >  1:1403742/1‑151 (MQ=255)          aaGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCcg                                                                                                                                         <  2:1137412/150‑1 (MQ=255)             gTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGa                                                                                                                                      >  1:396516/1‑150 (MQ=255)             gTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGa                                                                                                                                      >  1:25967/1‑150 (MQ=255)             gTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGa                                                                                                                                      >  2:1556433/1‑150 (MQ=255)              ttttCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACg                                                                                                                                    >  2:614501/1‑151 (MQ=255)                             cTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTAtt                                                                                                                      >  1:1031325/1‑150 (MQ=255)                             cTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTc                                                                                                                     <  2:249263/151‑1 (MQ=255)                               ccAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCAc                                                                                                                   <  1:31736/151‑1 (MQ=255)                               ccAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCAc                                                                                                                   <  1:318620/151‑1 (MQ=255)                                   aaaTCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAac                                                                                                               <  1:296894/151‑1 (MQ=255)                                   aaaTCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAac                                                                                                               <  1:674552/151‑1 (MQ=255)                                     aTCCGGTCAGGGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAacac                                                                                                             <  2:494682/151‑1 (MQ=255)                                     aTCCGGTCAGCGGGCTGGTTGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCATTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAacac                                                                                                             <  2:1461266/151‑1 (MQ=255)                                     aTCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAacac                                                                                                             <  1:705687/151‑1 (MQ=255)                                     aTCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAacac                                                                                                             <  2:1396762/151‑1 (MQ=255)                                     aTCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAacac                                                                                                             <  2:1151476/151‑1 (MQ=255)                                     aTCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAacac                                                                                                             <  2:1217476/151‑1 (MQ=255)                                     aTCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAacac                                                                                                             <  2:723612/151‑1 (MQ=255)                                     aTCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAacac                                                                                                             <  2:865938/151‑1 (MQ=255)                                          gTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGt                                                                                                        >  2:1543886/1‑151 (MQ=255)                                                gggCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGcc                                                                                                   >  1:98777/1‑150 (MQ=255)                                                         gCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTcc                                                                                          <  2:1155276/150‑1 (MQ=255)                                                         gCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTcc                                                                                          <  1:1570541/150‑1 (MQ=255)                                                                  aaaTAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTTTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTgcg                                                                                >  1:1178722/1‑151 (MQ=255)                                                                     tAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCgg                                                                             >  1:217774/1‑151 (MQ=255)                                                                     tAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCgg                                                                             >  1:1236253/1‑151 (MQ=255)                                                                     tAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCgg                                                                             <  1:118803/151‑1 (MQ=255)                                                                       aCTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGGCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCa                                                                           <  1:620712/151‑1 (MQ=255)                                                                                gCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGa                                                                  >  2:1596156/1‑151 (MQ=255)                                                                                gCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGa                                                                  >  1:745524/1‑151 (MQ=255)                                                                                 cccGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGa                                                                  <  2:1481876/150‑1 (MQ=255)                                                                                     gTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGc                                                              >  2:1094705/1‑150 (MQ=255)                                                                                     gTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGc                                                              >  1:889700/1‑150 (MQ=255)                                                                                     gTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGc                                                              >  2:1341789/1‑150 (MQ=255)                                                                                     gTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGc                                                              >  1:79101/1‑150 (MQ=255)                                                                                         ggACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATg                                                         >  1:17165/1‑151 (MQ=255)                                                                                                     cAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGTGCGGCACGGCCGGGAACGCCGATGCCAGCATCGACa                                             >  2:1220371/1‑151 (MQ=255)                                                                                                        tAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAacc                                          <  1:1623175/151‑1 (MQ=255)                                                                                                        tAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAac                                           <  2:746547/150‑1 (MQ=255)                                                                                                          ggCCTGCTTCTACCGCCAGACGAACATCCTTCGCCATGTTGACCGGGTCGTCCGCGACGGCACTCGCCGCATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCCGCATGGCTGGGAAGCCCGAGGCCAGCTTCGACTacca                                         >  2:1193301/1‑150 (MQ=255)                                                                                                           gCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAaccac                                        <  2:977319/150‑1 (MQ=255)                                                                                                            ccTGCTTCTCCCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGGCGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAaccacc                                       <  2:848272/150‑1 (MQ=255)                                                                                                            ccTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAaccacc                                       <  1:1442202/150‑1 (MQ=255)                                                                                                              tGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGc                                    <  1:227894/151‑1 (MQ=255)                                                                                                                cTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCa                                   <  1:860245/150‑1 (MQ=255)                                                                                                                cTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCa                                   >  1:830633/1‑150 (MQ=255)                                                                                                                cTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCAc                                  <  1:1016079/151‑1 (MQ=255)                                                                                                                cTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCAc                                  <  2:1233534/151‑1 (MQ=255)                                                                                                                cTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCAc                                  <  1:1474973/151‑1 (MQ=255)                                                                                                                cTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCAc                                  <  1:375394/151‑1 (MQ=255)                                                                                                                     aCCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTgg                             <  2:429599/151‑1 (MQ=255)                                                                                                                       cGCCAGACGAAATGCCTGCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGcc                           <  2:19502/151‑1 (MQ=255)                                                                                                                            gACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCctg                         >  2:804429/1‑148 (MQ=255)                                                                                                                            gACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCctg                         >  2:552956/1‑148 (MQ=255)                                                                                                                             actaaaTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCctgctg                      <  2:896521/147‑1 (MQ=255)                                                                                                                             aCGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCctgctg                      <  1:290197/150‑1 (MQ=255)                                                                                                                                aaaTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGat                   <  2:155883/150‑1 (MQ=255)                                                                                                                                aaaTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGat                   <  1:485282/150‑1 (MQ=255)                                                                                                                                aaaTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGat                   <  1:576252/150‑1 (MQ=255)                                                                                                                                     ccTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGAtaataa              <  2:1289866/150‑1 (MQ=255)                                                                                                                                     ccTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATaataat             >  2:134635/1‑151 (MQ=255)                                                                                                                                      cTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAAtc            <  1:714396/151‑1 (MQ=255)                                                                                                                                       ttCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAAtct           <  1:1212370/151‑1 (MQ=255)                                                                                                                                       ttCGCCATGTGGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTTCTGGATAATAAtct           <  1:1484475/151‑1 (MQ=255)                                                                                                                                        tCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAAtct           >  1:1615698/1‑150 (MQ=255)                                                                                                                                            cATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGc      <  1:1180276/151‑1 (MQ=255)                                                                                                                                             aTGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCa     >  2:1056411/1‑151 (MQ=255)                                                                                                                                              tGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCAt    <  1:565207/151‑1 (MQ=255)                                                                                                                                              tGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCAt    <  1:1305685/151‑1 (MQ=255)                                                                                                                                              tGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCAt    <  2:121071/151‑1 (MQ=255)                                                                                                                                               gTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATg   >  1:579797/1‑151 (MQ=255)                                                                                                                                                 tGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATgg  >  2:1515317/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGG  >  NC_000913/4191849‑4192142

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