Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,270,455	0	G	T	42.6%	65.1 / 25.2	54	intergenic (+853/+170)	yhaC/garK	pentapetide repeats‑related protein/glycerate kinase I
Reads supporting (aligned to +/- strand):  ref base G (11/20);  new base T (23/0);  total (34/20)
Fisher's exact test for biased strand distribution p-value = 3.13e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATA  >  NC_000913/3270308‑3270605                                                                                                                                                    |                                                                                                                                                            cAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGctc                                                                                                                                                          <  1:207561/151‑1 (MQ=255)    tACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTgtctc                                                                                                                                                        >  1:400525/1‑150 (MQ=255)     aCGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTg                                                                                                                                                      <  1:388497/151‑1 (MQ=255)         gCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGGTTTCTCTCtttt                                                                                                                                                    >  1:114264/1‑149 (MQ=255)          ccGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCtttt                                                                                                                                                    >  1:432778/1‑148 (MQ=255)              aCCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑gg                                                                                                                                            <  2:495438/151‑1 (MQ=255)                 ttATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGTCCT‑gg                                                                                                                                            >  1:114611/1‑148 (MQ=255)                     gAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑gg                                                                                                                                   <  1:224759/150‑1 (MQ=255)                     gAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑gg                                                                                                                                   <  1:673980/150‑1 (MQ=255)                            tATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTGTCTCTTTTGGACT‑GGT‑C‑GT‑GGG‑GTTcccc                                                                                                                          >  1:896930/1‑151 (MQ=255)                               ttGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTcccccc                                                                                                                        <  2:119672/150‑1 (MQ=255)                                tGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTccccccc                                                                                                                       <  2:954653/150‑1 (MQ=255)                                 ggCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCa                                                                                                                      <  2:198049/150‑1 (MQ=255)                                      tGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTTTTGCACT‑GGG‑C‑GG‑GGG‑TTTTCCCCCCCGGCggt                                                                                                                >  1:755428/1‑149 (MQ=255)                                        cTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTTTTGCACT‑GGG‑CGGG‑GGG‑TTTccccccccgcgctac                                                                                                               >  1:38591/1‑143 (MQ=255)                                         tCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACT‑GGT‑C‑GT‑GGG‑GTTTCCCCCCCCGCGTTTcc                                                                                                             >  1:231141/1‑151 (MQ=255)                                          ccGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGCGTTGTCTCTGTT‑CACTAGGC‑C‑AG‑GGG‑GTTTCCCCCCATGCGCTAcc                                                                                                             >  1:59093/1‑150 (MQ=255)                                                ggAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑GTTCCCCCCCCGGCGGTACCTCGCAcc                                                                                                      >  2:598916/1‑151 (MQ=255)                                                      tACCGTGCCACGGACTGTTACTAGCCGCGCGGTGCGCTCTTACCGCACCC‑TTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGccc                                                                                                <  2:426236/150‑1 (MQ=255)                                                        ccGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCCGGCGTTACCTGGCACCCTGccc                                                                                                >  2:445497/1‑149 (MQ=255)                                                          gTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACT‑GGT‑C‑GT‑GGG‑GTTTCCCCCCCCGCGGTTCCTGGCAcccccccc                                                                                                >  1:634281/1‑142 (MQ=255)                                                           tGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACT‑GGT‑C‑GT‑GGG‑GTTCCCCCCCCGGCGGTTCCCGGCAACCCGCCCTAttg                                                                                           >  2:479737/1‑149 (MQ=255)                                                             ccACGGACTGTTATCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAg                                                                                         <  2:794705/151‑1 (MQ=255)                                                                    cTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGGTTGCTCTCTGTTGCACT‑GGT‑C‑GTGGGG‑TTTCCCCCCCA‑GCGGTACCTGGCCCCCCGCCCCATGGAGccccg                                                                                    >  2:679162/1‑147 (MQ=255)                                                                       ttACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACC‑GGG‑C‑GT‑GGG‑GTTTCCCCCCAGGCGTTACCTGGCACCCCGCCCTATTGGGCCCCGGCttt                                                                               >  2:106427/1‑151 (MQ=255)                                                                           cAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTTTTTCTTTTGCACT‑GGT‑C‑GT‑GGG‑GTTCCCCCCCCGGCGTTACCCCGCACCCTGCCCCATGGAGCCCGGACCTtcctc                                                                           >  2:91333/1‑151 (MQ=255)                                                                           cAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTtcctc                                                                           <  1:732382/151‑1 (MQ=255)                                                                                gcgcgGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTTTGTTGCCCT‑GGT‑C‑GG‑GGG‑GTTCCCCCCCCCGCGGTTCCCGGCACCCCACCCCCTCGAGCCCGGGCTTTtcccccctc                                                                      >  2:951344/1‑151 (MQ=255)                                                                                      tGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGGTTTCTCTCTTTTTCACT‑GGT‑C‑GT‑GGG‑GTTCCCCCCCAGGCGGTACCTGGCACCCTGCCCTTTGGAGCCCGGACTTTTCTCcccccccccc                                                                 >  2:555969/1‑150 (MQ=255)                                                                                          ctctTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑GTTTCCCCCCCGGCGGTACCTGGCACCCCGCCCCATGGAGCCCCGGCTTTCCTCCCCCCCCCCCCtctc                                                            >  2:429277/1‑151 (MQ=255)                                                                                              tACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑GTTTCCCCCCCGGCGGTACCTGGCACCCTGCCCTATGGAGCCCCGACTTTTCTCCCCCCCCCCCGTCTccccc                                                        >  2:464962/1‑151 (MQ=255)                                                                                                   cACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCCCT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCCCGCGTTACCCGGCCCCCTGCCCCCTTGAGCCCCGGCTTTCCTcccctccccccccctccccc                                                        >  1:196073/1‑146 (MQ=255)                                                                                                    aCCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACT‑GGT‑C‑GT‑GGG‑GTTTCCCCCCCGGCGGTTCCTGGCCCCCCCCCCCATGGGGCCCCGGCTTTCCTCCCCCCCCCCCCGCTccccccca                                                     >  1:963590/1‑145 (MQ=255)                                                                                                    aCCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGATTTCTCCCTTTTTTCCT‑TGT‑C‑GT‑GGGCTTTCCCCCCCA‑GCTTTTCCTGTCCCCCCGTCCTATGGTGTCCGCACATTTCCCCCCTCCCCCCCGCTCCCCCgc                                                      >  2:655791/1‑146 (MQ=255)                                                                                                          tCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGgacga                                             <  1:279524/150‑1 (MQ=255)                                                                                                             cccTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGgacgacga                                          <  1:865821/150‑1 (MQ=255)                                                                                                              ccTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTCGTGCAAC‑GGT‑C‑GG‑GGG‑GTTCCCCCCCAGGCGGTACCCGGCCCCCCGCCCTATGGAGCCCCGACTTTTCCCCCCCCCCCCCCTCTCCCCCCCAAAGggcgacg                                           >  2:32881/1‑148 (MQ=255)                                                                                                                 tACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCCCT‑GGT‑C‑GT‑GGG‑GTTTCCCCCCCCGCGTTTCCTTGCACCCTGCCCTATGGAGGCCCGACTTTCCTCCCCCCCCCCCCTCTCCCCCCACGAGggcga                                             >  2:495836/1‑143 (MQ=255)                                                                                                                 tACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGaa                                      <  2:428742/150‑1 (MQ=255)                                                                                                                   ccTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTGTCTCTGTTTCACT‑GGT‑C‑GT‑GGG‑TTTTCCCCCCCCGCGTTACCCGGCACCCCGCCCTATGGAGGCCGGGCTTTTCTCCCCCCCCCCCCCCTCCCCCCCAGAGgacgacg                                           >  2:9922/1‑143 (MQ=255)                                                                                                                    cTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCCGTTTCACT‑TGTCC‑GG‑GGG‑TTTCCCCCCCCGCCG‑TAACTGGCACCCCGCCCCATCGAGGCCCGACCTTTCCCCCCCCCCCCCGTCTCCCCCCCAGAAGACGACGACCaa                                      >  2:376206/1‑147 (MQ=255)                                                                                                                      gATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAgcgg                                  <  1:464952/149‑1 (MQ=255)                                                                                                                       aTCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGGCCC‑TGT‑C‑GC‑GGG‑GTTTCCCCCCAAGCGGTACCCCGCCACCCGCCCCATCGAGACCCGGCCTTCCCCCCCCCCCCCCGCCTCCCCCCCAGAGGACGACGAAAAAacggcg                                >  2:831846/1‑150 (MQ=255)                                                                                                                            gCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACT‑GGT‑C‑GT‑GGG‑GTTCCCCCCCCGGGGTTACCCGGCACCCCGCCCCATGGGGCCCCGGCTTTTCTCCCCCCCCCCCCCCTCCCCCCCAGAGGGCGACGACGAAacgg                                  >  1:447609/1‑143 (MQ=255)                                                                                                                              ttGCGCGGGCCATCGGCGGTTTTCTCTCTTTTTCACT‑GGT‑C‑GT‑GGG‑GTTTCCCCCCCGGCGTTACCTGGCACCCCGCCCCATTGAGCCCCGGCTTTCCTCCCCCCCCCCCCTCCCCCCCCAAGAGGACGACGACGAAacgg                                  >  2:488505/1‑141 (MQ=255)                                                                                                                                  gcgGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCa                    <  2:843253/151‑1 (MQ=255)                                                                                                                                       ccATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTc               <  2:198068/151‑1 (MQ=255)                                                                                                                                        cATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCg              <  2:231145/151‑1 (MQ=255)                                                                                                                                           cggcggTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGgcg           <  1:445487/151‑1 (MQ=255)                                                                                                                                           cggcggTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGgc            <  2:114262/150‑1 (MQ=255)                                                                                                                                            ggcggTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGgcgc          >  2:500830/1‑151 (MQ=255)                                                                                                                                               ggTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGgcgcgc        >  1:907431/1‑150 (MQ=255)                                                                                                                                                gTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGgcgcgc        >  2:586180/1‑149 (MQ=255)                                                                                                                                                gTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGgcgcgc        >  1:779001/1‑149 (MQ=255)                                                                                                                                                    gCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGtata  <  2:896947/151‑1 (MQ=255)                                                                                                                                                    gCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGtata  <  1:951325/151‑1 (MQ=255)                                                                                                                                                    gCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGtat   <  2:672133/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                            CAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACT‑GGT‑C‑GT‑GGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATA  >  NC_000913/3270308‑3270605

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