Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,735	0	A	T	40.2%	132.9 / 46.2	82	Y2* (TAT→TAA)	ydaC	DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (5/44);  new base T (33/0);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 5.58e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413586‑1413885                                                                                                                                                      |                                                                                                                                                          aTCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTta                                                                                                                                                           <  1:97985/150‑1 (MQ=255) aTCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTta                                                                                                                                                           <  2:255126/150‑1 (MQ=255) aTCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTta                                                                                                                                                           <  1:1245644/150‑1 (MQ=255)     gTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAt                                                                                                                                                      <  2:787336/151‑1 (MQ=255)          aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTTTATGaa                                                                                                                                                   >  2:150871/1‑149 (MQ=255)              tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGaaaaa                                                                                                                                                >  2:32674/1‑145 (MQ=255)                 gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑c                                                                                                                                            <  1:369107/148‑1 (MQ=255)                 gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCat                                                                                                                                         <  1:809863/151‑1 (MQ=255)                 gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCa                                                                                                                                          <  2:758226/150‑1 (MQ=255)                 gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCa                                                                                                                                          <  2:138789/150‑1 (MQ=255)                    ctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCatt                                                                                                                                        <  1:1207682/149‑1 (MQ=255)                     tctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAAAACCC‑CCcttttt                                                                                                                                     >  1:355622/1‑151 (MQ=255)                      ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAACACCT‑CCaatatt                                                                                                                                     >  1:446491/1‑150 (MQ=255)                      ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCatta                                                                                                                                       <  2:176941/148‑1 (MQ=255)                        ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCattattatt                                                                                                                                  <  2:1494290/151‑1 (MQ=255)                        ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCattatta                                                                                                                                    <  1:833387/149‑1 (MQ=255)                            aCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCa                                                                                                                              <  2:1379112/151‑1 (MQ=255)                             cGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑g                                                                                                                            <  1:13028/151‑1 (MQ=255)                             cGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑g                                                                                                                            <  1:183343/151‑1 (MQ=255)                              gCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATTAAAAACT‑CCCTTAtttttt                                                                                                                                 >  1:114037/1‑146 (MQ=255)                              gCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑g                                                                                                                            <  2:1500505/150‑1 (MQ=255)                               ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAAAACCC‑CCCTTTTTATTTTCA‑Gt                                                                                                                           >  1:1487265/1‑150 (MQ=255)                               ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTg                                                                                                                          <  1:1354585/151‑1 (MQ=255)                                 tGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTAATTTTTTTATTAACACCT‑CCATTAtttttttctctggg                                                                                                                         >  2:1310513/1‑143 (MQ=255)                                  gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATAAAAACAC‑CCATTTTTTTTTTTA‑GTgg                                                                                                                         >  2:1150051/1‑149 (MQ=255)                                  gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGt                                                                                                                        <  2:606942/150‑1 (MQ=255)                                    ttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTACATGAACACCC‑CCCTtttttttttccctgggt                                                                                                                        >  2:421228/1‑140 (MQ=255)                                     tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTc                                                                                                                      <  1:1391631/149‑1 (MQ=255)                                      cttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAACCC‑CCCTTTTTAtttttcatggttcttt                                                                                                                    >  1:1012519/1‑138 (MQ=255)                                        tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGtt                                                                                                                   <  2:518993/149‑1 (MQ=255)                                         ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTa                                                                                                                 <  1:1197189/150‑1 (MQ=255)                                         ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTa                                                                                                                 <  2:347147/150‑1 (MQ=255)                                         ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTa                                                                                                                 <  2:718442/150‑1 (MQ=255)                                            tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTc                                                                                                              <  2:336944/150‑1 (MQ=255)                                             aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAAAACCC‑CCCTTATTTTTTTCC‑GTGGTTCGttttttt                                                                                                              >  1:881415/1‑148 (MQ=255)                                              ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAACACCT‑CCCTTATTTTTTTCA‑GTgggtgttttttttcat                                                                                                           >  2:1304877/1‑137 (MQ=255)                                              ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAACCC‑CCCTTATTTTTTTCA‑CTGGGTTGttttttt                                                                                                              >  1:202304/1‑147 (MQ=255)                                               cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTCAATTTTTACATGAAAACCC‑CCCTTATTATTTTTA‑GTGGGTCGTTTTTTTCAt                                                                                                           >  1:1003099/1‑150 (MQ=255)                                               cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTAAGAAAACCT‑CCATTTTTTTTTTTC‑GTGGTTTGtttttt                                                                                                               >  2:631452/1‑146 (MQ=255)                                                 gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATc                                                                                                          <  2:395428/149‑1 (MQ=255)                                                 gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCtt                                                                                                        <  2:630082/151‑1 (MQ=255)                                                 gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCt                                                                                                         <  1:835737/150‑1 (MQ=255)                                                   gagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATAAAAACCC‑CCCTTTTTTTTTTTA‑TGGGTTCGtttttttttttttt                                                                                                       >  2:752228/1‑142 (MQ=255)                                                   gagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTc                                                                                                      <  2:285535/151‑1 (MQ=255)                                                    agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAACACCT‑CCCTTATTATTTTCA‑GTGGGTCGTTTTTTTCCTCttt                                                                                                       >  1:1358127/1‑149 (MQ=255)                                                        ttttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATTAAAAACT‑CCCTTATTTTTTTTT‑GTGGGTCGtttttttccttttttgtg                                                                                                   >  2:1069864/1‑137 (MQ=255)                                                           tttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTATACATTAACACCC‑CCCTTATTTTTTTCA‑GTGGGTTGTTTTTTTCATCTTTCTAGTGctt                                                                                              >  1:900104/1‑151 (MQ=255)                                                                cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTATAATAAAAAAC‑CCATTATTTTTTTTA‑GGGGTTTGTTTATTTCTTATTTCGAGTGGTtttttt                                                                                         >  2:1228476/1‑151 (MQ=255)                                                                  ttCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGAAAACCC‑CCCTTATTTTTTTTA‑ATGGGTTGTTTTTTTCATCTTTTTTGGGCttttttttt                                                                                       >  1:835958/1‑150 (MQ=255)                                                                         ttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTACATTTTTACAATTATCCCC‑CCACTCTTATTTCCA‑TTGGGTTGTTTTTTTTTTCTTTTGAGCGATTCtttttttctt                                                                                   >  1:1398916/1‑143 (MQ=255)                                                                         ttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTcc                                                                                 <  1:1257174/150‑1 (MQ=255)                                                                           cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTccac                                                                               <  1:771100/150‑1 (MQ=255)                                                                            gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAATAAAACCTCCCATTTTTTTTTCC‑‑GTGGGTCGTTTTTTTTATCTTTCTTGTGGTTTTTTTTTTCTTtcacc                                                                              >  1:1385983/1‑150 (MQ=255)                                                                              tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGAAAACCC‑CCCTTATTTTTTTTA‑GTGGTTCGGTTTTTTCCTCTTTTGAGTGCTTTTTTTTTTCTTccccc                                                                              >  1:467717/1‑148 (MQ=255)                                                                               aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑aa                                                                         <  1:246713/151‑1 (MQ=255)                                                                                 tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAAAACCC‑CCACTATTTTTTTTA‑GGGGTTTTTTTTTTTTATCTTTTGCGAGGttttttttttct                                                                                    >  2:1387701/1‑135 (MQ=255)                                                                                   tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTAAATTAAAACTC‑CCCTTATTTTTTTCATGGGGTTGGTTTTTTTC‑TCTTTTTAGGTCTTTTTTTTTACTTTCACCCT‑aa                                                                         >  2:165669/1‑147 (MQ=255)                                                                                   tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAAGAAAACCC‑CCCTTTTTATTTTCA‑GTGGGTTGTTTTTTTCATCTTTTTAGTGGttttttttttttttccccctaa                                                                          >  2:1339772/1‑133 (MQ=255)                                                                                    aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGt                                                                    <  1:12867/151‑1 (MQ=255)                                                                                      aTTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTc                                                                  <  2:139351/151‑1 (MQ=255)                                                                                        tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATAAAAACCA‑CCCTTTTTTTTTTTA‑GTGGGTCGGTTTTTTTATTTTTTTTGTGGTTTTTTTTTACTTTCCCCCT‑AACCCGt                                                                    >  2:767440/1‑147 (MQ=255)                                                                                          tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAACACCC‑CCCTTATTTTTTTCC‑GTGGTTTGTTTTTTTCCTTTTTTTAGTGCTTTTTTTTTTCTTTCCCCCT‑AACCCGGTCggg                                                               >  1:948901/1‑150 (MQ=255)                                                                                           gAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGa                                                              <  1:153145/150‑1 (MQ=255)                                                                                             atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTc                                                           <  1:229179/151‑1 (MQ=255)                                                                                             atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTc                                                           <  2:435518/151‑1 (MQ=255)                                                                                                     tCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAAGAAAACCC‑CCCTTTTTATTTTCC‑GTGGGTCGTTTTTTTCCTCTTTTGAGTGCTTCTTTTTTTCTTTCCCCCT‑AAACGGTTCGGGGTTTGCATCCGt                                                   >  1:1167818/1‑151 (MQ=255)                                                                                                        aGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAATGAACAACC‑CCATTTTTTTTTTCT‑CTGGGTTGTTTTTTTCCTCTTTTTGGTGGTTCTTTTTTATTTCCCCCCT‑AAACGGTTCTGGGGTCGCCTcc                                                     >  2:719382/1‑146 (MQ=255)                                                                                                          ttAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTAAGAAAAACC‑CCCCTATTTTTTTCA‑GTGGTTTGTTTTTTTTCTCTTTCTTGAGGTTCTTTTTTTCTTTTCCCATAAACCGGTT‑GGGGTTTGCATCCGtt                                                  >  2:1387888/1‑147 (MQ=255)                                                                                                          ttAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAATTTTATTCATGAACACCT‑CCTCTATTATTTTCA‑CTGGTGCGTTTTTTTTCTCGTTTTAGGGTTTCTTTTTTTCTTTCTCCAC‑AACCCGATCGGGGtttttatctctt                                                  >  2:1341273/1‑137 (MQ=255)                                                                                                             gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTc                                           <  1:930450/151‑1 (MQ=255)                                                                                                                 gAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTc                                       <  2:734537/151‑1 (MQ=255)                                                                                                                 gAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTc                                       <  2:1344305/151‑1 (MQ=255)                                                                                                                  aGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCg                                      <  2:204710/151‑1 (MQ=255)                                                                                                                   gTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTCATGAAAAACC‑CCCTTTTTTTTTTCC‑CTGGGTGGTTTTTTCTATCTTTTCAGAGCTTCTTTTTTTTTTTCACCCT‑TACCGGGTTGTGGTTCTCATCCCCTTCACTTct                                          >  2:867337/1‑145 (MQ=255)                                                                                                                            tccttTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATc                              <  1:1494090/149‑1 (MQ=255)                                                                                                                             ccttTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCt                            <  1:423456/150‑1 (MQ=255)                                                                                                                             ccttTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTg                           <  1:124544/151‑1 (MQ=255)                                                                                                                             ccttTTCAATAGTGGCGGTAAATTTTTTAAAGAAAAACT‑CCCTTTTTATTTTCA‑GTGGGTCGGTTTTTTCATCTTTTTCGTGCTTTTTTTTTTCTTCCCCCAT‑AAACCGGTCGGGGGTCGCCACCCCTTCACTTCGCCCGGTAATATcc                             >  1:848169/1‑149 (MQ=255)                                                                                                                                tttCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGaa                         <  1:1211923/150‑1 (MQ=255)                                                                                                                                    actactGGCGGTAATTTTATCCATGAACCCTT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCGCTTACTCCAT‑AACCGGTGCGGGATTCGCATCCTTTCCACTTCGCCCGGTAATATCCTGACATgg                     <  2:607122/145‑1 (MQ=255)                                                                                                                                     aTAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCt                   <  2:36885/151‑1 (MQ=255)                                                                                                                                           gCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAcc             >  2:186508/1‑151 (MQ=255)                                                                                                                                           gCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAcc             >  1:1309141/1‑151 (MQ=255)                                                                                                                                            cGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCg            >  1:264736/1‑151 (MQ=255)                                                                                                                                            cGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTCTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCg            <  1:3125/151‑1 (MQ=255)                                                                                                                                                      ataCATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg  <  1:438746/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                          ATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CCATTATTATTTCCA‑GTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413586‑1413885

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