Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	44.8%	88.3 / 18.1	58	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (1/31);  new base G (26/0);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.22e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925754‑926050                                                                                                                                                    |                                                                                                                                                        ttCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTTA‑GGGGGGGGtt                                                                                                                                                       >  2:169346/1‑149 (MQ=255)        ccGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGGTTGAAacccc                                                                                                                                              >  1:1097177/1‑151 (MQ=255)          gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGa                                                                                                                                            <  2:471676/151‑1 (MQ=255)                             tcgcgggTAGCATACTGCACCCGCCGTCCATGTTCACTGCCGTACAGACAGAGAAAATGCGAAAAAAAAGCGCGCACTTTCGTACGAGCTCTTCTTTATATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGTTGCGTTGCCGtata                                                                                                                              <  1:340909/141‑1 (MQ=18)                                caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACttt                                                                                                                      <  1:462476/151‑1 (MQ=255)                                   gTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCa                                                                                                                   <  1:1146923/151‑1 (MQ=255)                                      aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGc                                                                                                                <  2:882465/151‑1 (MQ=255)                                       cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCg                                                                                                               <  2:1121452/151‑1 (MQ=255)                                         tACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTCGGACCCCCCATACGGTGCCGCATAAACACctt                                                                                                                      >  2:452238/1‑142 (MQ=21)                                             ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTcc                                                                                                         <  2:556958/151‑1 (MQ=255)                                                cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTc                                                                                                      <  1:679686/151‑1 (MQ=255)                                                 ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGGAACCCCCGTAACAAGGCCG‑ATACACACTTTTCAGGCGTGCTc                                                                                                          >  2:1149775/1‑146 (MQ=21)                                                 ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCGATACGTTGCCGGATACACACATTCCAGGCGTGGTCCtt                                                                                                       >  2:138498/1‑149 (MQ=18)                                                  cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGGACCCCCCATACGTTTCCGTATACACACTTTTCAAGCGGGCTccct                                                                                                       >  2:273719/1‑146 (MQ=18)                                                   gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGGAACCCCC‑CTAAGGTGCCGCATACACACTTTTCCGGGGTGCTcc                                                                                                         >  2:610221/1‑145 (MQ=18)                                                     ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGGTTGAAACCCCCCAACGGTGCCGGATACACACTTTTCAGGCGGGGTCCTTCCGcc                                                                                                  >  2:259367/1‑150 (MQ=18)                                                     ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGTTTGAAGCCCCCCAACGTTTCCGGGTAAACACCTTTCCGGTGTGCTCCTTccagc                                                                                                  >  1:416973/1‑146 (MQ=18)                                                       accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACt                                                                                               <  2:860256/151‑1 (MQ=255)                                                            gTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGac                                                                                          <  1:298744/151‑1 (MQ=255)                                                            gTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGac                                                                                          <  1:228152/151‑1 (MQ=255)                                                             ttCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCGATACGTTGCCCTTTACACACCTTCCACGCGTGCTCCCTccaccccccgggc                                                                                          >  2:50019/1‑138 (MQ=17)                                                              tCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTTGAACCCCCCATACGGTGCCGTGTACACACTTTCCCGGGGGGGCCCTTCAGCCAccccggca                                                                                         >  1:865282/1‑143 (MQ=17)                                                                   gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCCAAACCCCCCAAAGTTGGCGTATAAAAACTTTTCAGGCGTGCTCCCTCTGCCACCCCGGCA‑cccccc                                                                                  >  1:1001921/1‑147 (MQ=17)                                                                     cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCa                                                                                <  2:812845/151‑1 (MQ=255)                                                                        acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAAt                                                                             <  1:762357/151‑1 (MQ=255)                                                                             cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGGTTGAACCCCCCCTACGTTTCCGTATAAACACCTTCCCGGCGTGCTCCTTCACCCACCCCGGCA‑ccccccccaaattttt                                                                        >  1:923277/1‑137 (MQ=14)                                                                              agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCTAACCCCCCATACGTTGCCCTATACACACTTTTCCGGCGGGCTCCTTCAGCCCCTCGGGCA‑CCCCCCCCAAAtttttt                                                                       >  1:1003983/1‑151 (MQ=14)                                                                               gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTg                                                                      <  2:840767/151‑1 (MQ=255)                                                                                  aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCTGACCCCCCAAACGTTGCCGTATACACACCTTTCCGGCGGGCTCCCTCCGCCACCCGGACA‑CCCCCCCCAAAtttttt                                                                       >  2:1103535/1‑147 (MQ=14)                                                                                  aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTa                                                                   <  2:1132301/151‑1 (MQ=255)                                                                                    aaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGT‑GGGGGGGGTTTCAACCCCCGGTAAGGTTCCGGATACACAATTTTCAAGCGTGCCCCTTCCGCCACCCGGCGA‑CCTCACCCAAATGTTTTTCTa                                                                   >  2:23616/1‑149 (MQ=9)                                                                                         gAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCATAAGTTGCCGTTTACACACTTTCCAGGCGGGCTCCCTCAGCCACCCGGACA‑CCCCACCAAATTTTTTTTCTACCaaa                                                              >  2:45474/1‑149 (MQ=14)                                                                                          aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCt                                                           <  1:437161/151‑1 (MQ=255)                                                                                              aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGGACCCCCCCTACGTTTTCGGATACACACCTTCCACGCGGGGTCCCTCCGCCACCCGGGCC‑CCCCACCCAAATTTTTTTCTACCaaaaccc                                                          >  2:922552/1‑144 (MQ=9)                                                                                                   cTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGGAACCCCCATACGTTTGCGTATACACACTTTTCCGGCGGGCTCCTTCAGCCACCCGGGCA‑CCCCCCCAAAATGTTTTTTTTCCAAACCCCCTgggggg                                                  >  2:453078/1‑151 (MQ=2)                                                                                                     cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGc                                                 <  2:111153/150‑1 (MQ=255)                                                                                                       cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCCCAACCTTGCCGTATACAAACTTTTCAGGCGTGCTCCCTCTGCCACCCGGACA‑CCCCACCCAAATTTTTTTTTCCCAAACCTCATGGGTGGGaa                                               >  2:22946/1‑150 (MQ=9)                                                                                                       cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTgg                                                  <  2:355345/147‑1 (MQ=255)                                                                                                        aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGGTCGAACCCCCCATACGTTTCCGTATACACACATTTCCGGCGGGGCCCCTCCGCCCCCCCGGCACCCCCCCCAAATTTTTTTCTTACAAAAcccctcggggg                                                   >  1:1016259/1‑137 (MQ=0)                                                                                                        aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAAACCCCCCCAAGGTGCCGCATACAAACTTTTCACGGGGGCTCCTTCAGCCCCCCCGGCA‑CCCCACCAAAATGTTTTTTTACCAAACCTCACGGGGGGCaa                                               >  2:72792/1‑149 (MQ=9)                                                                                                        aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACg                                             <  2:203286/151‑1 (MQ=255)                                                                                                           ttCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACggg                                           <  2:40450/150‑1 (MQ=255)                                                                                                            tCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACgggg                                          <  1:832470/150‑1 (MQ=255)                                                                                                               tACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCCATACGTTGCCGTGTACACACCTTTCAGGCGGGGTCCCTCCGCCCCCCGGGCA‑CCCCACCAAAATTTTTTTCTACCAAACCTCATGGGTGGGAACgggggg                                        >  2:437007/1‑147 (MQ=12)                                                                                                                   aGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGctacta                                  <  1:87017/151‑1 (MQ=255)                                                                                                                         tcttTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCGCTACGTTGTCCTATATACACCCTCCAGGCGGGGTCCCTTCGCCAATCCGGCA‑CCTCACCACAATGTTTTGTTACCACAACCCATTGGGGGGAACAGGGGGCAACTCAAggg                             >  1:815295/1‑150 (MQ=2)                                                                                                                         tcttTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAg                               <  1:216350/148‑1 (MQ=255)                                                                                                                           tttAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCGATACGTTGGCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACCCGGGCA‑CCCCACCCAAATGTTTTTTTACCAAACCTCATGGGGGGGAACGGGGGGCTACTATAGGGAg                           >  2:3530/1‑150 (MQ=17)                                                                                                                            ttAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAg                           <  2:851803/149‑1 (MQ=255)                                                                                                                                  atGGCGGTGA‑GGGGGGGGTTCGAACCCCCGATACGGTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACCCGGACA‑CCTCACCAAATTTTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAg                     >  1:502015/1‑149 (MQ=255)                                                                                                                                  atGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                   <  2:651166/151‑1 (MQ=255)                                                                                                                                      cGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               <  2:665895/151‑1 (MQ=255)                                                                                                                                      cGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               <  1:400647/151‑1 (MQ=255)                                                                                                                                      cGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               >  1:195420/1‑151 (MQ=255)                                                                                                                                        gTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               <  2:724759/149‑1 (MQ=255)                                                                                                                                             gggggggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaa         <  2:957187/151‑1 (MQ=255)                                                                                                                                                   gATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAACCGGTCAAgaagaa   <  2:1063182/151‑1 (MQ=255)                                                                                                                                                   gATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAgaagaa   <  2:875194/151‑1 (MQ=255)                                                                                                                                                    aTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  2:978979/151‑1 (MQ=255)                                                                                                                                                    aTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  >  1:633453/1‑151 (MQ=255)                                                                                                                                                    aTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGCGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  >  1:939321/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                        TTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925754‑926050

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