Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	42.5%	122.1 / 24.2	74	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (2/40);  new base G (31/0);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.36e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925751‑926050                                                                                                                                                        |                                                                                                                                                            gggTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGAt                                                                                                                                                            <  2:429196/149‑1 (MQ=255)           ccGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCCTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AAccccc                                                                                                                                                 <  1:1042601/151‑1 (MQ=255)              tACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTG‑AAACCCCCAt                                                                                                                                              >  2:141046/1‑151 (MQ=255)               aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTG‑AAccccc                                                                                                                                                 >  1:266345/1‑147 (MQ=255)               aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGAt                                                                                                                                              <  1:836683/150‑1 (MQ=255)                 aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑g                                                                                                                                          <  1:297050/151‑1 (MQ=255)                  ggCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑Gt                                                                                                                                         <  2:900410/151‑1 (MQ=255)                            aaaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATac                                                                                                                               <  2:1090384/151‑1 (MQ=255)                              aTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGGTTG‑AACCCCCCAAAC‑GTTGCCGTATaca                                                                                                                              >  2:339536/1‑150 (MQ=25)                               ttCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGATTG‑AAacccc                                                                                                                                                 >  2:233379/1‑131 (MQ=255)                                 cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATacaca‑c                                                                                                                          <  2:1235066/150‑1 (MQ=255)                                 cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGTTTG‑AACCCCCCCTAC‑GTTGCCGCATaaaca‑c                                                                                                                          >  2:1044214/1‑150 (MQ=25)                                  acaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCG‑AAACCCCCAAAA‑GTTGGCGGATaca                                                                                                                              >  1:696372/1‑146 (MQ=25)                                    aGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGGTCG‑AACCCCCGCTAC‑GTTGCCGTATAAACA‑Cttt                                                                                                                       >  1:271249/1‑150 (MQ=255)                                       tAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTG‑AAACCCCCAAAC‑GATGCCGGAAACACA‑Cctt                                                                                                                       >  1:256869/1‑147 (MQ=25)                                         aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGc                                                                                                                 <  1:1283589/151‑1 (MQ=255)                                         aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCCTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGc                                                                                                                 <  1:202901/151‑1 (MQ=255)                                         aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAAGGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAgtc                                                                                                                 <  1:1013422/151‑3 (MQ=255)                                          cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGCAACCCCCCA‑AC‑GCTGGCCGATTCACC‑CTTTCCAgg                                                                                                                  >  2:253025/1‑149 (MQ=21)                                          cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCg                                                                                                                <  1:659155/151‑1 (MQ=255)                                           aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCG‑AAACCCCCCCAC‑GTTGCCGGATACACA‑CTTTcc                                                                                                                     >  1:450669/1‑145 (MQ=21)                                           aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGt                                                                                                               <  1:525734/151‑1 (MQ=255)                                             aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGc                                                                                                             <  2:955788/151‑1 (MQ=255)                                               tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTG‑AAACCCCCATAC‑GTTGTCGTATACACACCTTTCCAgcgtg                                                                                                               >  1:416081/1‑145 (MQ=21)                                                ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTcc                                                                                                          <  2:172161/151‑1 (MQ=255)                                                  aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCtt                                                                                                        <  1:122031/151‑1 (MQ=255)                                                   cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTG‑AAACCCCCAAAC‑GTTGGCGTATAAAAA‑CTTTCCCgggggggccc                                                                                                          >  2:1132581/1‑140 (MQ=21)                                                   cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTc                                                                                                       <  2:1008079/151‑1 (MQ=255)                                                           ccaTGTTCACTGCCGTACAGACAGATAAAATGCG‑‑AAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTc                                                                                               <  1:605412/150‑1 (MQ=255)                                                            caTGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTG‑AAACCCCCCTAC‑GGTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTTACCCACTc                                                                                               >  1:18213/1‑150 (MQ=18)                                                              tGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCgg                                                                                             <  1:784019/150‑1 (MQ=255)                                                               gTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCG‑AACCCCCCATAC‑GTTGCCGCATACACA‑CTTTTCAGGCGTGCTCCTTCAGCCACCCgg                                                                                             >  2:1163458/1‑149 (MQ=17)                                                               gTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGac                                                                                           <  1:1254145/151‑1 (MQ=255)                                                               gTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGATCG‑GACCCCCCCCAC‑GTTTCCGTATACACA‑CCTTCCAGGGGGGGGCCTTCAccccccccg                                                                                             >  1:295946/1‑143 (MQ=17)                                                                  cACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGATC‑GGGGGGGGCACC‑CACCCCCCACAA‑GTTCCCGGATTCACA‑AATTCCATGCGCGCTCCCTCAGTCAcaccgtca                                                                                          >  2:1005653/1‑142 (MQ=12)                                                                   aCTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCt                                                                                       <  2:1172400/151‑1 (MQ=255)                                                                         gTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCaa                                                                               <  2:1362440/151‑1 (MQ=255)                                                                           acagacagATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCG‑AAACCCCCATAC‑GCTTCCGTATACACA‑CTTTTCAGGGGGGCTCCTTCAGCCACTCGGGCAccccccccaa                                                                                >  2:351799/1‑142 (MQ=14)                                                                             agacagaTAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCG‑AACCCCCCCTAC‑GTTTCCGTATAAACA‑CTTTTCAGGCGGGCTCCTTCAGCCACCCGGGCACCT‑CA‑CCaaa                                                                              >  2:359490/1‑148 (MQ=14)                                                                              gacagaTAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGt                                                                          <  2:275224/151‑1 (MQ=255)                                                                              gacagaTAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGt                                                                          <  2:431972/151‑1 (MQ=255)                                                                               acagaTAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGtt                                                                         <  1:1146256/151‑1 (MQ=255)                                                                                cagaTAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGTTTG‑AACCCCCCACAC‑GTTGCCGTATAAACA‑ATTTCCAGGCGGGCTCCTTCCGCCACCCCGACACCC‑CA‑CCCAAATAttt                                                                        >  1:215299/1‑151 (MQ=14)                                                                                  gaTAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGATA‑GGGGGGGGTTTGTAACCCCCG‑TCA‑GGTGTCGTATTCACC‑CTTTTCAGCCGCGCGCCTTCAGCCACTCCGGAACCC‑CA‑CCAAAATAtttt                                                                       >  1:1225526/1‑150 (MQ=9)                                                                                   aTAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCG‑AACCCCCCCAAC‑GGTGCCGTATACAAA‑CTTTTCCGGCGGGCTCCCTCAGCCACCCGGACACCT‑CA‑CCCAAATGtttt                                                                       >  2:1235840/1‑149 (MQ=14)                                                                                     aaaaTGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGTTGG‑GAACCCCCCTAC‑CTTGGCGTATAAACA‑CTTTTCAGGCGGGGGCCTTCAGCCCCTCCGGCACCTCCC‑CCAAATTTTTTTTCt                                                                    >  1:327460/1‑151 (MQ=9)                                                                                       aaTGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCG‑GACCCCCCATAC‑GTTGCCGTATACACA‑CTTTTCAGGCGTGCTCCTTCAGCCACCCGGACCCCT‑CA‑CCCAAATGTTTTTCTAcc                                                                 >  2:1165087/1‑151 (MQ=14)                                                                                        aTGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATTGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATGC‑GTTGTCGTATACACT‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCa                                                                <  2:539721/151‑1 (MQ=255)                                                                                        aTGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTG‑GACCCCCCACACGGTTGCC‑TATACACA‑CTTTTCACGCGTGGGCCTTCAGCCCCCCCGACACCC‑CA‑CCCAAATGTTTTgtttcc                                                                 >  1:1368802/1‑145 (MQ=9)                                                                                            g‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGG‑GGGGGGGGATCG‑AAACCCCCAAAA‑GCTGCCGGTTAACCC‑CCTTCCCGGGGTGCTCCTTCAACCCCTCCGACACCC‑CA‑CCAAAATGTTTTTGTAACAAAcc                                                            >  2:480195/1‑151 (MQ=2)                                                                                               aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTc                                                          <  2:359879/151‑1 (MQ=255)                                                                                                         cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGc                                                 <  2:573252/150‑1 (MQ=255)                                                                                                          gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGAGGGGGGGGTTTTG‑AAACCCC‑CCCC‑CTTTTCGTAATCACC‑CCTTCCCGGCGCGCTCCCTCCGCCCCCCCGGCACCC‑CC‑CCCAAATGTTTTTCTAACAAACCCCCTGGGTg                                                   >  1:98448/1‑147 (MQ=2)                                                                                                            aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACg                                             <  2:1130016/151‑1 (MQ=255)                                                                                                                tCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTTA‑ACCCCCCAAAAC‑ATTG‑CGGATAACCA‑CTTTTCAGGCCTGCTCCCTCCACCACTCGGAAAACC‑CACCCAAATTTTTTTCTTACAAAACCCCCGGGGGGGCACCgg                                            >  1:101538/1‑149 (MQ=2)                                                                                                                    aCGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCG‑AAACCCCCAAAC‑GTTGCCGGATACACA‑CCTTCCAAGCGTGCTCCCTTAGCCACCCCGACACCC‑CC‑CCAAAATTTTTTTCTTCCAAAACTCATTGGTGGGAACcggggg                                        >  2:640424/1‑146 (MQ=9)                                                                                                                      gAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGctact                                   <  2:522843/151‑1 (MQ=255)                                                                                                                        gCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCG‑AACCCCCCATAC‑GTTGGCGTATACACA‑CTTTTCAGGCGTGGTCCTTCAGCCACTCGGGCACCT‑CA‑CCAAAATTTTTTTCTTCCCAACCTCCTTGGGGGCAACgggg                                          >  2:1121357/1‑142 (MQ=14)                                                                                                                        gCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCG‑AACCCCCCATAC‑GGTTCCGTATAAACA‑CTTTTCAGGCGGGCTCCTTCAACCCCTCGGACACCC‑CA‑CCAAAATTTTTTTTTACCAAAACTCCTGGGGGGCAACGGGGGGCTACtat                                 >  2:1322499/1‑151 (MQ=14)                                                                                                                         ctctTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTG‑AAACCCCCCTAC‑GTTTCCGGATACACA‑CTTTCCCGGGGTGCTCCCTCAGCCACTCGGACACCC‑CA‑CCCAAATGTTTTTCTTCCAAACCTCCTGGGGGGGAACGGGGGGcta                                     >  2:1004143/1‑146 (MQ=9)                                                                                                                             tcttTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAggg                             <  2:272760/150‑1 (MQ=255)                                                                                                                                 tAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                         <  2:834317/150‑1 (MQ=255)                                                                                                                                 tAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                         <  1:434092/150‑1 (MQ=255)                                                                                                                                   aaTATGGCGGTGA‑GGGGGGGGGTCG‑AACCCCCGATAA‑GTTTCCGTATAAACA‑CTTTTCAGGCGTGGGTCCTCCGCCACTCGGGCACCT‑CA‑CCCAAATTTTTTTCTACCCAAACCCATTGGGGGGAACCGGGGGGTTCTATATGGAggt                         >  1:760953/1‑146 (MQ=2)                                                                                                                                      atGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                   <  1:551753/151‑1 (MQ=255)                                                                                                                                        ggCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaa                  >  2:741759/1‑150 (MQ=255)                                                                                                                                            gTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACgg             <  2:142873/151‑1 (MQ=255)                                                                                                                                               a‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc           >  1:542583/1‑150 (MQ=255)                                                                                                                                               a‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          <  2:1042503/151‑1 (MQ=255)                                                                                                                                                  ggggggATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaag        <  2:752561/151‑1 (MQ=255)                                                                                                                                                   gggggATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga       <  2:1023702/151‑1 (MQ=255)                                                                                                                                                    ggggATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa      <  2:570628/151‑1 (MQ=255)                                                                                                                                                    ggggATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa      <  1:339535/151‑1 (MQ=255)                                                                                                                                                    ggggATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa      <  2:902935/151‑1 (MQ=255)                                                                                                                                                        aTTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  2:36078/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                            GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCG‑AACCCCCGATAC‑GTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCT‑CA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925751‑926050

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