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

TCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689369‑689661                                                                                                                                                    |                                                                                                                                                    tcaACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTg                                                                                                                                                >  1:25368/1‑151 (MQ=255)  caACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGc                                                                                                                                               <  3:1201054/151‑1 (MQ=255)       cagcagcagAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGcc                                                                                                                                          <  2:760042/151‑1 (MQ=255)        agcagcAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGcc                                                                                                                                          >  2:875491/1‑150 (MQ=255)         gcagcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCaa                                                                                                                                        <  3:198642/151‑1 (MQ=255)          cagcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAAt                                                                                                                                       >  2:514604/1‑151 (MQ=255)           agcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATc                                                                                                                                      <  2:677286/151‑1 (MQ=255)           agcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATc                                                                                                                                      <  2:357176/151‑1 (MQ=255)              aaaaCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCAcc                                                                                                                                   <  2:594902/151‑1 (MQ=255)              aaaaCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCAcc                                                                                                                                   <  2:1348168/151‑1 (MQ=255)                 aCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCAccgcc                                                                                                                                >  1:418560/1‑151 (MQ=255)                 aCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCAccgcc                                                                                                                                <  1:61561/151‑1 (MQ=255)                     aaCAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTa                                                                                                                             >  4:24557/1‑150 (MQ=255)                         aTGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGcc                                                                                                                        >  3:143164/1‑151 (MQ=255)                             cGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCgttg                                                                                                                    >  3:969655/1‑151 (MQ=255)                              gTCAGCCCCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGCCGGGTGGTGGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACGGCGGGATCATCGCCTGAA‑TCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCgttg                                                                                                                    <  4:304287/150‑1 (MQ=255)                                   cACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTACCCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtg                                                                                                              <  2:765260/151‑1 (MQ=255)                                          aGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTgcgc                                                                                                        >  1:407592/1‑150 (MQ=255)                                             ggcatgTTGCCGGTCGGTGCGTGT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTGGTGGTTAACCCCTCGGGGGTGAATTGGGGGTTCATCGCCTGAA‑TCCCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACa                                                                                                    <  4:25368/146‑1 (MQ=255)                                               ccAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGt                                                                                                  >  3:792080/1‑151 (MQ=255)                                                cAGTTGCCGGTACGTGCGTAT‑GGTGGGATTCGGGTGGTCGGCGCCACTTTAGTGGTTACCCCCCCGGGGGTGAACGGGGGGATCATCGCTTGAA‑TCCCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGt                                                                                                  <  3:542359/150‑1 (MQ=255)                                                   ttGCGGTTACGTGCGTAT‑GGTGTGAGTCGGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCACCGCCTGAA‑TCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTgg                                                                                                <  2:281588/149‑1 (MQ=255)                                                     gCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgc                                                                                            >  2:974133/1‑151 (MQ=255)                                                      ccGGTACGTTGGTAG‑GGTGTGGGTCGGGTGGGCGGCGTCACTTTAGTGTTTAACCCCCCGCGGGTGAACGGGGGGTTCACCGCCTGAA‑TCCCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgc                                                                                              <  3:541425/148‑1 (MQ=255)                                                       cggttggtgtgtATGGGTGGGAGCCCGGTGGTCGGCGTCGCGTAAATGGTTAACCCCTGGGGGGTAAATGGGGGTACCTGCGCTGTAT‑TCTCCCCCGGGGGCCGAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgc                                                                                              <  3:791029/142‑1 (MQ=255)                                                        ggTGCGTGGGTTG‑GGTGTGAGTCGGGTGGTGGGGGTCACGTTAGTGGTTAACCCCTGGCGGGTGAACTGGGGGATCACCGCCTGAA‑CCCCCCCCGGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgcgc                                                                                          <  3:1261204/150‑1 (MQ=255)                                                             gTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCCCTTTGGTGGTTACCCCCTCGCGGGTGAACGGCGGAATCACCGCCTGAA‑TCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTc                                                                                    <  2:765006/151‑1 (MQ=255)                                                                 gTAT‑GGTGTGGGTCGGGTGGCCGGCGTCCCGTTGGTGTTTACCCCCTCGGGGTGGAACTGCGGGATCCCCCCCTAAA‑TCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAgc                                                                                <  3:1067021/151‑1 (MQ=255)                                                                      ggtgtggtgccGTTGGTGGGCGCCCCGTTGGTGTTTCACACGCTGCGGGAGACCTTCGGGATCATCTCCCGAA‑TCTACCTCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAGCAGTGGGCGCGCCAGCTCCAGcgca                                                                             <  2:718961/141‑1 (MQ=255)                                                                         gggggTCCGGTGGTCGGGGTCACGTAATTGTTAAACCCCTCGGGGTTGACCTGGGGAATCATCGCCTAAA‑TCTCCCCCTGCGGGCCAAGCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgca                                                                             <  3:1035868/143‑1 (MQ=255)                                                                           gAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCCCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCAtt                                                                       <  2:131871/151‑1 (MQ=255)                                                                             gTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTc                                                                      <  4:322923/150‑1 (MQ=255)                                                                              tCCGGTGGTCGGCGTCCCGTTAGTGGTTAACCCCTCGGGGGTGAACGGCGGGATCATCGCCTGAA‑TCCCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcg                                                                     <  2:1189136/150‑1 (MQ=255)                                                                                 ggtggtCGGCGTCACTTTGGTGGTTAACCCCCGGCGGGTGAACTGCGGGATCTCCCCCTGAT‑TCCCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcg                                                                     <  3:4758/147‑1 (MQ=255)                                                                                  gtggtCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGcc                                                                 >  3:452508/1‑150 (MQ=255)                                                                                    ggtCGGCGTCACGTTAGTGGTTACCCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCAt                                                               <  4:156356/150‑1 (MQ=255)                                                                                     gtCGGCTTCACGTTAGTGTTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAC‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCt                                                             <  2:454533/151‑1 (MQ=255)                                                                                        ggCGTCTCGTTAGTGGTTAACCCCTCGCGGGTGACCTGCGGGATCATCGCCTGAA‑TCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTg                                                            <  2:1108027/149‑1 (MQ=255)                                                                                        ggCGTCCCGTTAGTGGTTACCACCCCGGGGGTGAACTGCGGGATCACCGCCTGAA‑TCTCCCCCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTg                                                            <  4:243682/149‑1 (MQ=255)                                                                                        ggCGCCACTTTAGGGGTTAACCCCCCGCGGGTGAGCTGGGGGATCATCGCCTGAA‑CCCCCCCCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGa                                                          <  3:761559/151‑1 (MQ=255)                                                                                            tCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtg                                                      >  3:187945/1‑151 (MQ=255)                                                                                                 tgaGTGTTTCACCCCCCGCGGGTGAACTGCGGGATCATAGCCTAAATTCTC‑CCCGGCGGCCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGaa                                                         <  3:813994/141‑1 (MQ=255)                                                                                                    gtgtTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCCCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATg                                              <  4:452039/151‑1 (MQ=255)                                                                                                      ggTTACCCCCTCGCGGGTGACCTGGGGGATCATCGCCTGAA‑CCCCCCCCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATg                                              <  3:251278/149‑1 (MQ=255)                                                                                                         tAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGAcc                                          >  1:185303/1‑150 (MQ=255)                                                                                                          aaCCCCTCGCGGGTGAACTGCGGGATCACCGCCTGAA‑TCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCa                                         <  2:183671/150‑1 (MQ=255)                                                                                                          aaCACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCa                                         >  2:1235906/1‑150 (MQ=255)                                                                                                            cacCGCGCGGGTAACCGGGGGGATCGTCGCCTGAA‑TCCACCCCTGCGGGCCACCCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCAt                                               <  3:1238864/142‑1 (MQ=255)                                                                                                             acCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAAt                                       >  3:1292051/1‑149 (MQ=255)                                                                                                             acCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAAt                                       >  3:562834/1‑149 (MQ=255)                                                                                                              ccccgcgGGTGAACTGCGGGATCATCGCCTGAA‑TCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGa                                    <  3:875526/148‑1 (MQ=255)                                                                                                                tCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCCACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAtt                                  >  4:456270/1‑151 (MQ=255)                                                                                                                  gggggTGAACTGCGGGATCACCGCCTGAA‑TCTCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTa                                <  3:514638/149‑1 (MQ=255)                                                                                                                  gcgGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAttt                                 >  4:402499/1‑150 (MQ=255)                                                                                                                      gTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCa                             >  1:179930/1‑150 (MQ=255)                                                                                                                         aaCTGCGGGATCACCGCCTGAA‑TCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAAc                          <  1:411917/150‑1 (MQ=255)                                                                                                                          aCTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAAcc                         >  3:356075/1‑150 (MQ=255)                                                                                                                             gCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGc                     >  1:262529/1‑151 (MQ=255)                                                                                                                             gCGGGATCATCGCCTGAA‑CCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGc                     <  1:278670/151‑1 (MQ=255)                                                                                                                               gggATCATCGCCTGAA‑TCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCAt                   <  4:262526/151‑1 (MQ=255)                                                                                                                               gggATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCa                        <  1:394760/146‑1 (MQ=255)                                                                                                                                  atcatcGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGt                >  4:223475/1‑151 (MQ=255)                                                                                                                                    catcGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTg              >  3:1241826/1‑151 (MQ=255)                                                                                                                                     atcGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             >  2:1199547/1‑151 (MQ=255)                                                                                                                                     acctccTGAA‑CCACCCCCCGCGGGCCAAACTCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             <  3:200515/147‑1 (MQ=255)                                                                                                                                      tcGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             >  3:1240373/1‑150 (MQ=255)                                                                                                                                       cGCCTGAA‑TCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGga            <  1:56737/150‑1 (MQ=255)                                                                                                                                       cGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             <  2:969616/149‑1 (MQ=255)                                                                                                                                         ccTGAA‑CCTCCCCCGGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGAt         <  2:143148/151‑1 (MQ=255)                                                                                                                                            gAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGt      >  2:225131/1‑151 (MQ=255)                                                                                                                                                  tcACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTcagc  >  4:278238/1‑150 (MQ=255)                                                                                                                                                    aCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTcagc  >  3:695186/1‑148 (MQ=255)                                                                                                                                                    aCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTcagc  >  3:1016247/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTAT‑GGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAA‑TCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689369‑689661

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