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.9%	66.5 / 29.1	56	intergenic (+853/+170)	yhaC/garK	pentapetide repeats‑related protein/glycerate kinase I
Reads supporting (aligned to +/- strand):  ref base G (8/24);  new base T (24/0);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 4.75e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATA  >  NC_000913/3270312‑3270605                                                                                                                                                |                                                                                                                                                             aCGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGctctc‑t                                                                                                                                                       <  2:540643/150‑1 (MQ=255)  cGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑Tg                                                                                                                                                      <  2:623488/150‑1 (MQ=255)       cGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCCC‑TTTTGccc                                                                                                                                                >  1:600718/1‑149 (MQ=255)        gTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACt                                                                                                                                               <  1:989454/151‑1 (MQ=255)           ccTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TTTTGCACTgg                                                                                                                                             >  2:877274/1‑150 (MQ=255)           ccTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGCTTTcttcctatttccctcg                                                                                                                                              >  1:282517/1‑136 (MQ=255)               aTGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TGTTTCACTgg                                                                                                                                             >  1:688283/1‑146 (MQ=255)                tGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTTGCACTgcgt                                                                                                                                           >  1:857903/1‑145 (MQ=255)                 gAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑ggg                                                                                                                                   <  2:505242/151‑1 (MQ=255)                    ccccTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGttt                                                                                                                                <  1:385177/151‑1 (MQ=255)                     cccTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TTTTGCACATGT‑CG‑T‑GGGgtt                                                                                                                                >  2:312944/1‑150 (MQ=255)                       cTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TGGTGCCCCGGG‑CG‑T‑GGGgtttc                                                                                                                              >  1:30618/1‑148 (MQ=255)                             ggCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCa                                                                                                                        <  2:718348/150‑1 (MQ=255)                              gCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAgg                                                                                                                      <  2:278562/151‑1 (MQ=255)                                  tGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGt                                                                                                                   <  1:331404/150‑1 (MQ=255)                                      ccGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TGTTGGACCGGT‑CG‑T‑GGGGTTTCCCCCCCGGCGCT‑Acc                                                                                                              >  1:877158/1‑150 (MQ=255)                                      ccGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TTTTGCACTGGG‑CG‑T‑GGGGTTTCCCCCCCGGCGGT‑Accc                                                                                                             >  1:421814/1‑150 (MQ=255)                                        gggTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TGTTTCCCTCGT‑CG‑T‑GGGTTTTCCCCCCCGGCgcttcc                                                                                                               >  1:803411/1‑143 (MQ=255)                                              aGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCt                                                                                                     <  2:29793/151‑1 (MQ=255)                                                   aCCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTTTCTC‑TTTTTTACTGGT‑CG‑T‑GGGGTTCCCCCCCCGGCGTT‑TCCtcggcccccccccc                                                                                                >  1:749248/1‑138 (MQ=255)                                                       tGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑G‑GGGGTTCCCCCCCCCGCGGT‑ACCCTGCACCCCGccccctg                                                                                             >  1:865362/1‑146 (MQ=255)                                                        gCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATgg                                                                                            <  2:1009449/150‑1 (MQ=255)                                                         ccACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAg                                                                                          <  1:523502/151‑1 (MQ=255)                                                         ccACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAg                                                                                          <  1:335575/151‑1 (MQ=255)                                                          cACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TTTTTCACCGGT‑CC‑T‑GGGGTTCCCCCCCCCGCGTT‑AACCTGCCCCCCGCCCCCTGGAg                                                                                          >  2:25658/1‑150 (MQ=255)                                                          cACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGGTTTCTTTC‑TTTTGCACCGGG‑CC‑T‑GGGGTTCCCCCCCCAGCGGT‑ACCTGGGCCCCCGccccttcgcg                                                                                          >  1:289403/1‑143 (MQ=255)                                                           aCGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TTTTTCCCTTGT‑CG‑T‑GGGGTTTCCCCCCCCGCGGT‑ACCTCGCAccccccccctttggggc                                                                                        >  2:193162/1‑137 (MQ=255)                                                               aCTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTTCACTGGT‑CG‑T‑GGGGTTTCCCCCCCGGCGTT‑ACCTGGCACCCCGCCCTATGGAG‑CCCgg                                                                                    >  2:687438/1‑150 (MQ=255)                                                                 tGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TGTTGGACTTGG‑CG‑T‑GGGGTTTCCCCCCCGGCGGT‑TCCCCGCACCCCCCCCCTTTGAG‑ccccggcc                                                                                 >  2:476788/1‑146 (MQ=255)                                                                  gTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACt                                                                                 <  1:476768/150‑1 (MQ=255)                                                                       cAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTTTC‑TGTTTCACCGGT‑CG‑T‑GGGGTTCCCCCCCCCGCGGT‑TCCTTGGACCCCTCCCTCTGGAG‑TCCgggcct                                                                                >  2:798713/1‑142 (MQ=255)                                                                          ccGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTcctcccc                                                                        <  1:345656/151‑1 (MQ=255)                                                                             cgcgGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTcctcccctcc                                                                     <  2:693046/151‑1 (MQ=255)                                                                             cgcgGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTcctcccctcc                                                                     <  2:845476/151‑1 (MQ=255)                                                                              gcgGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCg                                                                    <  1:325184/151‑1 (MQ=255)                                                                                          tACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TGTTTCACTCGG‑CG‑T‑GGGGTTTCCCCCCCCGCGTT‑TCCTGGCCCCCCCCCCCATTGAG‑GCCGGGCTTTTCCCCCCCCCCCCCCTCcccccc                                                        >  2:361329/1‑151 (MQ=255)                                                                                          tACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCGCCCGTCTcccc                                                         <  2:707367/150‑1 (MQ=255)                                                                                                       cACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TTTTTCACTGGT‑CG‑T‑GGGGTTTCCCCCCCGGCGTT‑TCCCGGCACCCCGCCCCCTGGAG‑GCCGGACTTTCCTCCCCCCCCCCCCTCCCCCCCGCAGAGgac                                               >  1:854997/1‑147 (MQ=255)                                                                                                        aCCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTGTCCC‑CGGTTCACTGGTCCG‑G‑GGGTTTCCCCCCC‑CGCGTT‑TCCTGGCACCCCGCCCCAATGAG‑CCCCGAACTTTCTCCCCTCCCCCCCGCTCCCCCGGAGAAggcga                                             >  2:192709/1‑148 (MQ=255)                                                                                                        aCCCTTACCTGATCCCGCTT‑‑GCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGgacgacgac                                         <  2:263196/150‑1 (MQ=255)                                                                                                         cccTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TTTTTCACTGGT‑CG‑T‑GGGTTTTCCCCCCCCGCGTTAACCTGGC‑CCCCCCCCTATGGGGCCCCGGCCTTTCC‑CCCCCCCCCCCCCCTCCCCCCCAGAgggc                                               >  2:249103/1‑143 (MQ=255)                                                                                                             tACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTTTCTC‑TTTTTCACTCGG‑CG‑T‑GGGGTTTCCCCCCCCGCGTT‑ACCTTGCACCCTCCCCTCTTGGG‑CCCGGACCTTTCTCCCCTCCCCCCCCCTCCCCCCAAGAGgacgac                                            >  1:894042/1‑144 (MQ=255)                                                                                                              aCCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAg                                     <  1:207872/150‑1 (MQ=255)                                                                                                               ccTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTTTCTC‑TTTTGCCCTTGG‑CG‑C‑GGGGTTTCCCCCCACGCGTT‑TACCGGCCCCCCGGCCCATGGAG‑CCCGGACTTTCCTTCCCTCCCCCCCTCTCCCCCCCAGAGGGCGACGACCaa                                      >  1:761525/1‑148 (MQ=255)                                                                                                                cTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TGTTGGACTTGT‑CG‑TGGGGTTTCCCCCCCA‑GCGGT‑ACCTGGCCCCCCGCCCCATGGGG‑CCCGGGCCTTCCTTCCCCCCCCCCGCCTCCCCCCAAGAGGGCCACCACGAAacgg                                  >  2:350194/1‑151 (MQ=255)                                                                                                                  gATCCCGCTTGCGCGGGCCATCGGCGGTTTTGTCTC‑TGTTGCACCGGG‑CG‑T‑GGGTTTTCCCCCCCGGCGTT‑ACCTGGCACCCCGCCCCATGGGG‑CCCGGGCTTTCCTCCCCCCCCCCCCTCTCCCCCGCAGAGGGCGACGACGaa                                      >  2:363368/1‑145 (MQ=255)                                                                                                                  gATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAgcgg                                  <  1:117925/149‑1 (MQ=255)                                                                                                                  gATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCAACGGG‑CG‑G‑GGGTTTTCCCCCCCGGCGGT‑TCCTGGCACCCCCCCCCATTGGG‑CCCCGGCTTTTCTCCCCCCCCCCCGTCCCCCCCCGACAGgacgaccaagcaacgg                                  >  1:937412/1‑139 (MQ=255)                                                                                                                        gCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TTTTTCACTGGG‑CG‑T‑GGGGTTTCCCCCCCGGCGTT‑ACCTGGCACCCCGCCCCATGGAG‑CCCCGGCTTTTCTCCCCCCCCCCCCTCTCCCCCCCAGAGGACGACGACGAAGGGGCGCCTGt                           >  2:277518/1‑150 (MQ=255)                                                                                                                         cTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TGTTGCCCCTGT‑CGTT‑GGGTTTCCCCCCC‑CGCGTT‑TCCTGGCACCCTGCCCCCTGGAG‑CCCGGGCTTTTCTCCCCCCCCCCCCTCCCCCCCGAAGAGGACGACGACCAAGCGGCGACCGt                           >  1:942615/1‑149 (MQ=255)                                                                                                                             cgcgGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTc                     <  1:901238/151‑1 (MQ=255)                                                                                                                              gcgGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCa                    <  2:484072/151‑1 (MQ=255)                                                                                                                                gggCCATCGGCGGTTTGCTCTC‑TTTTGCCCTGGT‑CG‑T‑GGGTTTTCCCCCCCGGCGTT‑ACCTGGCACCCCGCCCCATGGAG‑CCCGGGCTTTCCTCCCCCCCCCCCCTCTCCCCCCCAGAGGGCGGCGACGAAGCGGCGGCTGTCTTGTccg                   >  1:323125/1‑148 (MQ=255)                                                                                                                                    cATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGGTTCCCCCCCCGGCGGT‑ACCTGGCACCCTGCCCCATGGAG‑CCCGGGCTTTCCCCCCCCCCCCCCCCCCCCCCCCCAGAGGGCGACGACGAAGCGGCGACTGTCTGGTCAGCt                 >  1:271614/1‑148 (MQ=255)                                                                                                                                       cggcggTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGgc            <  1:446498/150‑1 (MQ=255)                                                                                                                                               tGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTT‑CCCCCCAGGCGTT‑AACTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCGCCCGTCCCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGGCAGCTTCGGCGCGCAg      >  2:562253/1‑147 (MQ=255)                                                                                                                                                gCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGtata  <  1:406384/151‑1 (MQ=255)                                                                                                                                                gCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGtat   <  1:172185/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                             ACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGT‑CG‑T‑GGGTTTCCCCCCCAGGCGTT‑ACCTGGCACCCTGCCCTATGGAG‑CCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATA  >  NC_000913/3270312‑3270605

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