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

CATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689367‑689661                                                                                                                                                        |                                                                                                                                                    catcaACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑TCGGGGGCGGC‑GGCATCATAGAGGATAACACCTCGCGGGGGAAC‑TGCGGGATCATCGCCTGAATCTCAcc                                                                                                                                                   >  2:1377566/1‑150 (MQ=255)      acacAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCGCCCTGCgg                                                                                                                                             <  1:1302565/151‑1 (MQ=255)      acacAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTAGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCgg                                                                                                                                             <  2:1213123/151‑1 (MQ=255)         cagcagCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGcc                                                                                                                                          <  2:1308571/151‑1 (MQ=255)           gcagcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCa                                                                                                                                         <  2:953873/150‑1 (MQ=255)            cagcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCaa                                                                                                                                        <  1:673327/150‑1 (MQ=255)               caAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACCCCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAAGCTCCCCCTGCGGGCCAATCAc                                                                                                                                    <  2:935197/151‑1 (MQ=255)               caAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCAc                                                                                                                                    <  1:1522559/151‑1 (MQ=255)                 aaaCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCAccg                                                                                                                                  <  1:949183/151‑1 (MQ=255)                 aaaCCAAACAATGCCGGCATCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCGCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCAccg                                                                                                                                  <  2:1120289/151‑1 (MQ=255)                       aaCAAGCCCGTCAGCACCCACAGCGGCCAGTTGCCGGTGCGTGCGTATGTTGT‑GAGCC‑GGGTGGTCGGC‑GCCACGTTAGTGGTTACCCCCCCGCGGGTGAAC‑GGCGGGATCATCGCTTGAACCTCCCCCTGCGGGCCAATCACCGCCGTaa                                                                                                                            <  2:202865/151‑1 (MQ=255)                            tGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCg                                                                                                                       >  1:1073909/1‑151 (MQ=255)                              ccGTCAGCACCCACAGGGGCTGGTTGCGGGTAGGTGCGTATGGTGG‑GGGTC‑CGGTGGTCGGC‑GCCACGTTAGTGGTTACCCCCTCGCGGGTGAAC‑TGCGGAATCACCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCgtt                                                                                                                     <  2:300006/151‑1 (MQ=255)                                  cAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑GGGTGGTCGGC‑GTCACGTTAGTGGTTAACCCCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGttgttg                                                                                                                 <  1:448460/151‑1 (MQ=255)                                  cAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑AAATC‑CGGGGGTCGGC‑GTCATGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCgttgtt                                                                                                                  >  2:1127032/1‑150 (MQ=255)                                   aGCCCCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACCCCTCGCGGGTGAAC‑TGGGGGATCATCGCCTGAATCTCCCCCTGCGGGCAAATCACCGCCGTAATGCCgttgtt                                                                                                                  <  1:1277354/149‑1 (MQ=255)                                    gcccccACGGCGGCCGTTTGCCGTTCCGTGCGTATGTGGG‑GGTTC‑CGGTGGTCGGC‑TCCACGTTAGTGGTTACCCCCTCGCGGGTAAAC‑TGCGGATTCATCGCCTGAACCCCCCCCTCCGGCCCAATCACCGCCGTAATGCCGTTGTTgg                                                                                                                <  2:1635773/147‑1 (MQ=255)                                     cccccGCAGCGGCCGGTGGCGGTTGCGTGGGTTTGGGGT‑AATTC‑CGGTGGTCGGG‑TTCCTGTTAGTGGTTACCCCCTCGCGGGTGACC‑TGGGGAATCATCTCGTGAACCTCCCCCTGCGGCCCATTCACCGCCGTAATGCCGTTGTTGGtg                                                                                                              <  2:1131755/149‑1 (MQ=255)                                     cACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑AAGTC‑CCGGGGGCGGC‑GGTACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGGTGTTGGt                                                                                                               >  2:29763/1‑150 (MQ=255)                                        ccACAGCGGCCAGTTGCGGGTAGGTGCGTAGGGTGT‑GATGC‑CGGTGGTCGGT‑GTCACGTTAGTGTTTACCCCCGCGCGGGTGAAC‑GGCGGGATCCTCGCCTAAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCgttgtt                                                                                                                  <  2:261432/144‑1 (MQ=255)                                          acaGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGGCGGC‑GGCACGTTAGTGGTTAACACATCGCGGGGGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTCGTTGGggctg                                                                                                           >  2:1291788/1‑149 (MQ=255)                                          acaGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑AGGTCTCGGGGGT‑GGC‑CGCATGATTGAAAATAACAACTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTgcg                                                                                                         >  2:549774/1‑151 (MQ=255)                                           caGCGGCCAGTTGCCGGTACGTGCGTATGGTGTAGATCT‑CGGTGG‑CGGCGGTATCATTA‑AAGATAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATATCATCCTGCGGGACAATCACCGCCGTAATGCAGTTGTTGGAGATgcgc                                                                                                        >  2:1403171/1‑151 (MQ=255)                                             gCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCa                                                                                                       >  1:339404/1‑150 (MQ=255)                                                gCCAGTTGCCGGTGCGTGCTTATGGTGG‑GGGTC‑CGGTGGTCGGC‑GCCACGTAAGTGGTTAACACCACGGGGGTGACC‑GGGGGGATCTTCGCTTGAACCTCACCCTCCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCa                                                                                                       <  2:626774/147‑1 (MQ=255)                                                  cAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTg                                                                                                 >  1:143158/1‑151 (MQ=255)                                                   aGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTgg                                                                                                >  1:1027768/1‑151 (MQ=255)                                                      tGCCGGTACGTGCGTATGGTGTAGATCT‑CGGTGGTCGCC‑GTATCATTA‑AAAAAAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCATGCGGGCCAATCACCGCCGAAATGCCGTTGTTGGTGCTGCGCCACAGt                                                                                                  >  2:1551920/1‑146 (MQ=255)                                                      tGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcg                                                                                             >  1:1043740/1‑151 (MQ=255)                                                        ccGGTCCGTGCTTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGGGGTTAACCCCGCGCGGGGGAAC‑TGGGGTATCACCGCCTGAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcg                                                                                             <  1:1252302/149‑1 (MQ=255)                                                        ccGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgcg                                                                                           >  1:197135/1‑151 (MQ=255)                                                               gTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTc                                                                                    >  1:397041/1‑151 (MQ=255)                                                                tGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTcc                                                                                   >  1:1474746/1‑151 (MQ=255)                                                                 gCGTATGGTGT‑GAGCC‑CGGTGGTCGGC‑GCCACGTTAGTGGTAAACGCCTCGCGGGTGAAC‑GGCGGGATCACCGCCGGAAGCCCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTc                                                                                    <  1:1030819/149‑1 (MQ=255)                                                                   ggaTGGTGT‑GACTC‑GGGTGTGGGGC‑GCCCCGTTAGTGTTAAACCCCCCGCGGGTGACC‑GGCGGGATCCTCGCCGTATTCCCCCCCGGCGGGCCAACCGCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCa                                                                                  <  2:1130109/147‑1 (MQ=255)                                                                     aTGGTGT‑GAGCC‑GGGTGGTCGGC‑GTCACGTTAGTGGTTAACCCCTCGCGGGTAAAC‑TGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAgcg                                                                               <  1:1168775/150‑1 (MQ=255)                                                                       ggTGT‑GAGTC‑GGGTGGTCGGC‑GTCACGTTAGTGTTTAACCCCTCGCGGGTGAAC‑TGGGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgca                                                                             <  1:765791/150‑1 (MQ=255)                                                                         tgt‑gAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGTGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcacgc                                                                          >  2:627218/1‑151 (MQ=255)                                                                           t‑gAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCAt                                                                        >  1:1244294/1‑151 (MQ=255)                                                                              aGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTc                                                                      >  2:1280818/1‑151 (MQ=255)                                                                                gtcggttggtggGC‑GCCATGTGAGGGGTAAACGCCTCGGGGGTGAAC‑TGCGGGTGCCTCGCCTTAATCTCCCCTGGCGGGCGACACTCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTc                                                                      <  2:737535/144‑1 (MQ=255)                                                                                        gtCGGC‑GTCGCTTTAGTGTTTAACACCTCGCGGGGGAAC‑GGCGGGTTCACCGCCGGAATCTCACCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCt                                                             <  2:1310570/151‑1 (MQ=255)                                                                                        gtCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCt                                                             >  1:363625/1‑151 (MQ=255)                                                                                                 cACGTTAGTGGTTAACCCCTCGCGGGTGACC‑TGCGGGATCATCGCCTGAATCCCCCCCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAg                                                        <  2:165770/148‑1 (MQ=255)                                                                                                 cACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtgt                                                     >  2:1147436/1‑151 (MQ=255)                                                                                                  aCGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGtt                                                    >  2:1603800/1‑151 (MQ=255)                                                                                                     ttAGTGGTTAACACCTAGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGcc                                                 >  2:1190280/1‑151 (MQ=255)                                                                                                         tGGTTAACACCTGGCGGGTGACC‑TGCGGGACCATCGCCTAAATCCCACCCTGCGGCCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATg                                              <  1:1532401/150‑1 (MQ=255)                                                                                                            ttagctccgcgcgGGTGAAC‑GGCGGCATCACCTCGGGAATCCCCCCCGGCGGCCCAACCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGaa                                                         <  2:1159674/127‑1 (MQ=255)                                                                                                             taccacCTCGCGGGTGAAC‑GGGGGGATCCTCGCCTGAATCTCGCCCTGCGGGACAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGAcc                                          <  1:1588822/147‑1 (MQ=255)                                                                                                              aaCCCCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGAcc                                          <  1:1490851/149‑1 (MQ=255)                                                                                                              aaCACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCa                                         >  2:34548/1‑150 (MQ=255)                                                                                                                ccccgcgcgGGGGAAC‑TGCGGGATCATCGCCTGAACCCCCCCCTGCGGGCCAACCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGAc                                           <  2:1385727/141‑1 (MQ=255)                                                                                                                ccccTCGGGGGTGAAC‑TGCGGGATCATCGCCGGAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAAt                                       <  1:1441264/148‑1 (MQ=255)                                                                                                                  ccTCGCGGTTGAAC‑TGGGGAATCATCCCCGGAATCTCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATa                                      <  1:341609/149‑1 (MQ=255)                                                                                                                  ccTCGCGGGTAAAC‑TGCGGGATCATCGCCTGAATCTCACCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAg                                     <  1:1514819/150‑1 (MQ=255)                                                                                                                    tCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAtt                                  >  1:1632379/1‑151 (MQ=255)                                                                                                                      gcgGGTGAAC‑TGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAt                                   <  1:702000/148‑1 (MQ=255)                                                                                                                      gcgGGTGAAC‑GGCGGTATCCTCGCGTGAATCCCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAtt                                  <  2:1626138/149‑1 (MQ=255)                                                                                                                      cggggTGGACTTGCGGGATCACCGCCTGAGTCCCCCCTGGCGGCCCAATCACCGCCGTAATGCCGTGGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCaa                                        <  2:1224981/142‑1 (MQ=255)                                                                                                                        ggctgcaC‑TGCGGTTTCATCGCCTGACTCTCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATGTa                                <  2:523642/143‑1 (MQ=255)                                                                                                                          gTGAAC‑TGCGGGATCATCGCCTGATTCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCaa                            <  2:1026934/151‑1 (MQ=255)                                                                                                                             aaC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAAcc                         >  2:1188398/1‑151 (MQ=255)                                                                                                                             aaC‑TGCGGAATCATCCCCTGAAGCTCCCCCTGCGGGCCAACCGCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCaa                            <  2:1073050/148‑1 (MQ=255)                                                                                                                                    gggATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCAt                   >  1:92979/1‑151 (MQ=255)                                                                                                                                    gggATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCAt                   >  1:681624/1‑151 (MQ=255)                                                                                                                                          atcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             >  2:616326/1‑151 (MQ=255)                                                                                                                                            cGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             >  1:1612296/1‑149 (MQ=255)                                                                                                                                              ccggAATCTCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGAt         <  1:778574/148‑1 (MQ=255)                                                                                                                                               ctaaaTCCCCCCTTGAGGGCCAATCTCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             <  1:799577/143‑1 (MQ=255)                                                                                                                                                       cACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTcagc  >  2:157082/1‑149 (MQ=255)                                                                                                                                                        aCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTcagc  >  1:283972/1‑148 (MQ=255)                                                                                                                                                        |                                                                                                                                                    CATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACGTTAGTGGTTAACACCTCGCGGGTGAAC‑TGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689367‑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: ‑