Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	689,514	0	A	C	44.9%	112.5 / 18.2	70	G456G (GGT→GGG)	lnt	apolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (31/7);  new base C (0/31);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 3.25e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689371‑689661                                                                                                                                                |                                                                                                                                                    aaCACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGc                                                                                                                                               <  3:1276438/150‑1 (MQ=255) aaCACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCg                                                                                                                                              >  3:914831/1‑151 (MQ=255)  acacAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCgg                                                                                                                                             >  3:770274/1‑151 (MQ=255)    acaGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGc                                                                                                                                           <  1:153088/151‑1 (MQ=255)      agcagcAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCa                                                                                                                                         <  4:545559/151‑1 (MQ=255)        cagcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAAt                                                                                                                                       >  3:853429/1‑151 (MQ=255)          gcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCa                                                                                                                                     <  3:1015573/151‑1 (MQ=255)            aaaaCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCAcc                                                                                                                                   <  3:89735/151‑1 (MQ=255)                ccAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCAccgccg                                                                                                                               <  1:1008929/151‑1 (MQ=255)                    aCAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAAt                                                                                                                           >  2:175215/1‑151 (MQ=255)                      aaTGCCGTCAGCCCCCACGGCGGCCAGTGGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGGCGGCGTCACTTTAGTGGTTACCCCCTCGCGGGTGACCTGGGGGATCACCGCCTGAATCCCCCCCTGGGGGCAAATCACCGCCGTAATGc                                                                                                                         <  1:914708/151‑1 (MQ=255)                          ccGTCAGCCCCCACAGGGGCCAGTTGCGGGTCCGTGGGTATGGTGTGAGTCCGGTGGTGGGCGTCACTTTAGTGTTAAACCCCCCGCGGGTAAACTGGGGGATCATCGCCTGAATCTCCCCCGGCGGGCCAATCACCGCCGTAATGCCg                                                                                                                       <  1:1175522/149‑1 (MQ=255)                                gCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTgg                                                                                                                <  2:460255/150‑1 (MQ=255)                                 cccccACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtg                                                                                                              <  2:328854/149‑1 (MQ=255)                                 cACCCACAGCGGCCAGTTGCCGGTACGTGCTTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCGCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCAAATCACCGCCGTAATGCCGttgttg                                                                                                                 <  3:171531/148‑1 (MQ=255)                                   cccACAGCGCCCAGTGGCCGGTACGTGCTTGTGGGGTGAGTCCGGTGGTCGGCGTCACGTTAGTGTTTAACCCCCCGCGGGTAAACGGCGGGATCATCGCCTGAATCCCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGt                                                                                                               <  1:1297868/148‑1 (MQ=255)                                    ccACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGGGGGCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtg                                                                                                              >  1:420706/1‑148 (MQ=255)                                         gCGGCCAGTGGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGCCACGTTAGTGGTAAACCCCTCGCGGGTGAACTGCGGTATCACCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTgcg                                                                                                         <  3:151783/148‑1 (MQ=255)                                             ccAGTTGCCGGTACGTGCGTGTGGTGTGAGTCCGGTGGTGGGGATCAGGTGATTGTTTAACCCCTCGCGGGTGAACGGCGGAATCAGCCCCGAATTCTCCCCCGGGGGTCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACa                                                                                                    <  3:232359/149‑1 (MQ=255)                                             ccAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTAGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTGATGCCGTTGTTGGTGATACGCAACAGt                                                                                                  >  1:202594/1‑151 (MQ=255)                                                  tGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcg                                                                                             >  1:290995/1‑151 (MQ=255)                                                   gCCGGTGCGTGCGTAGGGTGTGAGTCCGTTGGTCGGGGTCACGTTAGTGGTTAACCCCGCGCGGGTAACCTGGGGGATAATCGCCTGAATCTCCCCCTGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTgg                                                                                                <  4:276557/147‑1 (MQ=255)                                                      ggTACGTGCTGATGGTGTGAGCCGGGTGGGCGGCGTCACGTTAGTGGTTAACCCCTCGGGGGTGAACTGCGGGATCATCGCCTGAATCTACCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgcgc                                                                                          <  2:191080/150‑1 (MQ=255)                                                           gTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGGGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGc                                                                                      <  2:372674/149‑1 (MQ=255)                                                           gTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTc                                                                                    >  4:286726/1‑151 (MQ=255)                                                             gCGTATGGTGTGAGTCCGGTGTTGGGCGTCACGTTGGTGGTTAACCCCTCGCGGGTGAACTGGGGGTTCATCGCTTAATTTTCTCCCGGCGGCCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTc                                                                                    <  2:184687/149‑1 (MQ=255)                                                              cGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGGGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTcc                                                                                   <  2:416990/149‑1 (MQ=255)                                                                  tGGTGTGAGCCCGGTGGGCGGCACCATGTTAGTGGTAAACCCCTCGGGGGAAAACTTCGGGACCATCGCCTGAGCCGCCCCCTGCGGGCCAGCCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCt                                                                                     <  4:287559/143‑1 (MQ=255)                                                                    gtgtGAGTCGGGTGGTCGGGGCCACTTTAGTGTTTAACCCCTCGCGGGGAAGCTGGGGGATCAGCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAgcgc                                                                              <  3:110091/148‑1 (MQ=255)                                                                      gtgAGTCCGGTGGTCGGCGTCACGTTAGTGGTAAACCCCTCGCGGGTAAACTGCGGGATCTTCGCCTGAATCTCCCCCTGCGGCCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcacg                                                                           <  4:402591/149‑1 (MQ=255)                                                                       tgAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcacgca                                                                         >  3:163803/1‑150 (MQ=255)                                                                        gAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCAtt                                                                       >  2:403547/1‑151 (MQ=255)                                                                            ccGGTGGTCGGCGTCACGTTGGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTAAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcgcg                                                                   <  3:1046234/151‑1 (MQ=255)                                                                                tggtCGGCCTCACTTTAGTGTTTAACCCCTCGGGGGTGACCTGCGGGATCATCGCCTGATTCTCCCCCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGcc                                                                 <  3:104908/149‑1 (MQ=255)                                                                                  gtCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCt                                                             >  3:362656/1‑151 (MQ=255)                                                                                     ggCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTgg                                                           >  3:97189/1‑150 (MQ=255)                                                                                     ggCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTgg                                                           >  3:506719/1‑150 (MQ=255)                                                                                     ggCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTgg                                                           >  1:1201518/1‑150 (MQ=255)                                                                                     ggCGTCACGTTAGGGTTTAACCCCTCGCGGGGGAACTGGGGGATCCCCGCCTGAATCCCCCCCCGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTgg                                                           <  2:240296/150‑1 (MQ=255)                                                                                        gTCACGTTAGTGTTTACCACCTCGGGGGTGAACTGGGGGATCATCGCCTAAATCTCCCCTTGCGGGCCAATCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgt                                                       <  2:146802/151‑1 (MQ=255)                                                                                         tCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATc                                                              >  4:9070/1‑143 (MQ=255)                                                                                          cagcttggTGGTTAACACCTGGGGGGTGACCTGCGGGATATACGCATAAACCTCCCCCTGAGGGCAAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGaa                                                         <  3:453948/140‑1 (MQ=255)                                                                                          cACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATc                                                              <  2:9070/142‑1 (MQ=255)                                                                                           aCGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGtt                                                    >  3:1378244/1‑151 (MQ=255)                                                                                           aCGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGtt                                                    >  2:87534/1‑151 (MQ=255)                                                                                              ttAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGcc                                                 >  1:803302/1‑151 (MQ=255)                                                                                                 gTGGTTAACCCCTCGGGGGTGAACTGCGGGATCACCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCAt                                               <  2:33506/150‑1 (MQ=255)                                                                                                     ttAACCCCTGGCGGGGGACCGGGGGGATCATCGCGTGAATCCCCCCCGGGGGGCCACTCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGa                                            <  3:686532/149‑1 (MQ=255)                                                                                                     ttAACCCCTCGCGGGTGAACGGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGAc                                           <  3:649708/150‑1 (MQ=255)                                                                                                     taaaCCCCCCGCGGGTAAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGa                                            <  3:158796/147‑1 (MQ=255)                                                                                                             tCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGTGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAtt                                  >  1:799235/1‑151 (MQ=255)                                                                                                               gggggTGACCTGGGGGATCATCGCCTGAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAtt                                  <  1:279876/147‑1 (MQ=255)                                                                                                               gggggTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAttt                                 <  3:827527/148‑1 (MQ=255)                                                                                                               gcgGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAttt                                 >  4:227804/1‑150 (MQ=255)                                                                                                                cgGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTAc                               >  1:975466/1‑151 (MQ=255)                                                                                                                  ggTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTAcc                              >  3:108775/1‑150 (MQ=255)                                                                                                                    tGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCaa                            >  1:624595/1‑150 (MQ=255)                                                                                                                     gAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAAc                          >  1:1391650/1‑151 (MQ=255)                                                                                                                      aaCTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAAcc                         <  3:941189/151‑1 (MQ=255)                                                                                                                      aaCTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGACGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAAc                          >  2:389959/1‑150 (MQ=255)                                                                                                                         tGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACg                      >  4:107265/1‑151 (MQ=255)                                                                                                                         tGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACg                      >  1:440154/1‑151 (MQ=255)                                                                                                                          gggggATCATCGCCTGAATCTCCCCCGGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACg                      <  2:296773/148‑1 (MQ=255)                                                                                                                            ggtatcatcCCCTGAATCCCCCCCTGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCAt                   <  3:104614/148‑1 (MQ=255)                                                                                                                            gggATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCAt                   >  3:935086/1‑151 (MQ=255)                                                                                                                                tcatcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGtt               >  1:200657/1‑151 (MQ=255)                                                                                                                                 caccGCCTAAACCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGtt               <  3:88944/147‑1 (MQ=255)                                                                                                                                       cgaaaTCCCCCCCTGCGGACCAATCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATa        <  1:1378096/148‑1 (MQ=255)                                                                                                                                           aTCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGCGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTc     >  1:434251/1‑150 (MQ=255)                                                                                                                                              tcACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTcagc  >  3:127437/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                    AACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689371‑689661

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