Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	47.5%	86.2 / 16.6	60	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (3/28);  new base G (28/0);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 9.70e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAA  >  NC_000913/925753‑926049                                                                                                                                                     |                                                                                                                                                         gTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTc                                                                                                                                                       <  2:505387/151‑1 (MQ=255)        gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAAcc                                                                                                                                                  <  2:832631/149‑1 (MQ=255)              cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTTAACCCCCCATAc                                                                                                                                          >  2:1021154/1‑151 (MQ=255)                 gCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGtt                                                                                                                                       <  1:547317/151‑1 (MQ=255)                  cAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGAGGGGGGGGGTTTGAAccccccaacgtt                                                                                                                                        >  2:974307/1‑143 (MQ=37)                         gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGtat‑a                                                                                                                              <  2:50733/151‑1 (MQ=255)                             ttCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCCTACGTTTCCGTAT‑aca                                                                                                                            >  1:845275/1‑149 (MQ=255)                             ttCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGGACCCCCCCTACGTTTcccttt                                                                                                                                >  1:665459/1‑142 (MQ=25)                              tCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑acaca                                                                                                                          <  1:370653/150‑1 (MQ=255)                                  aGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTc                                                                                                                     <  1:428606/151‑1 (MQ=255)                                  aGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTc                                                                                                                     <  1:481384/151‑1 (MQ=255)                                    gTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑AGGGGGGGTTTGAACCCCCCCTACGTTGCCGGAA‑ACAAACttt                                                                                                                      >  1:904734/1‑148 (MQ=25)                                      aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAAACCCCGCTACGTTGCCGTAT‑ACACACctt                                                                                                                      >  2:720792/1‑146 (MQ=255)                                           aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCGCTACGTTGCCGTAT‑AAACACTTTTCCGGCggg                                                                                                             >  1:260719/1‑148 (MQ=21)                                               cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCTA‑GGGGGGGGTTTGAACCCCCCCAACGATGCCGTCT‑ACACACTTTTCACGCGCGCttc                                                                                                         >  2:289209/1‑148 (MQ=21)                                                 cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCtt                                                                                                       <  1:45488/150‑1 (MQ=255)                                                  ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCa                                                                                                     <  1:669122/151‑1 (MQ=255)                                                   cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAACCCCCCCACGGTGTCGTAT‑ACACACCTTTCACGGGTGGTcc                                                                                                         >  2:374679/1‑146 (MQ=18)                                                     cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGGACCCCCCATTCCTTTCCGTAT‑ACACACCTTTCAGGGGGGCTCCCTCAGcc                                                                                                  >  1:973267/1‑151 (MQ=18)                                                         ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCGCCACGCTTTCGTTT‑ACACAATTTTCCCGAGTGCCCCTTTAGCCACTc                                                                                              >  1:645340/1‑151 (MQ=18)                                                         ccaTGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTc                                                                                              <  2:315959/150‑1 (MQ=255)                                                                cACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAACCCCCCATACGTTGCCCTAT‑ACACACTTTCCAGGCGGGCTCCTTCACCCACCCGGACA‑cc                                                                                      >  2:363229/1‑151 (MQ=17)                                                                 aCTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCACACGTTGCCGGAT‑ACACACCTTTCAGGCGTGCTCCCTCAGCCACCCCGACA‑ccc                                                                                     >  1:908670/1‑150 (MQ=17)                                                                     ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTTGGACCCCCCCTACGTTGCCGTAT‑ACACACCTTTCCGGCGGGGTCCTTCAACCACCCCGACA‑ccccccc                                                                                 >  1:330461/1‑146 (MQ=17)                                                                      cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑TGGGGGGGTTCGAAACCCCCCCACGCTTCCCCAT‑ATAAACCTTCCGCGCGTGGTTCTTCTGCCACACGGAAA‑CCT‑CAcc                                                                                >  1:949134/1‑150 (MQ=12)                                                                             acagaTAAAATGCGAAAAAAAAGCTCGCACCTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGtt                                                                        <  1:530795/151‑1 (MQ=255)                                                                                gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAAA‑AGGGGGGAATCTGAAACCCGAAAAGTTTTGCTAT‑ACACAAATTCCAAAGGTGCCCCTTTAAACACCCGGCCA‑CCA‑CACCAAAAtttttt                                                                      >  2:77118/1‑150 (MQ=9)                                                                                 aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGc                                                                    <  1:555809/151‑1 (MQ=255)                                                                                  tAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGGTCGAACCCCCCCTACGGTTCCGTTTAAAAAACTTCCCGGGGGGGTCCCTTCGGCCCCT‑GGGCA‑CCT‑CCCCAAAATGttttttt                                                                   >  1:99218/1‑148 (MQ=9)                                                                                          gAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAACCCCCCATAAGTTGCCGCAT‑ACAAACTTTCCCGGGGGGCGCCTTCTGCCCCTCCGACA‑CCT‑CACCCAAATGTTTTTCTTCCAAAcc                                                           >  1:416948/1‑151 (MQ=14)                                                                                            aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTc                                                         <  1:1003955/151‑1 (MQ=255)                                                                                            aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGTTTGAAACCCCCAAAAGTTGACGTAT‑AAACAAATTCCAGGGGGGGCCCTTTCACCACTCGGACA‑CCC‑CACCAAAATTTTTTTTTACCaaaacc                                                          >  2:748119/1‑147 (MQ=2)                                                                                             aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCGATACGGTGCCGTGT‑ACACACTTTTCCGGCGGGGTCCTTCAGCCACTCGGGCACCCC‑CCCCAAATtttttttttcccaacc                                                            >  1:124913/1‑138 (MQ=11)                                                                                               aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATg                                                      <  2:973222/151‑1 (MQ=255)                                                                                                aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTTGAACCCCCCCTACGTTTCCGTTT‑ACACACTTTTCCGGCGGGCGCCTTTAGGCACTCGGGCA‑CCTCCCCCAAATttttttttatccaaa                                                             >  1:227655/1‑135 (MQ=2)                                                                                                 aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGGAACCCCCCGAAGTTTGCGTAT‑ACACACTTTTCAGGCGGGGTCCTTCTGCCCCCCCGGCA‑CCC‑CACCCAAATGTTTTTCTACCAAAACTccttgg                                                    >  2:879415/1‑146 (MQ=9)                                                                                                  aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCCTACGTTGCCGTAT‑ACACACTTTTCAGGCGTGCTCCTTTCGCCCCCCGGGCA‑CCC‑CACCCAAATTTTTTTTTTCCAAACCCCCTggg                                                    >  1:541549/1‑150 (MQ=11)                                                                                                            ttCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCCTACGTTGCCGTAT‑ACACACTTTTCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAAATGTTTTTCTACCAAACCTCATGGGGGGCAACgggg                                         >  1:780332/1‑151 (MQ=17)                                                                                                            ttCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACggg                                          <  2:224176/150‑1 (MQ=255)                                                                                                                 aCGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGcta                                    <  2:105398/151‑1 (MQ=255)                                                                                                                   gAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGGACCCCCCCCACGTTTCCGTAT‑AAAAACTTTTCAGGCGGGGTCCTTCCGCCCCCCCGACC‑CCC‑CACCCAAATGTTTTTCTACCAAAACTCATGGGTGGCAAAGGGGGGcta                                    >  2:283480/1‑149 (MQ=12)                                                                                                                   gAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCCCTACCTTGTCGTAT‑AAACACCTTCCAGGCGGGGTCCTTCCGCCCCCCGGACA‑CCC‑CCCCCAAATGTTTTTCTACCAAACCTCATGGGGGGCGACGGGgcg                                       >  2:277804/1‑146 (MQ=12)                                                                                                                   gAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAAACCCCCATACGTTGCCGTAT‑ACACACTTTTCCGGCGGGGTCCTTCCACCCCCCGGACA‑CCC‑CCCCCAATTGTTTTTCTTCCAAAACCCCTGGGGGGGAAAGGGgcg                                       >  2:576252/1‑146 (MQ=2)                                                                                                                      ctctTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CTCCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACtat                                <  1:941621/150‑1 (MQ=255)                                                                                                                       tcttctTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCATACCTTTCCGTAT‑TCACACTTTTCCGGGGGGCTCCTTCAGCCCCCCGGACA‑CCC‑CACCCAAATGTTTTTCTTCCCAACCCCATGGGGGGGAAAGGGgcgc                                      >  2:210663/1‑143 (MQ=9)                                                                                                                          tcttTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAg                              <  2:322977/148‑1 (MQ=255)                                                                                                                             ttAAATATGGCGGTGA‑GTTTGCTCTTCGATCCCCCTATACGTTGCCTTAT‑ACTCTCTTTCCTGGCGTTCTCCTTCATCCGCTCTTTCA‑CCT‑CACCAAATTGTTTTGCGTCCAAACCTCGTGGGTGTCATCTGCGCACTCCtatatttagt                         <  2:968893/150‑7 (MQ=255)                                                                                                                             ttAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGa                           <  1:622286/148‑1 (MQ=255)                                                                                                                                 atatGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAg                    <  1:657184/151‑1 (MQ=255)                                                                                                                                  tatGGCGGTGA‑GGGGGGGGGTCGAACCCCCCATACGTTGCCGTAT‑ACACACTTTTCCGGCGGGGTCCTTCCGCCACCCGGACA‑CCC‑CCCCCAAATGTTTTGCTACCAAACCCCATGGGGGGGAACGGGGCGCTACTATAAGGGGGTGGAg                    >  2:703548/1‑150 (MQ=18)                                                                                                                                       cGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc              <  2:416038/151‑1 (MQ=255)                                                                                                                                        ggTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc              >  1:344487/1‑150 (MQ=255)                                                                                                                                        ggTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg             <  2:1241/151‑1 (MQ=255)                                                                                                                                          tGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt           >  2:226367/1‑151 (MQ=255)                                                                                                                                            a‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa         <  2:503590/151‑1 (MQ=255)                                                                                                                                              gggggggATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa         <  1:862334/150‑1 (MQ=255)                                                                                                                                                 ggggATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa     <  1:315411/151‑1 (MQ=255)                                                                                                                                                 ggggATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa     <  1:926650/151‑1 (MQ=255)                                                                                                                                                   ggATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAagaaga   <  2:990505/151‑1 (MQ=255)                                                                                                                                                    gATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAgaagaa  <  2:330443/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                         GTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTAT‑ACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCT‑CACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAA  >  NC_000913/925753‑926049

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