Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,735	0	A	T	50.6%	57.5 / 133.5	79	Y2* (TAT→TAA)	ydaC	DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (4/35);  new base T (40/0);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 2.53e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-11

CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAT  >  NC_000913/1413588‑1413883                                                                                                                                                    |                                                                                                                                                         cAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTtata                                                                                                                                                        <  1:615306/150‑1 (MQ=255) cAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTtata                                                                                                                                                        <  1:427861/150‑1 (MQ=255)   gTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAt                                                                                                                                                     >  1:60527/1‑151 (MQ=255)        aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAac‑a                                                                                                                                              <  2:635611/151‑1 (MQ=255)          tttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑GGac‑a                                                                                                                                              >  1:276996/1‑149 (MQ=255)            tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTGAATTTTTTCAT‑TAAC‑A‑CCTc                                                                                                                                         >  1:29765/1‑151 (MQ=255)               gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTc                                                                                                                                         <  1:295004/148‑1 (MQ=255)               gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCat                                                                                                                                      <  1:65491/151‑1 (MQ=255)                  ctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCatt                                                                                                                                     <  1:551369/149‑1 (MQ=255)                   tctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAT‑GAAA‑A‑CCTCCcttatt                                                                                                                                  >  2:725578/1‑151 (MQ=255)                    ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCattat                                                                                                                                   <  1:260642/149‑1 (MQ=255)                    ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCattat                                                                                                                                   <  2:46254/149‑1 (MQ=255)                      ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAA‑TTAA‑A‑CCTCCcttattttt                                                                                                                               >  1:199726/1‑151 (MQ=255)                        aaaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTc                                                                                                                             <  2:816622/151‑1 (MQ=255)                          aCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCa                                                                                                                           <  1:876453/151‑1 (MQ=255)                             ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑g                                                                                                                       <  2:1048254/151‑1 (MQ=255)                              cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGGATTTTTTTAAT‑AAAA‑A‑CCCCCCTCATTTTTTc                                                                                                                             >  2:347842/1‑145 (MQ=255)                                gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGtt                                                                                                                    <  1:1029056/151‑1 (MQ=255)                                   tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCg                                                                                                                  <  2:591880/150‑1 (MQ=255)                                   tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTATAT‑GAAA‑A‑CCTCCCTTATTATTTTCCGT‑GGGTCGt                                                                                                                 >  2:956006/1‑151 (MQ=255)                                    cttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAATTTTATTAAA‑AAAC‑A‑ACACCAATATAATTTTCatggggtgggt                                                                                                                 >  1:998964/1‑141 (MQ=255)                                     ttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAT‑AAAA‑A‑CCCCCCTttttttttttaggtggtcggt                                                                                                                 >  1:325123/1‑137 (MQ=255)                                      tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑TAAA‑A‑CCTCCACTATTTTTTTCTCT‑GGTTCGgtt                                                                                                               >  2:202252/1‑150 (MQ=255)                                      tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGtt                                                                                                                <  1:1116126/149‑1 (MQ=255)                                         ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGCAAATTTTTAAAT‑CAAA‑A‑CCTTCCTTATTTTTTcctctgggccgcttttt                                                                                                             >  1:781059/1‑135 (MQ=255)                                          tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTc                                                                                                           <  1:393431/150‑1 (MQ=255)                                            ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑TAAC‑A‑ACTCCCTTATTTTTTTCAAT‑GGtttttttttt                                                                                                            >  1:842456/1‑139 (MQ=255)                                            ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑GAAC‑A‑CCCCCCTTATTTTTTTCAGT‑GGTTCGtttttt                                                                                                            >  2:220988/1‑147 (MQ=255)                                               gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATc                                                                                                       <  2:620403/149‑1 (MQ=255)                                               gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCtt                                                                                                     <  2:609107/151‑1 (MQ=255)                                                  agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAA‑GAAA‑ACCCCCCATTTTTTTTTCCA‑T‑GGGTCGtttttttccttttt                                                                                                    >  1:1055153/1‑141 (MQ=255)                                                  agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTTAT‑TAAA‑A‑CCCCCCTTTTTTTTTTCCGT‑GGTTCGTTTTTTTCATCttt                                                                                                    >  1:675861/1‑149 (MQ=255)                                                   gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑GAAA‑A‑CCTCCATTATTATTTTCCGT‑GGTTTCTTTTTTTCATCttt                                                                                                    >  1:1088546/1‑148 (MQ=255)                                                     aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAA‑TAAC‑A‑CCCTCCTTTTTTTTTTCAGT‑GGTTCGTTTTTTTTATCttttgtg                                                                                                >  1:197991/1‑146 (MQ=255)                                                       tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAT‑GAAA‑A‑CCCCCCTTTTTTTTTTCAAT‑GGGTTGttttttttctctttttagg                                                                                               >  1:457605/1‑136 (MQ=255)                                                         tttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTAATTTTTAAATCGAAA‑A‑CTTCCATT‑TTTTTTCTTCT‑GGTTTGttttttttattttttca                                                                                                 >  1:39184/1‑134 (MQ=255)                                                             gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAA‑GTAC‑A‑ACCCCCATATTTTTTTTACT‑GGGTTGTTTTTTTCATCttttttggggtttttt                                                                                       >  1:553988/1‑138 (MQ=255)                                                             gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCttt                                                                                       <  2:10718/151‑1 (MQ=255)                                                              cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAT‑TGAC‑A‑CCCCCCTTATTTTTTTTAAT‑GGGTCGTTTTTTCCATCTTTTTGGTGCTTCtttt                                                                                      >  1:610317/1‑151 (MQ=255)                                                              cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCtttt                                                                                      <  2:184695/151‑1 (MQ=255)                                                               aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTTCAT‑GTAC‑A‑CCTCCCCTTATATTTTCATT‑GGTGTGTTTTTTCCCCCTTTTCAGGGCTTCttttt                                                                                     >  2:940102/1‑151 (MQ=255)                                                                ttCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCAT‑TAAA‑A‑CCCCCCTTTTTATTTTCAGT‑GGTTTGTTTTTTTCATCTTTTTTGTGCTTCtttttt                                                                                    >  1:66392/1‑150 (MQ=255)                                                                 tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCAT‑GATT‑A‑ACCCCCTTTTTATTTTCTGG‑GGTTTGCTTTTTTCATTTTTTGTGTTCTTCtttttt                                                                                    >  1:639640/1‑149 (MQ=255)                                                                   tGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑TAAA‑A‑CCTCCCTTATTTTTTTCAGT‑GGGTCGTTTTTTTCCTCTTTTTAGTGCtttttttt                                                                                     >  2:307597/1‑147 (MQ=255)                                                                   tGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAAA‑TAAC‑A‑CCCCCCTTTTTTTTTTCCGG‑GGGTCGTTTTTTTCCTCTTTCTAGGGCttttttttttct                                                                                 >  2:340241/1‑147 (MQ=255)                                                                    gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAT‑GAAA‑A‑CCCCCATTATTATTTTCAGT‑GGGTCGTTTTTTTCATCTTTTTTGTGCTTTTTTTTTACtt                                                                                >  2:817561/1‑151 (MQ=255)                                                                    gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑TAAA‑A‑CCCCCCTTATTTTTTTCAGT‑GGGTTGTTTTTTTCATTTTTTTCGTGCTTCtttttt                                                                                    >  1:1102542/1‑146 (MQ=255)                                                                    gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACtt                                                                                <  2:277043/151‑1 (MQ=255)                                                                        tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTcca                                                                             <  1:938732/150‑1 (MQ=255)                                                                         cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑GAACAA‑CCCCCTTTTTTTTTTCCAG‑‑GGTTCTTTTTTTTCCTCTTTTTTGTGCttttttttttctttccc                                                                            >  1:231961/1‑141 (MQ=255)                                                                         cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTccac                                                                            <  1:288873/150‑1 (MQ=255)                                                                            tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAt                                                                         <  1:480617/150‑1 (MQ=255)                                                                             aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATa                                                                        <  2:897948/150‑1 (MQ=255)                                                                                ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAT‑GAAC‑A‑CCCCCCTTATTATTTTTCGT‑GGGTGGTTTTTTTTCTCTTTCTAGTGCTTCTTTTTTCTTTCCCCCATaa                                                                       >  1:770561/1‑148 (MQ=255)                                                                                 tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑GAAA‑A‑CCTCCCTTTTTATTTCCCCA‑GGGTCGTTTATTCCTTATTTCGTGGGGTTTTTTTTTTCTtcctccccaaacgg                                                                   >  2:998002/1‑143 (MQ=255)                                                                                    aTTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTc                                                                <  2:27216/151‑1 (MQ=255)                                                                                     ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑GAAC‑A‑ACCCCCTTTTTATTTTCCAT‑GGGTCGTTTTTTTCATCTTTCGAGTGGTTCTTTTTTACTTTCACCCTaaac                                                                     >  2:1030422/1‑143 (MQ=255)                                                                                      tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑TAAC‑A‑CCTCCCTTATTTTTTTCCGT‑GGTTCGTTTTTTCCATCTTTTTAGTGCTTTTTTTTTACTTTCACCATAAACCGtt                                                                 >  2:206754/1‑148 (MQ=255)                                                                                         gAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGCAATTTTTTTCAT‑GAAA‑A‑CCTCCCTTTTTTTTTTCCTT‑GGGTTGGTTTTTTCATCTTTCCAGTGGTTTTTTTTTACTTCCTCCCTAACCTGTTCggggt                                                           >  1:753098/1‑149 (MQ=255)                                                                                                  gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTGCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATcc                                                   <  1:908140/150‑1 (MQ=255)                                                                                                    ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑GAAC‑A‑ACCCCCTTATTATTTTCAGT‑GGTTCGGTTTTTTCCTCTTTCTAGTGCTTCTTTTTTACTTCCACCCTAAACGGTTCGGGATTCGCATcc                                                   >  1:412842/1‑148 (MQ=255)                                                                                                    ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGtt                                                <  1:713303/151‑1 (MQ=255)                                                                                                    ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGt                                                 <  2:514303/150‑1 (MQ=255)                                                                                                      aGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAA‑GAAC‑A‑CCCCCCTTATTATTTTCAGT‑GGTTTGTTTTTTTCATCTTTTGAGTGCTTCTTTTTCACTTTCACCATAAACGGTTTGGGATTTGGATCCCtt                                                >  1:84570/1‑149 (MQ=255)                                                                                                        ttAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAAT‑GAAA‑A‑CCTCCCTTATTTTTTTCAGT‑GGGTCGTTTTTTTCATCTTTTTAGTGCTTTTTTTTTTCTTCCCCCATAACCGGTTCGGGGTTCGCATCCCtt                                                >  2:743510/1‑147 (MQ=255)                                                                                                         tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTacac                                            <  1:32882/150‑1 (MQ=255)                                                                                                           gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                         <  2:496210/151‑1 (MQ=255)                                                                                                           gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                         <  2:652392/151‑1 (MQ=255)                                                                                                            tagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                         <  1:1493/150‑1 (MQ=255)                                                                                                              ggAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAT‑GAAC‑A‑CCCCCCTTTTTATTTTCAGT‑GGGTCGTTTTTTTCATCTTTTTAGTGCTTTTTTTTTTCTTTCCCCATAAACGGGTTGGGGTTTGCCTCCCTTTCACtttc                                        >  2:235438/1‑147 (MQ=255)                                                                                                                    ccACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑GAAC‑A‑CCTCCCTTTTTATTTTCCGT‑GGTTCGTTTTTTTCATCTTTCTAGTGCTTTTTTTTTTCTTTCCCCCTAAACGGTTCTGGGTTTGCATCCGTTTCACtttcctc                                     >  1:551416/1‑141 (MQ=255)                                                                                                                        cttccttTTCAATAGTGGCGGTAATTTTTTACAT‑AAAA‑A‑CCTCCCTTATTTTTTTCCGT‑GGGTTGTTTTTTTCATCTTTTTGGTGCTTTTTTTTTTCTTTCACCCTAAACGGTTCGGGGTTTGCCTCCCTTtcaccttcctccggta                                >  2:1089109/1‑135 (MQ=255)                                                                                                                        cttccttTTCAATAGTGGCGGTAATTTTTTAAAA‑GAAA‑A‑CCTCCCTTATTTTTTTCCGTGGGTTGGTTTTTTCC‑TCTTTTTGGTGCTTTTTTTTTACTTTCACCCTAACCCGTTCGGGGTTTGCCTCCCTTTCACTTTCCTCCGTtata                              >  1:896449/1‑149 (MQ=255)                                                                                                                         ttccttTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCCTTATTTTTTCCAGT‑GGGTTGTTTTTTTCATCTTTTGAGTGCTTCTTTTTTACTTTCACCATAAACCGGTCGGGGTTCGCATCCGCTACACTTTGCTCggga                                 >  2:1035572/1‑143 (MQ=255)                                                                                                                           ccttTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCGGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTg                         <  1:183924/151‑1 (MQ=255)                                                                                                                               ttCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGaaa                      <  2:488244/150‑1 (MQ=255)                                                                                                                                   aTAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCt                 <  2:66397/151‑1 (MQ=255)                                                                                                                                         gCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAcc           >  1:364685/1‑151 (MQ=255)                                                                                                                                            gTAATTTTATACAT‑GAGC‑G‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCGCCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCGCGGTAATATCCTGAAATGGCTTTCTCCGtt        <  1:656483/151‑1 (MQ=255)                                                                                                                                              aaTTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTa       <  2:505300/150‑1 (MQ=255)                                                                                                                                                  ttATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt  >  1:993231/1‑151 (MQ=255)                                                                                                                                                    ataCAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt  <  2:502306/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                         CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAT  >  NC_000913/1413588‑1413883

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