Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	45.4%	150.8 / 29.7	97	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (8/45);  new base G (44/0);  total (52/45)
Fisher's exact test for biased strand distribution p-value = 9.11e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAG  >  NC_000913/925751‑926047                                                                                                                                                       |                                                                                                                                                        gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGAt                                                                                                                                                        <  1:973057/151‑1 (MQ=255)      cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGaa                                                                                                                                                   <  2:1748189/151‑1 (MQ=255)      cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGa                                                                                                                                                    <  2:10191/150‑1 (MQ=255)        cTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCGTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGtttg                                                                                                                                                     >  2:87329/1‑145 (MQ=255)        cTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGaagc                                                                                                                                                 >  1:504003/1‑149 (MQ=255)               aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATa                                                                                                                                          <  1:1073888/151‑1 (MQ=255)                cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATAc                                                                                                                                         <  1:314871/151‑1 (MQ=255)                     aGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGAAAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGAGTTGAAcaccccaacagtt                                                                                                                                      >  2:7856/1‑141 (MQ=37)                           gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAACCCCCGCTACGTTTCCGtat                                                                                                                               >  2:739901/1‑150 (MQ=255)                             aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATaca                                                                                                                            <  1:965556/151‑1 (MQ=255)                                tCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAAACCCCCATACCGTGcc                                                                                                                                   >  2:1546345/1‑141 (MQ=37)                                 cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTTGAAACCCCCAAACGtttgcggataaa                                                                                                                            >  1:1728750/1‑137 (MQ=25)                                 cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCCTACGTTGCCGGATaaac‑a‑c                                                                                                                       >  1:1188932/1‑150 (MQ=25)                                  acaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑Ctt                                                                                                                     <  1:1019722/151‑1 (MQ=255)                                   caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGATCGAACCCCCCCTACGATTCCGGTTACAC‑A‑Attt                                                                                                                    >  1:1642918/1‑151 (MQ=25)                                     ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTcc                                                                                                                  <  1:92214/151‑1 (MQ=255)                                     ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTcc                                                                                                                  <  1:705376/151‑1 (MQ=255)                                      gTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCa                                                                                                                 <  2:1212285/151‑1 (MQ=255)                                         aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGc                                                                                                              <  1:1219907/151‑1 (MQ=255)                                         aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGc                                                                                                              <  1:688278/151‑1 (MQ=255)                                         aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGc                                                                                                              <  1:629522/151‑1 (MQ=255)                                          cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCAAACGTTGCCGTATAAAC‑A‑CTTTCCCGGCg                                                                                                             >  2:1541705/1‑151 (MQ=255)                                          cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCg                                                                                                             <  2:37482/151‑1 (MQ=255)                                           aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGt                                                                                                            <  2:1420632/151‑1 (MQ=255)                                            tACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCATACGTTGCCGTATACAC‑A‑CTTTCCAGGCggg                                                                                                           >  1:1523736/1‑149 (MQ=255)                                             aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAAACCCCCCTACGTTTCCGTATACAC‑A‑CCTTCCCGGCGTg                                                                                                           >  1:600861/1‑150 (MQ=21)                                                  aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCCTACCTTGCCGGATACAC‑A‑CTTTTCAGGCGTGGTcc                                                                                                       >  2:215800/1‑149 (MQ=21)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAAACCCCCATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTccct                                                                                                     >  2:651167/1‑147 (MQ=255)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCa                                                                                                   <  1:1242349/151‑1 (MQ=255)                                                     cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCGCTACGTTGGCGGATACAA‑A‑CTTTCCACGCGGGGTCCtt                                                                                                     >  1:271786/1‑148 (MQ=18)                                                        ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAAACCCCCCTACGTTGCCGTATACAC‑A‑CTTTCCCGGCGGGCTCCTTCCGcc                                                                                                >  1:478862/1‑150 (MQ=18)                                                         caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAACCCCCCACACGTTTGCGTATACAC‑A‑CTTTTCAGGCGGGCTCCTTCCGCCAc                                                                                              >  2:161351/1‑151 (MQ=18)                                                          accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGGTCGAACCCCCCATACGTTGCCGTATACAC‑ACCTTTCCA‑GCGGGCTCCTTCCGCCAc                                                                                              >  1:1302631/1‑150 (MQ=18)                                                          accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGGAACCCCCATACGGTGGCCTATACACAA‑CTTTCC‑GGGGGGGTCCTTCAccccc                                                                                              >  2:182301/1‑148 (MQ=18)                                                          accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACt                                                                                             <  1:698825/151‑1 (MQ=255)                                                           ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTc                                                                                            <  2:718113/151‑1 (MQ=255)                                                             aTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGGAACCCCCCTACCTTGCCGTATACAC‑A‑CTTTTCAGGCGTGGTCCtttcgcccccccg                                                                                          >  2:790325/1‑140 (MQ=18)                                                               gTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAAACCCCGATACGTTTCCGGATAGAA‑A‑CCTTTCCGGGGGGCTCCtt                                                                                                     >  1:905057/1‑138 (MQ=18)                                                               gTTCACTGCC‑TACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGac                                                                                        <  1:1434361/150‑1 (MQ=255)                                                                 tCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGGTTGAAACCCCCCTACGTTTCCGTATACAC‑A‑CTTTTCCGGCGGGCTCCTTCAGcc                                                                                                >  2:1563427/1‑141 (MQ=17)                                                                 tCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTCGAAACCCCCCCACGTTGTCCTGTACAC‑A‑CTTTTCCGGCGTGCTCCCTccacccccccggca                                                                                       >  2:396591/1‑137 (MQ=17)                                                                  cACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGaca‑c                                                                                     <  1:1656291/150‑1 (MQ=255)                                                                   aCTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTTGAACCCCCCCCACGTTGCCGTATACAC‑A‑CTTTCCAGGGGGGCGCCTTCCGcccc                                                                                              >  2:592858/1‑139 (MQ=17)                                                                      gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCCTACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACCCCGACA‑cccccc                                                                                >  1:307018/1‑147 (MQ=17)                                                                         gTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAACCCCCCCTACGTTGCCGGATACAC‑A‑CCTTTCCGGCGGGCTCCTTCAACCCCCCGGGCA‑CCCCA‑cc                                                                              >  1:964526/1‑149 (MQ=17)                                                                         gTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTTGAAACCCCCATACGTTGCCGGATACAC‑A‑CTTTTCAGGCGGGCTCCCTCAGCCCCCCGGGCA‑ccccccc                                                                               >  2:1377180/1‑144 (MQ=17)                                                                            cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAAACCCCCATACGGTGCCGGATACAC‑A‑CTTTTCCGGCGTGCTCCTTCAGCCACCCGGACA‑CCCCA‑CCaaa                                                                           >  1:1765859/1‑149 (MQ=14)                                                                              gacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCGATACGTTGCCTTATAAAA‑A‑CTTTCCCGGCGGGCTCCCTCAGCCACCCGGACC‑CCCCACCCAAAttt                                                                        >  2:72410/1‑150 (MQ=14)                                                                               acagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGtt                                                                      <  1:739934/151‑1 (MQ=255)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCCCCACGTTGCCGTATAAAA‑A‑CTTTCCAGGCGGGGTCCTTCAGCCCCCCCGACA‑CCCCA‑CCCAAATGtttt                                                                    >  2:1704001/1‑151 (MQ=14)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAAACCCCCCTACGTTTCCGTATACAC‑A‑CATTCCAGGCGGGGTCCCTCCGCCACCCCGACA‑CCCCA‑CCCAAATGtttt                                                                    >  2:1242373/1‑151 (MQ=14)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGtttt                                                                    <  1:660691/151‑1 (MQ=255)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGtttt                                                                    <  2:574674/151‑1 (MQ=255)                                                                                   aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGc                                                                  <  1:345279/151‑1 (MQ=255)                                                                                              aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTc                                                       <  2:736758/151‑1 (MQ=255)                                                                                               aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTCAACCCCCCATAACTTGCCCTATAACC‑A‑CTTTTCCGGCGGGGTCCTTCAGCCACCCGGGCACCCCCC‑CCAAATTTTTTTCTTCCCAAAc                                                          >  2:306707/1‑148 (MQ=2)                                                                                                 aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATg                                                    <  1:1243222/151‑1 (MQ=255)                                                                                                 aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTCGAACCCCCCCTTAGTTTGCGTATAAAA‑A‑CTTTTCCGGCGGGGGCCCTCAGCCACCCCGGCC‑CCCCC‑CCCAATTGTTTTTTTACCCAACCTc                                                       >  2:247834/1‑148 (MQ=9)                                                                                                  aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCCTACGGTGCCGGATAAAC‑A‑CTTTCCAGGCGGGCTCCCTCAGCCACTCCGACA‑CCCCA‑CCCAAATGTTTTTTTACCAAAACTCCTg                                                    >  1:1650216/1‑150 (MQ=14)                                                                                                  aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATgg                                                   <  2:656574/151‑1 (MQ=255)                                                                                                   aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAAACCCCCATACGGTGCCGGATTAAA‑A‑CCTTTCCGGGGGGCTCCTTCAGCCACCCGGGGC‑CCTCA‑CCAAAATTTTTTTTTACCCAACCTccacgg                                                  >  1:466639/1‑146 (MQ=2)                                                                                                   aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAAACCCCCATAAGGTTTCGGATACAC‑A‑CCTTTCAGGCGGGGTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAAATGTTTTTTTACCCAACCCCCTggg                                                  >  1:1025681/1‑151 (MQ=9)                                                                                                    aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGt                                                 <  1:899457/151‑1 (MQ=255)                                                                                                     gCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTCAAACCCCGCTAAGTTGGCGGTTACAA‑A‑CTTTTCCGGCGTGGTCCTTCAGCCACCCGGGCA‑CCTCA‑CCAAATTGTTTTTGTACCAAACCTccttgggg                                                >  1:1579707/1‑144 (MQ=9)                                                                                                     gCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTg                                                <  1:898436/151‑1 (MQ=255)                                                                                                         gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGT‑GGGGGGGGGTTGAAACCCCCGTACGTTGTCTTATAAAC‑A‑CCTTCCAGGCGGGGGCCCTCACCCCCCCGCGCA‑CCCCA‑CCAAAATGTTTTGTTACCaaaacccactggggg                                               >  2:1399843/1‑136 (MQ=2)                                                                                                         gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAAACCCCCACACGTTGCCGTATACAC‑A‑CCTTTCCGGCGGGCTCCCTCCGCCCCCCCGGCA‑CCCCC‑CCAAAATGTTTTTCTTCCCAAcccccttggggg                                               >  1:1094226/1‑138 (MQ=2)                                                                                                          cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAAc                                           <  1:733622/151‑1 (MQ=255)                                                                                                             tttCGTACGAGCTCTTCTTTAAATATGGCGGGGG‑GGGGGGGATTCTAACACCCCCTACCTTGCCGTTTGCAA‑A‑CATTCCAGGAGTGGTCCTTCCCCCCCCCGGAGA‑CACCA‑CCAAAATTTTTTGTGTCCAAAcccccttgggg                                                >  2:364385/1‑134 (MQ=2)                                                                                                              ttCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTGGAACCCCCCCTACCTTGCCGTATAAAC‑A‑CTTTTCCGGCGGGGTCCTTCCGCCCCCCCGGCA‑CCCCC‑CCCAAATGTTTTTCTTCCAAACCCCCTGGGGGGGAAAgggg                                       >  2:1392252/1‑151 (MQ=2)                                                                                                                cGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCATACGTTTGCGGATAAAC‑A‑CTTTCCAGGCGGGGTCCCTCAACCACCCCGGCA‑CCCCC‑CCCAAATTTTTTTTTACCAAAACTCATGGGTGGGAACgggg                                       >  2:1453767/1‑149 (MQ=9)                                                                                                                cGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCTGACCCCCCCTACGTTGCCGTATACCA‑A‑CCTTTCAGGGGGGGTCCTTCAGCCACCCCGACA‑CCCCC‑CCAAAATGTTTTTTTACCAAAACCCATGGGGGGGGACgggggg                                     >  2:1371023/1‑149 (MQ=2)                                                                                                                 gTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAAACCCCCATACGTTGGCGGATACAC‑A‑CTTTCCAGGCGGGCTCCCTCAGCCACTCGGGCA‑CCCCA‑CCCAAATGTTTTTTTTCCAAACCCCATGGGTGGCAACgggggg                                     >  1:1608374/1‑148 (MQ=14)                                                                                                                 gTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGATTCCAACCCCCGCTACGTTGGCGTTTACAC‑A‑CTTTTCAAGCGCGCCCCTTTTACCACTCGGAAA‑CCACA‑CCACAATGTTTTTTTAACCAACCCCTACGGTGGCAACGGGgcg                                     >  1:797933/1‑150 (MQ=2)                                                                                                                    cGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATCCGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGcta                                  <  1:907445/150‑1 (MQ=255)                                                                                                                     gAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGctact                                <  2:596043/151‑1 (MQ=255)                                                                                                                      aGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCATACGTTGCCGGATACAC‑A‑CTTTCCCGGCGGGCTCCTTCAGCCACCCGGACA‑CCCCA‑CCAAAATGTTTTTCTTCCAAACCCCCTGGGGGGGAACGGGGGGctac                                 >  1:764133/1‑149 (MQ=14)                                                                                                                         tcttctTTAAATATGGCGGTGAGGGGGGGGGTTGGAACCCCC‑CTACGGTGCCGTATACAC‑A‑CCTTTCAGGGGGGGCCCTTCCGCCACCCGGGCA‑CCCCA‑CCAAAATGTTTTTCTACCAAACCTCATAGGTGGCAACGGGGCGCTACtat                              >  1:919003/1‑149 (MQ=16)                                                                                                                         tcttctTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAACCCCCATACGTTGCCGTATACAC‑A‑CTTTCCCGGCGGGCTCCTTTAGCCACCCGGACA‑CCCCA‑CCCAAATGTTTTTCTACCCAACCCCCTGGGGGGGAACGGGGGGCTACtat                              >  1:992035/1‑149 (MQ=14)                                                                                                                            tcttTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAg                            <  2:406726/148‑1 (MQ=255)                                                                                                                              tttAAATATGGCGGTGA‑GGGGGGGGGTCGAACCCCCCATAAGTTGCCGTATACAA‑A‑CTCTCCAGGCGGGCTCCCTCCGCCCCCCCGACA‑CCCCC‑CCCAAATTTTTTTCTACCAAACCCCCTGGGGGGCAACGGGGCGCCCCTATAggggg                        >  2:960234/1‑148 (MQ=12)                                                                                                                               ttAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAg                        <  2:1265924/149‑1 (MQ=255)                                                                                                                               ttAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGt                       <  2:8514/150‑1 (MQ=255)                                                                                                                                tAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                      <  1:758233/150‑1 (MQ=255)                                                                                                                                   atatGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTgg                    >  1:1354396/1‑149 (MQ=255)                                                                                                                                   atatGGCGGTGA‑‑GGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGt                 <  2:1029493/151‑1 (MQ=255)                                                                                                                                       ggCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaa              <  1:896062/151‑1 (MQ=255)                                                                                                                                         cGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc            <  1:449614/151‑1 (MQ=255)                                                                                                                                         cGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc            >  2:1433749/1‑151 (MQ=255)                                                                                                                                          ggTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTTA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg           <  2:1642686/151‑1 (MQ=255)                                                                                                                                          ggTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg           <  1:431209/151‑1 (MQ=255)                                                                                                                                          ggTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg           >  2:959076/1‑151 (MQ=255)                                                                                                                                             gA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt         <  1:1692685/150‑1 (MQ=255)                                                                                                                                                   ggggATTCGAACCCCCGATACGTTGCCGTATACAC‑C‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa   <  1:1442614/151‑1 (MQ=255)                                                                                                                                                   ggggATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa   <  2:503982/151‑1 (MQ=255)                                                                                                                                                    gggATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag  <  1:1532350/151‑1 (MQ=255)                                                                                                                                                     ggATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga    >  1:1279095/1‑148 (MQ=255)                                                                                                                                                       |                                                                                                                                                        GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAG  >  NC_000913/925751‑926047

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: ‑