Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	44.4%	115.9 / 36.2	81	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (10/35);  new base G (36/0);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 3.96e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGA  >  NC_000913/925754‑926048                                                                                                                                                    |                                                                                                                                                           ttCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGT‑GGGGGGGGtttg                                                                                                                                                        >  1:485547/1‑149 (MQ=255)    aCTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AAc                                                                                                                                                   <  2:43214/151‑1 (MQ=255)       gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AAcccc                                                                                                                                                <  2:1303280/151‑1 (MQ=255)        ccGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTGT‑G‑AAccccc                                                                                                                                               >  2:9748/1‑151 (MQ=255)          gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGa                                                                                                                                             <  2:529876/151‑1 (MQ=255)          gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGGTG‑A‑AAccccc                                                                                                                                               >  2:395280/1‑149 (MQ=255)             cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATAc                                                                                                                                          <  2:1302615/151‑1 (MQ=255)                  aGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTT‑G‑AAccccccccacgtt                                                                                                                                       >  1:767132/1‑141 (MQ=37)                  aGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGATA‑TGGGGGGGTTT‑G‑AAACCCCGagtccgttc                                                                                                                                     >  2:1459682/1‑143 (MQ=37)                    cTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTT‑G‑AACCCCCCCTACGTTTcc                                                                                                                                    >  2:1012615/1‑150 (MQ=37)                    cTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTT‑G‑GAccccccccaccgttcc                                                                                                                                    >  1:397375/1‑139 (MQ=37)                      tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTC‑G‑AAaccccccaaactt                                                                                                                                       >  1:1600122/1‑137 (MQ=37)                      tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCg                                                                                                                                   <  2:1452529/149‑1 (MQ=255)                       aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGtat                                                                                                                                <  1:664134/151‑1 (MQ=255)                       aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGATG‑GGGGGGGGGTC‑G‑GACCCCCCACACgcttcc                                                                                                                                    >  1:37570/1‑142 (MQ=37)                         aaaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTT‑G‑AACCCCCCCCACGTTGCCGtgt                                                                                                                                >  1:705135/1‑147 (MQ=37)                                        aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGt                                                                                                              <  2:792195/151‑1 (MQ=255)                                         tACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCg                                                                                                               <  1:255896/149‑1 (MQ=255)                                          aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTT‑G‑AACCCCCCCCACGTTGCCGGATAAAAA‑CTTTTCAGGCggg                                                                                                             >  2:1337020/1‑148 (MQ=21)                                          aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTg                                                                                                             <  1:1345665/150‑1 (MQ=255)                                             ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTcc                                                                                                         <  1:868363/151‑1 (MQ=255)                                              cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCt                                                                                                        <  1:534541/151‑1 (MQ=255)                                                 ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTC‑G‑CAACCCCCATACCTTTCCCCATACAAA‑CCTTCCacg                                                                                                                 >  1:555785/1‑137 (MQ=21)                                                 ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTC‑G‑AACCCCCCACACGTTGCCGCATACACA‑CCTTTCCGGCGTGGTCCtt                                                                                                       >  1:158617/1‑149 (MQ=18)                                                      caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCAc                                                                                                <  2:1671306/151‑1 (MQ=255)                                                      caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCCCTACGTTGCCGTATACACA‑CTTTCCAGGCGGGCTCCCTCAGcc                                                                                                  >  2:529467/1‑149 (MQ=255)                                                        ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTc                                                                                              <  1:1238042/151‑1 (MQ=255)                                                            gTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGac                                                                                          <  2:232462/151‑1 (MQ=255)                                                                    ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTT‑G‑AACCCCCCAAACGTTTCCGGATACACA‑CATTTCCCGCGGGGTCCCTCTGCCACCCCGGCA‑ccccccc                                                                                 >  2:1496060/1‑146 (MQ=12)                                                                    ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGTTT‑G‑AAACCCCCATACGTTGCCGTAAAAACA‑CTTTCCAGGCGCGCTCCTTCCACCACTCGGACA‑CCCCaac                                                                                 >  2:687849/1‑149 (MQ=17)                                                                     cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCa                                                                               <  2:892274/151‑1 (MQ=255)                                                                      gTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGTTT‑G‑GAACCCCCCTACGTTTCCGTATAAACA‑CTTTCCCGGCGGGCTCCTTCCGCCCCTCCGACA‑CCCCA‑c                                                                                 >  1:257673/1‑148 (MQ=17)                                                                         cagacagcTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑AGGGGGGAATT‑G‑AACACCCGCCACGATTCCGTATAAAAA‑CTTTTCAAACGTGCGCCTTTCACCACTCGCACA‑CCGCA‑CCaaa                                                                             >  1:320871/1‑149 (MQ=9)                                                                         cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAAtt                                                                           <  1:1164057/151‑1 (MQ=255)                                                                          agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGATC‑G‑GACCCCCCCTAAGTTTCCGTAAAAAAA‑ACTTTCAAGCGGGCCCCTTCACCCACCCCGGCA‑CCCCC‑CCaaa                                                                             >  1:758186/1‑148 (MQ=9)                                                                           gacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTC‑G‑AACCCCCCCTACGTTGCCGCATACACA‑CTTTCCACGCGGGCTCCTTCCGCCCCCCCGGCA‑CCCCA‑CCcaa                                                                             >  2:407012/1‑147 (MQ=14)                                                                             cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTC‑GCAACCCCCG‑TAAGTTTCCGCATACACA‑CTTTTCCGGCGTGCGCCCTCCCCCCCCCGGGCA‑CCCCA‑CCCAAATGttt                                                                       >  2:270030/1‑151 (MQ=9)                                                                               gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTg                                                                     <  1:397258/151‑1 (MQ=255)                                                                                aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑AGGGGGGGTTC‑G‑AAACCCCCAAACGTTTCCGTATACAAA‑CCTTCCAGGGGTGCTCCCTCAGCCAACCGGACA‑CCTCA‑CCACAATGtttt                                                                      >  1:680914/1‑149 (MQ=14)                                                                                  aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTT‑G‑AACCCCCCCTACGTTGCCGTATACACA‑CCTTCCCGGCGTGGTCCTTTCTCCACCCCGGCA‑CCCCC‑CCAAAATTTTTTTCt                                                                   >  1:531666/1‑150 (MQ=14)                                                                                  aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGTTT‑G‑AACCCCCCGCAAGTTGGCGTATTAACACCTTTCCAGG‑GTGCTCCCTCCGCCCCCCCGGCA‑CCCCC‑CCCAAATGtttttttt                                                                  >  1:574269/1‑147 (MQ=9)                                                                                  aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTT‑G‑GACCCCCGAAAAGATTCCGGATACACA‑CTTTTCCGGGGTGGTCCCTCCGCCACCCCGACA‑CCCCC‑CCCAAATGTTTTTCtt                                                                  >  1:1584590/1‑150 (MQ=9)                                                                                   aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTT‑G‑AAACCCCCATAACGTGCCGTATACAAA‑CTTTTCAAGCGTGCTCCCTCCGCCACACCGAAA‑ACCCC‑CCCAAATGTTTTGCt                                                                   >  2:1136481/1‑149 (MQ=9)                                                                                   aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGATT‑G‑AAACCCCCCTACGCTGCCGTATAAACA‑ATTTCCAGGCGGGCTCCTTCAGCCCCTCCGACA‑CCCCA‑CCAAAATGttttttta                                                                  >  1:125206/1‑146 (MQ=14)                                                                                    aaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTT‑G‑AAACCCCCCTAAGTTGCCGTATAAAAA‑CTTTTCAAGCGGGGTCCCTCACCCCCTCCGGCA‑CCCCA‑CCCAAAtttttt                                                                      >  1:1453060/1‑145 (MQ=2)                                                                                     aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTC‑G‑GAACCCCCAGACGTTGCCGTCTATACACCTTTCCGGGGGGGTCCCTCCACCCAC‑CGGGCA‑CCCCACCCAAATtttttttttccc                                                                >  1:1341247/1‑145 (MQ=2)                                                                                       gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑aaa                                                            <  1:1690511/151‑1 (MQ=255)                                                                                         gAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAAcc                                                          <  1:757141/151‑1 (MQ=255)                                                                                          aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCt                                                         <  2:220664/151‑1 (MQ=255)                                                                                           aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTc                                                        <  2:1405759/151‑1 (MQ=255)                                                                                           aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAAAGGGGGGGGATCG‑G‑AACCCCGGCAACGTTGCC‑TATACACA‑CTTTCCCCGGGGGCGCCCTCAGCCACACCGACA‑CCCCC‑CCAAAATGTTTTTTTTCC‑aaa                                                            >  1:1217921/1‑147 (MQ=2)                                                                                             aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTA‑G‑AAACCCCCCTACGTTGCCGTATACAAA‑CCTTTCTGGCGGGCCCCTTCTTCCACCCCGGCACCCCCC‑CCAAATttttttcttcc                                                                 >  2:835790/1‑137 (MQ=2)                                                                                               aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTT‑G‑ACCCCCCG‑GAAGTTTCCGTGTAACAA‑ATTTTCAGGGGGGCTCCTCTAACCACCCGCGCA‑CCCCA‑CCACAATGTTTTTTTACC‑CAAACTc                                                        >  1:1609802/1‑147 (MQ=2)                                                                                                aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGG‑GGGGGGGGGTC‑A‑AACCCCCCAAATGGTGTCCTATGATCA‑AATTTCCGGCGGGCTCCTTCACCCCCCCCCCCA‑CCCCA‑CCCAACTGTTTTTTTACC‑AAACCCCaaagg                                                   >  1:841856/1‑147 (MQ=9)                                                                                                     cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCa                                              <  2:1274953/151‑1 (MQ=255)                                                                                                      gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCaa                                             <  2:3475/151‑1 (MQ=255)                                                                                                             cGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTT‑G‑AAACCCCCCCACGTTTCCGTATAAACA‑CTTTCCCCGCGGGCTCCCTCCCCCACCCCGGCA‑CCCCC‑CCCAAATGTTTTTTTTCC‑CAACCCCATGGGGGGGAACgggg                                        >  1:1168029/1‑149 (MQ=2)                                                                                                             cGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGgcg                                      <  2:962883/151‑1 (MQ=255)                                                                                                              gTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGGTCTG‑AACCCCC‑CCACGTTTTCCTATACACACCTTTCCCGCGGGGCCCCTTC‑GCCCCACCGGCA‑CCCCA‑CCACAATTTTTTTTTACC‑CAAACCCATGGGGGGGAAAgggg                                        >  1:829082/1‑148 (MQ=2)                                                                                                               tACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTT‑G‑AAACCCCCCTACATTGCCGGATAAAAA‑CATTCCCAGCGCGATCCCTCCACCACACCGGCA‑CCCCA‑CCCAAATGTTTTTCTACC‑AAACCCCCAGGGTGGGAAAgggg                                        >  2:1106767/1‑147 (MQ=9)                                                                                                               tACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTC‑G‑AAACCCCCACACGTTTCCGTATACACA‑CCTTCCAGGGGGGGTCCCTCCACCCCCCCGGCA‑CCCCA‑CCCAAATGTTTTTTTACC‑AAAACTCATTGGGGGGCACgggg                                        >  1:463710/1‑147 (MQ=2)                                                                                                                 cGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTC‑G‑AACCCCCCCTACGTTTCCGCATAAACA‑CCTTTCACGCGGGCTCCCTCCGCCCCCCCGGCA‑CCCCC‑CCCAAATGTTTTTCTACC‑AAAACTCATGGGGGGGAACGGGGGGccac                                  >  1:1658905/1‑151 (MQ=9)                                                                                                                   aGCTCTTCTTTAAATATGGCGGCGA‑GGGGGGGATTT‑G‑AACCCCCCATACGTTTTCGTATACACA‑CTTTCCCGGCGGGCTCCTTCAGCCCCCCCGGCA‑CCCCA‑CCCACATTTTTTGCTATC‑CAACCTCCAAGGTGGCGAAgggg                                        >  1:548140/1‑143 (MQ=2)                                                                                                                      tcttctTTAAATATGGCGGTGA‑GGGGGGGAGTT‑G‑GACCCCCCCAAAGTTGGCGTATAAACA‑CTTTCCCGGCGGGGCCCCTTAGCCACTCCGACA‑CCTCA‑CCAAAATGTTTTTTTACCAAAACCCCTGGGGGGGCAAGgggggcca                                    >  2:1081098/1‑141 (MQ=2)                                                                                                                         tcttTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAggg                           <  2:1648379/150‑1 (MQ=255)                                                                                                                            ttAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGt                        <  2:1172421/150‑1 (MQ=255)                                                                                                                              aaaTATGGCGGTGA‑GGGGGGGGGTC‑G‑AACCCCCCACACGTTGCCGTATAAACA‑ATTTTCAGGCGGGCTCCTTCAACCCCCCCGACA‑CCCCA‑CCAAAATGTTTTTCTACC‑AAACCCCCAGGGGGGGAACGGGGGGCTCCTATAGGGGGt                        >  1:312472/1‑148 (MQ=9)                                                                                                                                atatGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTgg                     >  2:1024520/1‑149 (MQ=255)                                                                                                                                  atGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGt                  <  2:1522609/150‑1 (MQ=255)                                                                                                                                  atGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGt                  <  2:1149022/150‑1 (MQ=255)                                                                                                                                  atGGCGGCGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTTCCGTATAAACA‑CCTTCCAGGCGGGCTCCTTCAGCCACCCGGGCA‑CCTCA‑CCAAATTTTTTTGCTACC‑AAACCTCATGGGGGGCAACGGGGCGCCACTATAGGGAggt                       >  1:105481/1‑143 (MQ=18)                                                                                                                                   tGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                 >  2:646499/1‑150 (MQ=255)                                                                                                                                    ggCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaa               <  1:630561/151‑1 (MQ=255)                                                                                                                                    ggCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaa               <  2:531683/151‑1 (MQ=255)                                                                                                                                    ggCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaa                >  2:1463738/1‑150 (MQ=255)                                                                                                                                                ggggATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa    <  1:1511974/151‑1 (MQ=255)                                                                                                                                                ggggATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa    <  2:1381812/151‑1 (MQ=255)                                                                                                                                                 gggATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag   <  1:254250/151‑1 (MQ=255)                                                                                                                                                 gggATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag   <  1:1545827/151‑1 (MQ=255)                                                                                                                                                  ggATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga     >  1:434523/1‑148 (MQ=255)                                                                                                                                                   gATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAagaaga  <  1:529450/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                           TTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACC‑AAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGA  >  NC_000913/925754‑926048

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