Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,523,438	G→C	9.8%	R219R (CGC→CGG)	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,523,438	0	G	C	9.8%	202.6 / 12.4	61	R219R (CGC→CGG)	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base C (4/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

CTCCCGCCACACAGGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCA  >  NC_000913/4523291‑4523583                                                                                                                                                    |                                                                                                                                                  ctcCCGCCACACAGGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGgc                                                                                                                                                  <  3:248922/149‑1 (MQ=255)  tcCCGCCACACAGGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGgcg                                                                                                                                                 >  1:332340/1‑149 (MQ=255)   cccGCCACACAGGCTGACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCCCCTGCGATGTCCCTGCAGTGTCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGgcgt                                                                                                                                                <  2:56336/149‑1 (MQ=255)                 tCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAgaga                                                                                                                                  >  2:574351/1‑149 (MQ=255)                     cGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAAt                                                                                                                              >  2:485349/1‑149 (MQ=255)                       gAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTc                                                                                                                            >  4:561528/1‑149 (MQ=255)                               cAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTgg                                                                                                                    <  1:574351/149‑1 (MQ=255)                                aGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGc                                                                                                                   >  1:507763/1‑149 (MQ=255)                                         tgtgGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGTGTGCGTCCTg                                                                                                                                        <  2:558101/119‑1 (MQ=255)                                         tgtgGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGTGTGCGTCCTg                                                                                                                                        >  1:558101/1‑119 (MQ=255)                                         tgtgGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTgc                                                                                                          <  1:485349/149‑1 (MQ=255)                                         tgtgGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCAATGGCGaattatct                                                                                                          <  4:241832/149‑7 (MQ=255)                                               aTTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAAt                                                                                                    <  4:378703/149‑1 (MQ=255)                                                ttGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATg                                                                                                   <  2:475817/149‑1 (MQ=255)                                                 tGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGCACGCTTCGCCAATgg                                                                                                  <  4:25123/149‑1 (MQ=255)                                                  gCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGt                                                                                                 <  2:151046/149‑1 (MQ=255)                                                        cccAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCAc                                                                                           >  4:293450/1‑149 (MQ=255)                                                         ccAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCAcc                                                                                          >  4:443358/1‑149 (MQ=255)                                                         ccAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCAcc                                                                                          >  3:15337/1‑149 (MQ=255)                                                           aaTCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCtt                                                                                        >  1:144276/1‑149 (MQ=255)                                                              ccTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGcc                                                                                     <  1:10333/149‑1 (MQ=255)                                                                ttCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGc                                                                                             >  4:106890/1‑139 (MQ=255)                                                                ttCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGc                                                                                             <  3:106890/139‑1 (MQ=255)                                                                ttCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGt                                                                                   >  3:89019/1‑149 (MQ=255)                                                                     ccaccaCCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtct                                                                              <  4:365515/149‑1 (MQ=255)                                                                      caccacCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCGCCTTGCCGTTGtctt                                                                             <  1:86902/149‑1 (MQ=255)                                                                       accaccTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtcttc                                                                            <  3:537068/149‑1 (MQ=255)                                                                          accTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTcc                                                                         <  1:460031/149‑1 (MQ=255)                                                                           ccTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCa                                                                        <  3:443358/149‑1 (MQ=255)                                                                                    tCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCAt                                                                       <  2:498868/141‑1 (MQ=255)                                                                                    tCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCAt                                                                       >  1:498868/1‑141 (MQ=255)                                                                                     cccGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCgg                                                              <  4:209108/149‑1 (MQ=255)                                                                                      ccGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCggg                                                             >  4:266398/1‑149 (MQ=255)                                                                                       cGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGc                                                            >  1:259380/1‑149 (MQ=255)                                                                                                ccAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTcgcg                                                   >  2:582823/1‑149 (MQ=255)                                                                                                ccAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTcgcg                                                   >  4:245929/1‑149 (MQ=255)                                                                                                            aGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGAcc                                       <  4:205734/149‑1 (MQ=255)                                                                                                            aGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGAcc                                       >  3:174342/1‑149 (MQ=255)                                                                                                                gccgccgGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAAGGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCa                                   <  4:460396/149‑1 (MQ=255)                                                                                                                     cgGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGtt                              <  2:332340/149‑1 (MQ=255)                                                                                                                      ggAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTg                             <  2:263796/149‑1 (MQ=255)                                                                                                                       gAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGc                            >  3:352118/1‑149 (MQ=255)                                                                                                                       gAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGc                            >  3:308581/1‑149 (MQ=255)                                                                                                                          gagGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGAt                         <  1:582823/149‑1 (MQ=255)                                                                                                                                     tttACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTgcg              >  4:137623/1‑149 (MQ=255)                                                                                                                                       tACAGGACGCACGCCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGcgcg            >  4:83067/1‑149 (MQ=255)                                                                                                                                        acaggacGCACGCCGTGCGGGCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCgcgg           >  3:41078/7‑149 (MQ=255)                                                                                                                                         caggacGCACGCCGTGCGGGCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCgcggc          >  2:584095/6‑149 (MQ=255)                                                                                                                                          aGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTcgc                                                    >  3:344881/1‑106 (MQ=255)                                                                                                                                          aGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTcgc                                                    <  4:344881/106‑1 (MQ=255)                                                                                                                                           ggacGCACGCCGTGCGGGCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGAtgtt                      <  4:155184/132‑1 (MQ=255)                                                                                                                                           ggacGCACGCCGTGCGGGCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGAtgtt                      >  3:155184/4‑135 (MQ=255)                                                                                                                                           ggacGCACGCCGTGCGGGCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGa        <  1:584095/146‑1 (MQ=255)                                                                                                                                            cgcGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAg       >  2:474070/3‑149 (MQ=255)                                                                                                                                            cccGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCgg                                            <  2:527528/112‑1 (MQ=255)                                                                                                                                            cccGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCgg                                            >  1:527528/1‑112 (MQ=255)                                                                                                                                             ccGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTcgc                                                    <  3:200249/103‑1 (MQ=255)                                                                                                                                             ccGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTcgc                                                    >  4:200249/1‑103 (MQ=255)                                                                                                                                             ccGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGt      >  4:303013/1‑149 (MQ=255)                                                                                                                                             ccGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGt      <  3:561528/149‑1 (MQ=255)                                                                                                                                              cGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTg     <  4:268114/149‑1 (MQ=255)                                                                                                                                               gCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGc    >  2:188487/1‑149 (MQ=255)                                                                                                                                               gCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGc    >  2:225537/1‑149 (MQ=255)                                                                                                                                                 aCGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCca  >  4:97655/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CTCCCGCCACACAGGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCA  >  NC_000913/4523291‑4523583

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