Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,751	0	A	C	40.5%	105.1 / 29.8	75	N59K (AAT→AAG)	ralR	Rac prophage; restriction alleviation protein
Reads supporting (aligned to +/- strand):  ref base A (9/35);  new base C (30/0);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 5.87e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATG  >  NC_000913/1413602‑1413899                                                                                                                                                         |                                                                                                                                                      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:652392/151‑1 (MQ=255)                                                                                             gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                                         <  2:496210/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)                                                                                                                                        aCAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGAc              <  1:955884/151‑1 (MQ=255)                                                                                                                                         cAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACa             <  1:571736/151‑1 (MQ=255)                                                                                                                                             gAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAAcc          <  1:435824/151‑1 (MQ=255)                                                                                                                                                     cTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACaaa     >  2:936973/1‑150 (MQ=255)                                                                                                                                                       ccATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAAt   <  2:198022/150‑1 (MQ=255)                                                                                                                                                       ccATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATg  <  1:486496/151‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                      GCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAAC‑A‑CCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATG  >  NC_000913/1413602‑1413899

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