Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,735	0	A	T	45.7%	71.5 / 53.3	70	Y2* (TAT→TAA)	ydaC	DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (8/30);  new base T (32/0);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.34e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTG  >  NC_000913/1413594‑1413881                                                                                                                                              |                                                                                                                                                       tACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑Gaa                                                                                                                                               <  1:267232/150‑1 (MQ=255)      tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCAT‑GAACACCTc                                                                                                                                         >  1:516154/1‑151 (MQ=255)         gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAT‑TAAAACcccc                                                                                                                                        >  2:1077331/1‑149 (MQ=255)         gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTc                                                                                                                                         <  2:1164301/148‑1 (MQ=255)             tctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑TAACACCCCCattttt                                                                                                                                  >  1:405420/1‑151 (MQ=255)               tGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑TAACACCCCCctttttat                                                                                                                                >  1:249708/1‑151 (MQ=255)                ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCattattatt                                                                                                                               <  1:852165/151‑1 (MQ=255)                ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCattatt                                                                                                                                  <  1:881301/148‑1 (MQ=255)                ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTACAT‑GAACAACCCCctttttatt                                                                                                                               >  1:385584/1‑151 (MQ=255)                   aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTATATTTTTACAT‑TAAAACCACCCTTCTTAtttt                                                                                                                             >  1:796696/1‑149 (MQ=255)                     cGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTATACAT‑TAACAACTCCCTTAtttttt                                                                                                                              >  1:944177/1‑147 (MQ=255)                       ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑GAACACCCCCATTAtttttt                                                                                                                              >  2:345838/1‑145 (MQ=255)                       ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑g                                                                                                                       <  1:837978/151‑1 (MQ=255)                          gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGtt                                                                                                                    <  2:53807/151‑1 (MQ=255)                          gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTATAT‑TGACAGATCCTTTATTAATTCCaattggt                                                                                                                     >  1:1095735/1‑145 (MQ=255)                            ttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑TAACACCCCCCTTATTATTTTCAGT‑gggt                                                                                                                    >  1:339558/1‑147 (MQ=255)                               ttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAT‑TAAAACCCCCCTTTTTTTTTTCCGT‑GGGTTGttt                                                                                                               >  1:535212/1‑151 (MQ=255)                                tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAT‑TAACACCCCCCTTTTTTTTTTCAGT‑GGGTTGtttt                                                                                                              >  2:348376/1‑150 (MQ=255)                                tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGttt                                                                                                               <  2:1071333/150‑1 (MQ=255)                                tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGtt                                                                                                                <  1:90260/149‑1 (MQ=255)                                 ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTa                                                                                                              <  1:313730/150‑1 (MQ=255)                                  tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTATAT‑GAAAACCCCCCTTTTTATTTTCAAT‑gggttttttttt                                                                                                            >  2:972385/1‑141 (MQ=255)                                    tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTc                                                                                                           <  1:131321/150‑1 (MQ=255)                                       cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑GAAAACCTCCCTTATTTTTTTCCGT‑GGGTCGtttttttc                                                                                                          >  1:1131518/1‑146 (MQ=255)                                          agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAT‑TAAAACCCCCCTTTTTTTTTTCCGG‑GGTTTGTTTTTTTCATCttt                                                                                                    >  2:980864/1‑151 (MQ=255)                                                ttttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTGATTTTTTTAAT‑GACAACCCCCATTTTTTTTTTCAGT‑GGTTCGTTTTTTTCATtttttccggt                                                                                              >  2:845532/1‑145 (MQ=255)                                                 tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAT‑GAACACCCCCCTTTTTATTTTCAGT‑GGTTCGTTTTTTTCATCTTTTGAggg                                                                                              >  1:184164/1‑148 (MQ=255)                                                 tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTg                                                                                              <  1:1002020/150‑1 (MQ=255)                                                      cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAAT‑GAACACCCCCCTTATTTTTTTCCGT‑GGGTTGTTTTTTTCATCTTTTGAGTGcttttt                                                                                        >  1:605743/1‑151 (MQ=255)                                                       gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCttt                                                                                       <  1:878361/151‑1 (MQ=255)                                                         aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt                                                                                     <  2:401192/151‑1 (MQ=255)                                                           tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc                                                                                    <  1:22316/150‑1 (MQ=255)                                                            cTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAAATTTTTATAT‑GAAAACCACCCCTTTTTTTTTTAGT‑GGGTTGTTTTTTTCATCCTTTTTGTGCTTCttttt                                                                                     >  2:681175/1‑148 (MQ=255)                                                              gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAT‑GAAAACCTCCCTTATTATTTTCCGT‑GGGTCGTTTTTTTCATCCTTTCAGTGCTTCTTTTTTAC‑tt                                                                               >  2:101892/1‑151 (MQ=255)                                                               aTTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTc                                                                              <  2:472947/151‑1 (MQ=255)                                                                    gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTAT‑GAACACCCCCCTTTTTTTTTTCAGT‑GGTTCGTTTTTTTCATTTTTTGCGTGGttttttttttctt                                                                                >  2:582674/1‑140 (MQ=255)                                                                    gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGCAATTTTTAATATGGAAAACCCCCATTTTTTTTTTCACT‑‑GATGGCTTTTTTTATCTTTTTGGTGCTTCttttttgctt                                                                                >  2:874002/1‑140 (MQ=255)                                                                     ttAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑GAACACCCCCCTTTTTTTTTTCAGT‑GGGTTGTTTTTTTCCTTTTTTTAGTGCTTTTTTTTTTC‑TTtcac                                                                           >  2:501967/1‑148 (MQ=255)                                                                      tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTAT‑GAACACCCCCCTTTTTTTTTTCCTT‑GGTTTGTTTTTTTTATCTTTTTGGTGCTTCttttttattttcccc                                                                           >  2:877747/1‑138 (MQ=255)                                                                       aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCAT‑TAAAACCCCCCTTTTTTTTTTCACT‑GGTTTGTTTTTTTTATCTTTTTAGGGCTTTTTTTTTTC‑TTccacc                                                                          >  2:731446/1‑147 (MQ=255)                                                                        gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAc                                                                     <  2:393429/151‑1 (MQ=255)                                                                        gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAAATTTTTATAA‑TAAAACCTCCCCTTTTTTTTTTCGT‑GGTGTGTTTTTTTCATCTTTTTGGGGGTTCtttttttttt                                                                                >  2:695697/1‑136 (MQ=255)                                                                          ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAca                                                                    <  2:993109/150‑2 (MQ=255)                                                                            aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAGTTTTTTAAAT‑GAACAACCCCCTTTTTATTTTCAGG‑GGGTTGTTTTTTTCATCTTTTTGGGGGTTTTTTTTTAC‑TTTCACCATAAC‑Cgg                                                                 >  1:823974/1‑150 (MQ=255)                                                                            aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑Cgg                                                                 <  2:504079/150‑1 (MQ=255)                                                                             gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGtt                                                               <  1:233030/151‑1 (MQ=255)                                                                                tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAT‑GAACACCCCCCTTTTTTTTTTCAGT‑GGGTTGTTTTTTTCATTTTTTTAGTGCTTTTTTTTTTCTTTCCCCCTCAAC‑CGGtt                                                               >  2:978553/1‑149 (MQ=255)                                                                                tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTAGTTAT‑TAACCCCTCCCTTATTTTTTTCCGT‑GGTTTGTTTATTTCATTTTTCTTTGGCTTCTTTTTTAA‑ATTCACCAAAAA‑Cgg                                                                 >  1:799874/1‑146 (MQ=255)                                                                                 ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCgg                                                            <  2:1028285/150‑1 (MQ=255)                                                                                    aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAA‑GAAAACCTCCCTTTTTATTTTCATT‑GGTTTGTTTTTTTCATCTTTTTTGTGCTTTTTTTTTTC‑TTTCACCCTAAA‑CCGTTTGGGGtt                                                        >  2:97730/1‑151 (MQ=255)                                                                                       atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑TAAAACCCCCCTTTTTTTTTTCAGT‑GGGTTGTTTTTTCCCTCTTTTTGGGGGTTCTTTTTTTC‑TTCCCCCCTAAC‑CCGGTCGGGGttt                                                       >  1:111114/1‑148 (MQ=255)                                                                                        tCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑TAACACCCCCCCTTTTTTTTTCAGTGGGTTGGTTTTTTCC‑TTTTTCCGGGGGTTCTTTTTTAC‑TTTCCCCCTAACACGGTTggggtt                                                         >  2:1158488/1‑145 (MQ=255)                                                                                            gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATcc                                                 <  1:1109264/150‑1 (MQ=255)                                                                                              ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑GAACACCTCCATTTTTTTTTTCAAA‑GGTTCGTTTATTTCATCTTTCGAGTGGTTCTTTTTTAC‑TTTCCCCATAAA‑CCGGTCGGGGTTCGCATCCGtt                                              >  1:506618/1‑151 (MQ=255)                                                                                              ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑GAACACCCCCCTTATTTTTTTCCGT‑GGGTCGGTTTTTTCCTCTTTTTAGTGCTTCTTTTTTTC‑TTCCCCCACAAA‑CGGGTCGGGGTTTTCATCCGtt                                              >  2:340272/1‑151 (MQ=255)                                                                                                aGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTAATATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTa                                             <  1:752812/150‑1 (MQ=255)                                                                                                     gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTc                                       <  1:567304/151‑1 (MQ=255)                                                                                                             gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGt                                <  2:176822/150‑1 (MQ=255)                                                                                                               cACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAat                             <  2:484248/151‑1 (MQ=255)                                                                                                                 ccttccttTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAatat                           <  1:795270/151‑1 (MQ=255)                                                                                                                   ttccttTTCAATAGTGGCGGTAATTTTTTAAAA‑GAAAACCTCCCTTTTTATTTTCCGG‑GGGTTGTTTATTTCATCTTTTGTGTGGTTTTTTTTTAC‑TTCCCCCATAAA‑CGGGTTGGGGGTCGCCTCCGTTACACTTCCCTCGGGaa                              >  1:286315/1‑146 (MQ=255)                                                                                                                    tccttTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATcc                         <  2:54771/150‑1 (MQ=255)                                                                                                                     ccttTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCGC‑TTCCGCCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTg                       <  1:822161/151‑1 (MQ=255)                                                                                                                          tCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAAt                   <  1:798149/150‑1 (MQ=255)                                                                                                                            aaTAGTGGCGGTAATTTTTTAAAT‑GAACACCCCCATTATTATTTTCAGT‑GGGTTGTTTTTTTCATCTTTCGAGTGCTTCTTTTTTTC‑TTTCACCATAAC‑CCGGTTGGGGTTTTCCTCCCTTTCACTTCGCTCCGTGATATTCTGaaa                    >  2:441941/1‑147 (MQ=255)                                                                                                                                     ggTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTCTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt       >  1:805115/1‑151 (MQ=255)                                                                                                                                      gTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt      >  2:243409/1‑151 (MQ=255)                                                                                                                                        aaTTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTAc    >  1:246020/1‑151 (MQ=255)                                                                                                                                         aTTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACt   <  2:545084/151‑1 (MQ=255)                                                                                                                                           tttATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTg  >  1:100097/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                       TACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAAC‑CGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTG  >  NC_000913/1413594‑1413881

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