Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	144,184	0	T	G	42.1%	92.0 / 15.8	57	G161G (GGT→GGG)	yadH	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (12/21);  new base G (24/0);  total (36/21)
Fisher's exact test for biased strand distribution p-value = 1.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGCGTGGCGCGTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGTT  >  NC_000913/144035‑144325                                                                                                                                                      |                                                                                                                                                ggTGGCGTGGCGCGTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGgtg                                                                                                                                                <  1:430894/151‑1 (MQ=255)      cgtggcgCGTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTTCCTTGCGGGGTTGTTGAAAGGGGTGttt                                                                                                                                           >  1:588078/1‑151 (MQ=255)            gcgTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTTTTCTCCCTTTCGGGGTTGCTGAACGGGGTGTTTGCCaa                                                                                                                                      >  1:109483/1‑150 (MQ=255)             cgTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCTTTCGGGTTTTCTGGACGGTGGGTTTggccaa                                                                                                                                     >  2:93508/1‑145 (MQ=255)             cgTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTTCCTTTCGGGTTTTCTGAACGGGGTGTTTGCCaaa                                                                                                                                     >  2:1021603/1‑150 (MQ=255)              gTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAAAGGGGTGTTTTCCcaaa                                                                                                                                    >  1:857033/1‑150 (MQ=255)              gTGGTCTGTTTGTTGGCATTCTGGCGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTTCCTTTTGGGGTTGGTGCACGGGGTGTTTGTCaaaa                                                                                                                                    >  2:410477/1‑150 (MQ=255)                  tCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTTTTCTCCCCTGCGGGTTTGGTGAAAGGGGTGTTTTGCAAAACggt                                                                                                                                >  1:817687/1‑148 (MQ=255)                     gtttgttGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTTCCTTGCGGGGTTTCTGTACGGTGTGTTTGCCAAAACGTTC‑GAt                                                                                                                           >  1:589348/1‑151 (MQ=255)                       ttgttgGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATg                                                                                                                          <  1:315085/150‑1 (MQ=255)                         gttgGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGAc                                                                                                                        <  1:147144/150‑1 (MQ=255)                           tgGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCTTTCGGGTTTTCTGGACGGGGTGTTTTCCCAAACGTTt‑gatgaact                                                                                                                      >  2:493505/1‑147 (MQ=255)                           tgGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACAt                                                                                                                      <  1:214602/150‑1 (MQ=255)                                       gACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCa                                                                                                         <  1:768167/151‑1 (MQ=255)                                        aCGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGGTGTTTTTCCTTGCGGGTTTGCTGAACGGGGGGTTTTCCAAAACGTTC‑GATGACAACAGCCTGGTGcc                                                                                                          >  2:610805/1‑149 (MQ=255)                                              aTTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCTTTCGGGGTTGCTGAACGGGGTGTTTGCCAAAAAGTTT‑GATGACATCAGCCTGGGGGCAA‑CCttt                                                                                                  >  2:827605/1‑150 (MQ=255)                                                 tCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTtgtg                                                                                               <  2:827213/150‑1 (MQ=255)                                                            tgtgCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCt                                                                                                                <  2:160970/123‑1 (MQ=255)                                                            tgtgCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGGTTGCTGAACGGGGTGTTTTCCAAAACGTTC‑GATGACATCAGcc                                                                                                                 >  1:160975/1‑122 (MQ=255)                                                               gCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGGTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCCA‑CCTTTGTGTTTACGCCACTCAc                                                                                 >  1:886127/1‑150 (MQ=255)                                                                ccATTTCTGGTGCTTTCGTGGGTATTCGTTGCCTTAGCGCTGGTACTCACGGGGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCGCTCAAGt                                                                               <  2:621044/151‑1 (MQ=255)                                                                ccATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGt                                                                               <  2:126505/151‑1 (MQ=255)                                                                ccATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCCTGCGGGGTTTCTGAACGGGGTGGTTTCCACAACGTTC‑GATTACATCAGCCCGGGGCCCA‑ACTTTTTGTTTACGCCACTccc                                                                                 >  2:845492/1‑147 (MQ=255)                                                                 cATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTa                                                                              <  2:989371/151‑1 (MQ=255)                                                                  aTTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGGTTGCTGAACGGGGTGTTTTCCCAAAAGTTC‑GATGACATCAGCCTGGTGCCCA‑CCTTTGTGTTTACGCCCCTCCCGTAt                                                                             >  2:677004/1‑151 (MQ=255)                                                                     tCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCTTGGGGGTTTTCTGAAACGGGTGTTTTCCCAAAAGGTT‑GCTGGCATTTGCCCTGGGCCAA‑CCTTTGTGTTTAAGCCCCTCActgtttt                                                                           >  2:267303/1‑144 (MQ=255)                                                                     tCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGGGTTTTTTCTCCTTGCGGGTTTGATGAACAGGGGGTTTTGCCAAACGTTC‑GGTGAAATCAGCCTGGTGCCAA‑ACCTTTTGGTAAAGCCACTCACGTCttt                                                                           >  1:70874/1‑150 (MQ=255)                                                                      cAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTg                                                                          <  1:840837/150‑1 (MQ=255)                                                                      cAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTg                                                                          <  2:312478/150‑1 (MQ=255)                                                                       aGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCTTTCGGGGTTGCTGAAAGGGGGGTTTTGCAAAAAGTTC‑GATGACATCAGCCCGGGGCCAA‑CCTTTGTGTTAACGCCCCTCACGTATTTggg                                                                        >  1:301281/1‑151 (MQ=255)                                                                          tGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTTCGGGGTTGCTGAACGGGGTGTTTGCCCAAACGGTT‑GATGACATCAGCCTGGGGCCAA‑CCTTTGTTTTTACGCCACTCACGTATTTGGGCgg                                                                     >  2:655908/1‑151 (MQ=255)                                                                            cATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCgggg                                                                   <  2:577059/151‑1 (MQ=255)                                                                             aTTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCTTTCGGGGTTGGTGAACGGGGGGGTTGCCAAAACGTTC‑TGTGACATCAGCCTGGTGCCCA‑CCTTTGTGTTTACGCCCCCCACCTATTTgggggggg                                                                   >  2:472216/1‑150 (MQ=255)                                                                              ttCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGGTTGCTGAACGGGGTGTTTTCCAAAACGTTC‑GATGACATTAGCCTGGTGCCAA‑CCTTTTTGTTAACGCCACTCACGTATTTGGGCgggg                                                                   >  1:995469/1‑149 (MQ=255)                                                                               tCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGGTTGCTGAACGGGGGGTTTTCCAAAACGTTC‑GATGACATCAGCCTGGTGGCAA‑CCTTTTTGTTTACGCCACTCACGTATTTTGGCGGGGGCt                                                                >  1:701527/1‑151 (MQ=255)                                                                                cGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTTGTGTTTTCCCCTGCGGGGTTTCTTAAAAGGGGGTTTGCCAAAACGGTCTGTTGAAAT‑AGGCGTGTGCCAACCCTTTGTTTTAA‑GCCCCCCACATATTTTGGCGGGGGCtt                                                               >  1:924911/1‑151 (MQ=255)                                                                                       aTTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCa                                                        <  2:128042/151‑1 (MQ=255)                                                                                        ttCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGGTGAACGGGGGGGTTGCCCAAAAGTTC‑GATTGAATCAGCCTGGGGCCCA‑CCTTTTTGTTTACGCCACTCCCGTATTTTGGCGGGGTCTTTTtctc                                                         >  2:433117/1‑149 (MQ=255)                                                                                          cGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTTCCTTGCGGGGTTTCTGAACGGGGTGTTTGCCAAAAAGGTC‑GGTGACATCCGCCTGGTGCCAA‑CCTTTGTGTTTACGCCACTCACGTATTTTGGCGGGGGCTTTTTCt                                                          >  2:778104/1‑146 (MQ=255)                                                                                              gCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTTCGGGGTTTCTGAACGGGGTGTTTTCCAAAACGTTC‑GATGACATCAGCCTGGGGCCCA‑CCTTTGTGTTTACGCCACTCACGTATTTTGGGGGGGTCTTTTTCTCACTGACtt                                                 >  1:1038779/1‑151 (MQ=255)                                                                                               ccTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCCTGCGGGGTTTCTGAACGGGGTGTTTGCCAAAACGTTC‑GATGACATCAGCCCGGGGCCAA‑CCTTTGTGTTTACGCCACTCACGTATTTGGGCGGGGGCTTTTACTCACTGACttt                                                >  2:12536/1‑151 (MQ=255)                                                                                                  tAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTttgt                                              <  2:500920/150‑1 (MQ=255)                                                                                                  tAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTttgt                                              <  1:499283/150‑1 (MQ=255)                                                                                                   aaCGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGCCTttgttg                                            <  2:730761/151‑1 (MQ=255)                                                                                                     cGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGGTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGGCTTTTACTCACTGACTTTGTTgcc                                          >  1:46462/1‑151 (MQ=255)                                                                                                         ggTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTgccgc                                        <  1:441055/149‑1 (MQ=255)                                                                                                                  ggcggTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCa                              >  1:336199/1‑150 (MQ=255)                                                                                                                  ggcggTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCa                              >  2:877406/1‑150 (MQ=255)                                                                                                                   gcggTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCa                              >  1:328620/1‑149 (MQ=255)                                                                                                                     ggTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAggg                          >  2:862468/1‑151 (MQ=255)                                                                                                                     ggTGTCGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAggg                          <  2:997667/151‑1 (MQ=255)                                                                                                                        gttgttCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTg                       >  1:897743/1‑151 (MQ=255)                                                                                                                           gttCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCg                    >  2:728587/1‑151 (MQ=255)                                                                                                                           gttCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCg                    >  2:291638/1‑151 (MQ=255)                                                                                                                                   ttGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGa             >  1:472176/1‑150 (MQ=255)                                                                                                                                    tGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGaa            <  1:783396/150‑1 (MQ=255)                                                                                                                                             gCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGtt  <  2:322150/151‑1 (MQ=255)                                                                                                                                             gCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGt   <  1:827636/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                GGTGGCGTGGCGCGTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTC‑GATGACATCAGCCTGGTGCCAA‑CCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGTT  >  NC_000913/144035‑144325

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