Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,523,445	T→G	8.2%	D217A (GAC→GCC) ‡	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,523,445	0	T	G	8.2%	192.4 / 10.2	61	D217A (GAC→GCC) ‡	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base G (3/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-02

TCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGC  >  NC_000913/4523307‑4523593                                                                                                                                           |                                                                                                                                                     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                                                                                 >  1:498868/1‑141 (MQ=255)                                                                    tCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCAt                                                                                 <  2:498868/141‑1 (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:308581/1‑149 (MQ=255)                                                                                                       gAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGc                                      >  3:352118/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)                                                                                                                                      gtgcgtCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCAtg       <  4:174342/149‑1 (MQ=255)                                                                                                                                           tCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGTGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGc  <  1:553372/149‑1 (MQ=255)                                                                                                                                           tCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGc  <  1:557681/149‑1 (MQ=255)                                                                                                                                           tCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGc  <  3:266398/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                     TCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGC  >  NC_000913/4523307‑4523593

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: ‑