Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,523,446	C→G	7.4%	D217H (GAC→CAC) ‡	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,523,446	0	C	G	7.4%	216.4 / 6.8	65	D217H (GAC→CAC) ‡	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand):  ref base C (29/31);  new base G (3/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

TCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCC  >  NC_000913/4523307‑4523594                                                                                                                                            |                                                                                                                                                     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                                                                                                     >  3:15337/1‑149 (MQ=255)                                         ccAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCAcc                                                                                                     >  4:443358/1‑149 (MQ=255)                                           aaTCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCtt                                                                                                   >  1:144276/1‑149 (MQ=255)                                              ccTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGcc                                                                                                <  1:10333/149‑1 (MQ=255)                                                ttCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGc                                                                                                        <  3:106890/139‑1 (MQ=255)                                                ttCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGc                                                                                                        >  4:106890/1‑139 (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                                                  >  3:174342/1‑149 (MQ=255)                                                                                            aGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGAcc                                                  <  4:205734/149‑1 (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                                 >  3:155184/4‑135 (MQ=255)                                                                                                                           ggacGCACGCCGTGCGGGCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGAtgtt                                 <  4:155184/132‑1 (MQ=255)                                                                                                                           ggacGCACGCCGTGCGGGCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGa                   <  1:584095/146‑1 (MQ=255)                                                                                                                            cgcGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAg                  >  2:474070/3‑149 (MQ=255)                                                                                                                            cccGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCgg                                                       >  1:527528/1‑112 (MQ=255)                                                                                                                            cccGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCgg                                                       <  2:527528/112‑1 (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:225537/1‑149 (MQ=255)                                                                                                                               gCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGc               >  2:188487/1‑149 (MQ=255)                                                                                                                                 aCGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCca             >  4:97655/1‑149 (MQ=255)                                                                                                                                      gtgcgtCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCAtg        <  4:174342/149‑1 (MQ=255)                                                                                                                                           tCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGTGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGc   <  1:553372/149‑1 (MQ=255)                                                                                                                                           tCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGc   <  3:266398/149‑1 (MQ=255)                                                                                                                                           tCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGc   <  1:557681/149‑1 (MQ=255)                                                                                                                                            ccTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGcc  <  3:303013/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                     TCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCC  >  NC_000913/4523307‑4523594

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