Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	689,514	0	A	C	39.0%	102.7 / 11.0	59	G456G (GGT→GGG)	lnt	apolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (32/4);  new base C (0/23);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.22e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689364‑689661                                                                                                                                                         |                                                                                                                                                    aCTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCa                                                                                                                                                     <  2:910220/151‑1 (MQ=255)   tcatcaACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGGGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCAcc                                                                                                                                                   >  2:179707/1‑151 (MQ=255)                 gcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTC‑CGGGGGTCGGCGTCA‑CGTTAGAGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATc                                                                                                                                      >  1:94977/1‑150 (MQ=255)                      aCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCAccgcc                                                                                                                                <  1:129887/151‑1 (MQ=255)                             aaTGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGc                                                                                                                         >  1:19524/1‑151 (MQ=255)                               tGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTC‑CGGTGGTCGGCGGCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCg                                                                                                                       >  2:961441/1‑151 (MQ=255)                                 ccGTCAGCCCCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTC‑CGGTGGTGGGCGTCA‑CGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCgtt                                                                                                                     <  2:17521/151‑1 (MQ=255)                                        cgcccACAGCGGCCAGTTGCCGTTACGTGCGTATGGGGTGAGCC‑GGGTGGCCGGCGTCA‑CGTTAGTGGTTAACCCCTCGCGGGTGACCGGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGttgttg                                                                                                                 <  1:228590/146‑1 (MQ=255)                                          ccaccggCGGCCGTTGGCCGGTACGTGCTGATGGTGTGAGTC‑GGGTGGTCGGCTTCC‑CGTTAGTGTTAACCACCCCGCGGGTGACCTGCGGAATCATCGCCTGCATCCCCCCTTGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGt                                                                                                               <  1:298518/142‑1 (MQ=255)                                            cacaGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtgctgc                                                                                                          >  2:116395/1‑151 (MQ=255)                                             acaGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTgcg                                                                                                         >  2:575763/1‑151 (MQ=255)                                              caGCGGCCAGTTGCCGGTACGTGCGTATGGTGTAGGTCTCGGCGG‑CCGCCCCA‑CCATTGAGACTAACACCTCGCGGGTGAACTGCGCGATCACCGCCTGACTCTCACCCTGCGCGCCAATCACCCCCGTACCGCCGTCGTCGGtgctcc                                                                                                          >  1:13763/1‑147 (MQ=255)                                                gCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCa                                                                                                       <  2:986366/150‑1 (MQ=255)                                                  ggCCAGTTGCCGGTACGTGCGTATGGTGTAGGTCTCGGTGGTCGCCGT‑A‑CGTTTAAAGAAAACAAATCGCGGGGGAAATGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCCACGCCGTAATGCCGTTGTTGGTGCTGCGCa                                                                                                       >  1:530254/1‑148 (MQ=255)                                                    ccAGTTGCCGGTACGGGCGTATGGTGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGt                                                                                                  <  2:672388/151‑1 (MQ=255)                                                          gCCGGTACGTGCGTATGGTGTGAGTC‑CGGTGGTCGGCGTCC‑CGTTAGTGGTTACCACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgc                                                                                              <  2:599907/149‑1 (MQ=255)                                                           ccGGTACGTGCGTATGGTGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgc                                                                                              <  2:656500/148‑1 (MQ=255)                                                           ccGGTACGGGCGTATGGTGTGAGTC‑GGGTGGTCGGCGTCA‑CGTTAGTGTTTACCCCCTCGCGGGTGACCTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgc                                                                                            <  2:241940/150‑1 (MQ=255)                                                            ggtgaCGTGCGTATGGTGTGTGCC‑CGGTGGTGGCCGAAA‑CGTTATGGGTAAACACCGCGTGGGTGAACGGCGGCATCACCCCCTGAATCTCCCCCGGGGGTCAAATCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgcg                                                                                           <  1:642141/146‑1 (MQ=255)                                                            cGGTACGTGCGTATGGTGTGAATC‑CGGTGGGCGGCGTCA‑CGTTAATGGTTAACACCTCGCTGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTGATGCCGTTGTTGGTGCTGCGCAACAgcg                                                                                                 >  1:63993/1‑142 (MQ=255)                                                               tACGTGCGTATGGTGTGGGTC‑TGGGGGGCGGCGGCA‑CGATAGAGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGcc                                                                                         >  1:754899/1‑149 (MQ=255)                                                                aCGTGCGTATGGTGTGAGTC‑CGGGGGGCGGCGTCA‑CCTTAATGGATAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGc                                                                                      >  1:398107/1‑151 (MQ=255)                                                                  gTGCGTATGGTGTAAGTC‑TGGTGGGCCGCCGCATCATTA‑TGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCCGCTc                                                                                    >  1:706361/1‑151 (MQ=255)                                                                    gCGTATGGTGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTcc                                                                                   <  2:980064/150‑1 (MQ=255)                                                                      gtttGGTGTGAGTC‑GGGTGGTGGGCGTCC‑CGTTAGTGTTTAACCCCTCGCGGGTGAACTGCGGGATCTCCGCCTAAACCTCCCCCGGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAg                                                                                       <  2:1035277/141‑1 (MQ=255)                                                                          ggtttgatgc‑cGGGGGTCGGGTGCA‑CGTTAGGGGTTACCCCCTCGCGGGTGAACTGGGGGATCATCGCCTGAATCTCACCCTGCGGCCCAATCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgca                                                                             <  1:812830/141‑1 (MQ=255)                                                                          ggTGTGAGTC‑GGGTGGTCGGGGTCA‑CGTTAGTGGTTAACCCCCCGGGGGTAAACTGCGGGTTCACCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAgc                                                                                <  2:1118133/147‑1 (MQ=255)                                                                          ggTGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcac                                                                            <  2:1054153/151‑1 (MQ=255)                                                                             gtgAGTC‑CGGTGGTGGGCGTCA‑CGTTAGTGGTTAACCCCTCGCGGGTGAACTGGGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcacgc                                                                          <  2:671623/150‑1 (MQ=255)                                                                                    cgggtggtCGGCGTCA‑CTTTAGTGGTTAACCCCCCGCGGGGAAACGGCGGGATCACCCCCTAAATCCCACCCTGGGGGCCAACCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcg                                                                     <  1:720077/147‑1 (MQ=255)                                                                                       gtggtCGGCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGcc                                                                 >  1:991910/1‑150 (MQ=255)                                                                                        tggtCGGCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCAt                                                               >  2:702508/1‑151 (MQ=255)                                                                                              gCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGa                                                          >  2:1114095/1‑150 (MQ=255)                                                                                               cgccA‑CGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCACCGCCTGAATCTCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAg                                                        <  1:476196/148‑1 (MQ=255)                                                                                               cGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGa                                                          >  1:594278/1‑149 (MQ=255)                                                                                               cGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGa                                                          >  2:207817/1‑149 (MQ=255)                                                                                                   a‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGtt                                                    >  2:1043663/1‑151 (MQ=255)                                                                                                      gTTAGTGGTTAACCCCTCGGGGGTGAACTGCGGGTTCACCGCCTGAATCTCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGc                                                  <  2:1046498/151‑1 (MQ=255)                                                                                                       ttAGTGGTTAACCCCTCGCGGGGAAACTGCGGGATCATCGCCGGAATCTCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGtt                                                    <  2:312080/148‑1 (MQ=255)                                                                                                       ttAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGcc                                                 >  2:612501/1‑151 (MQ=255)                                                                                                       ttAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGcc                                                 >  2:7185/1‑151 (MQ=255)                                                                                                          gTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATg                                              >  2:796699/1‑151 (MQ=255)                                                                                                           tGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATgg                                             >  1:970350/1‑151 (MQ=255)                                                                                                            ggTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATgg                                             <  2:633607/150‑1 (MQ=255)                                                                                                                  ccccTCGCGGGTGAACTGCGGGATCATCGCCTGAATCCCCCCCTGCGGGCCAATCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATa                                      <  1:373330/149‑1 (MQ=255)                                                                                                                   acCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAAt                                       >  1:787458/1‑149 (MQ=255)                                                                                                                          gggTAACCTGCGGAATCACCGCCGTAATCTCCCCCGGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAttt                                 <  1:702178/148‑1 (MQ=255)                                                                                                                             tGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCaa                            >  2:93977/1‑150 (MQ=255)                                                                                                                                 cTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCAc                       >  1:258841/1‑151 (MQ=255)                                                                                                                                   gCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGc                     >  2:255231/1‑151 (MQ=255)                                                                                                                                   gCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGc                     >  2:823910/1‑151 (MQ=255)                                                                                                                                   gCGGGATCACCGCCGGAATCCCCCCCTGGGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAGCAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCa                        <  1:739813/148‑1 (MQ=255)                                                                                                                                     gggATCACGGCCTGATTCTCCCCCAGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCAt                   <  1:1046351/151‑1 (MQ=255)                                                                                                                                      ggaccgccGCCTGAATCTCCCCCTGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATc                  <  1:532785/144‑1 (MQ=255)                                                                                                                                      ggATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATc                  >  1:604298/1‑151 (MQ=255)                                                                                                                                            tcGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTg              <  1:701228/149‑1 (MQ=255)                                                                                                                                              gCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             >  2:432206/1‑148 (MQ=255)                                                                                                                                               ccTGAATCTAACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             >  1:743146/1‑147 (MQ=255)                                                                                                                                                    aTCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTc     >  2:959321/1‑150 (MQ=255)                                                                                                                                                        cACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTcagc  >  2:907256/1‑149 (MQ=255)                                                                                                                                                         |                                                                                                                                                    ACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTC‑CGGTGGTCGGCGTCA‑CGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689364‑689661

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