Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,123	G→T	84.2%	R34S (CGC→AGC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,123	0	G	T	84.2%	227.4 / 41.2	87	R34S (CGC→AGC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base G (9/5);  new base T (44/29);  total (53/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

TCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTC  >  NC_000913/4116979‑4117265                                                                                                                                                 |                                                                                                                                               tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGc                                                                                                                                            <  1:288087/149‑1 (MQ=255) tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGc                                                                                                                                            <  1:364846/149‑1 (MQ=255) tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGc                                                                                                                                            <  1:401129/149‑1 (MQ=255)   gTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGa                                                                                                                                          >  3:189754/1‑149 (MQ=255)        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGc                                                                                                                                     >  1:138183/1‑149 (MQ=255)         ccGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCt                                                                                                                                    >  2:224105/1‑149 (MQ=255)                 tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTaa                                                                                                                                  >  2:642497/1‑143 (MQ=255)                 tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTaa                                                                                                                                  <  1:642497/143‑1 (MQ=255)                 tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat                                                                                                                            >  4:531735/1‑149 (MQ=255)                  gCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATAtt                                                                                                                           <  4:26263/149‑1 (MQ=255)                  gCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATAtt                                                                                                                           >  1:84316/1‑149 (MQ=255)                          aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGt                                                                                                                   >  4:42082/1‑149 (MQ=255)                           tCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTg                                                                                                                  >  3:71306/1‑149 (MQ=255)                             ggAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGAt                                                                                                                >  3:663945/1‑149 (MQ=255)                                  tGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCAtt                                                                                                           <  3:42082/149‑1 (MQ=255)                                   gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTa                                                                                                          <  2:84316/149‑1 (MQ=255)                                    atatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTAc                                                                                                         <  4:477166/149‑1 (MQ=255)                                    atatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTAc                                                                                                         >  4:302225/1‑149 (MQ=255)                                    atatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTAc                                                                                                         <  2:559550/149‑1 (MQ=255)                                      atCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGa                                                                                                                                          <  3:24605/114‑1 (MQ=255)                                      atCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGa                                                                                                                                          >  4:24605/1‑114 (MQ=255)                                       tCGGCTTTCGCCAGCACTTCTACGAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTCTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATCTTGGCATCGTGATTCATTACGAc                                                                                                      >  1:291269/1‑149 (MQ=255)                                         ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcg                                                                                                    >  1:74872/1‑149 (MQ=255)                                            tttCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcgcgc                                                                                                 >  3:588868/1‑149 (MQ=255)                                                        ttCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtgg                                                                                     >  1:402970/1‑149 (MQ=255)                                                        ttCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtgg                                                                                     >  2:494116/1‑149 (MQ=255)                                                         tCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcgcgcga                                                                                               <  4:447113/138‑2 (MQ=255)                                                         tCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCgg                                                                                               >  3:447113/1‑138 (MQ=255)                                                           tACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGc                                                                                  >  4:461044/1‑149 (MQ=255)                                                              cAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCt                                                                               <  3:302225/149‑1 (MQ=255)                                                              cAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCt                                                                               <  2:74872/149‑1 (MQ=255)                                                                 cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGt                                                                            <  3:531735/149‑1 (MQ=255)                                                                      aGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgc                                                                       <  4:92150/149‑1 (MQ=255)                                                                         tttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCa                                                                    >  3:631541/1‑149 (MQ=255)                                                                         tttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCa                                                                    >  3:658425/1‑149 (MQ=255)                                                                                ggCCCAGATTTCCCTTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGttat                                                             <  4:588868/149‑1 (MQ=255)                                                                                ggCCCAGATTTCCATTGGGTCGTGTTCTACACAGCCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGGCACGCTAATGAGATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGatat                                                             <  2:291269/149‑1 (MQ=255)                                                                                   ccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGGTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtt                                                          >  2:591874/1‑149 (MQ=255)                                                                                   ccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACTATATAtt                                                          >  3:594581/1‑149 (MQ=255)                                                                                   ccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtt                                                          >  3:143819/1‑149 (MQ=255)                                                                                   ccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtt                                                          >  3:46513/1‑149 (MQ=255)                                                                                                   tCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGtt                                          >  2:365189/1‑149 (MQ=255)                                                                                                   tCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGtt                                          >  3:179716/1‑149 (MQ=255)                                                                                                      tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTaa                                       <  1:224105/149‑1 (MQ=255)                                                                                                        ttCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAAtt                                     >  1:363073/1‑149 (MQ=255)                                                                                                        ttCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTTGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAAtt                                     >  4:537615/1‑149 (MQ=255)                                                                                                           tACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTc                                  <  4:179716/149‑1 (MQ=255)                                                                                                           tACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTc                                  <  2:282499/149‑1 (MQ=255)                                                                                                            aCCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                                 >  4:665028/1‑149 (MQ=255)                                                                                                            aCCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                                 >  4:385008/1‑149 (MQ=255)                                                                                                             cccAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTccc                                <  2:402970/149‑1 (MQ=255)                                                                                                             cccAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTccc                                >  3:154057/1‑149 (MQ=255)                                                                                                              ccAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCg                               >  3:564878/1‑149 (MQ=255)                                                                                                              ccAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCg                               >  2:599565/1‑149 (MQ=255)                                                                                                                aagttGGTGTTAGCTAGCTTGGCTCCAATTCGCTTTCCGACACGCTAATGATATTGGCATCGTTATCCATTACGACCGCGCGGGAGCTGGTGGTGGCCTGGTCGAGAGCAACGATATAGTTTTTTTCAGTCATGTTTAATGGTCTCGTa                             >  2:106720/5‑149 (MQ=255)                                                                                                                aaCCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTa                             <  3:385008/149‑1 (MQ=255)                                                                                                                aaCCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTa                             >  3:683249/1‑149 (MQ=255)                                                                                                                aaCCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTa                             <  2:561017/149‑1 (MQ=255)                                                                                                                  ccTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGt                           <  4:189754/149‑1 (MQ=255)                                                                                                                   cTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                                 >  2:22601/1‑142 (MQ=255)                                                                                                                   cTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                                 <  1:22601/142‑1 (MQ=255)                                                                                                                   cTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                                 >  1:693481/1‑142 (MQ=255)                                                                                                                   cTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                                 <  2:479458/142‑1 (MQ=255)                                                                                                                   cTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                                 >  1:479458/1‑142 (MQ=255)                                                                                                                   cTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                                 <  2:693481/142‑1 (MQ=255)                                                                                                                    tGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCa                         >  1:74371/1‑149 (MQ=255)                                                                                                                     ggTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCat                        <  2:138183/149‑1 (MQ=255)                                                                                                                       ttttGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCatat                      >  4:353657/1‑149 (MQ=255)                                                                                                                       ggttGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCatat                      <  1:264386/147‑1 (MQ=255)                                                                                                                           gggTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACa                  >  1:339722/1‑149 (MQ=255)                                                                                                                             gTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGc                <  1:591874/149‑1 (MQ=255)                                                                                                                             gTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGc                >  1:52217/1‑149 (MQ=255)                                                                                                                             gTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGc                <  4:143819/149‑1 (MQ=255)                                                                                                                                 aTTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAg            >  1:643509/1‑149 (MQ=255)                                                                                                                                 aTTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAg            >  3:234415/1‑149 (MQ=255)                                                                                                                                  tttGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGc           >  1:623477/1‑149 (MQ=255)                                                                                                                                  tttGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGc           >  1:104145/1‑149 (MQ=255)                                                                                                                                     gCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttt        <  1:365189/149‑1 (MQ=255)                                                                                                                                     gCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttt        <  4:71306/149‑1 (MQ=255)                                                                                                                                       tCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt      >  4:37398/1‑149 (MQ=255)                                                                                                                                       tCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt      >  1:62563/1‑149 (MQ=255)                                                                                                                                       tCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt      <  3:37398/149‑1 (MQ=255)                                                                                                                                       tCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt      >  4:666197/1‑149 (MQ=255)                                                                                                                                       tCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt      >  1:41962/1‑149 (MQ=255)                                                                                                                                       tCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt      >  2:333931/1‑149 (MQ=255)                                                                                                                                         aaaTTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGt    >  3:203683/1‑149 (MQ=255)                                                                                                                                           aTTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTc  >  3:111628/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                               TCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTC  >  NC_000913/4116979‑4117265

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