Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,487,150	0	G	A	41.4%	94.0 / 29.6	72	intergenic (+187/‑85)	hrpA/ydcF	putative ATP‑dependent helicase/DUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base G (13/28);  new base A (29/0);  total (44/28)
Fisher's exact test for biased strand distribution p-value = 1.07e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTTAGCGCCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGAC  >  NC_000913/1487001‑1487298                                                                                                                                                       |                                                                                                                                                           cATTTAGCGCCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAgg                                                                                                                                                           <  2:244348/151‑1 (MQ=255)      aGCGCCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑a                                                                                                                                                     <  1:831654/151‑1 (MQ=255)       gcgcCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAAGACaaa                                                                                                                                                      >  2:823705/1‑148 (MQ=255)        cgcCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACg                                                                                                                                                   <  2:1175887/151‑1 (MQ=255)        cgcCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGaaaaa                                                                                                                                                      >  1:360557/1‑145 (MQ=255)            gggTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGCTAAAGaaaaacaaatt                                                                                                                                                >  1:734826/1‑141 (MQ=255)              gTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGAAAT‑AAG‑Aatt                                                                                                                                              >  2:1176010/1‑149 (MQ=255)                   ttattaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGAAAAA‑ACC‑Attt                                                                                                                                              >  2:936091/1‑144 (MQ=255)                     attaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCttt                                                                                                                                      <  1:429980/150‑1 (MQ=255)                     attaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCttt                                                                                                                                      <  1:165857/150‑1 (MQ=255)                     attaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCttt                                                                                                                                      <  1:185871/150‑1 (MQ=255)                      ttaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTc                                                                                                                                     <  1:304192/150‑1 (MQ=255)                       taATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGAAAA‑ACG‑ATTTTCTTCttt                                                                                                                                      >  1:540593/1‑148 (MQ=255)                         aTTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAAGACAA‑AAG‑ATTTTCTTCTTTCag                                                                                                                                   >  1:9011/1‑149 (MQ=255)                         aTTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCagag                                                                                                                                 <  1:178669/151‑1 (MQ=255)                         aTTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑AAA‑AATTTCTCCttt                                                                                                                                      >  2:192006/1‑146 (MQ=255)                         aTTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGAAAA‑AAG‑ATTTTCTTCttt                                                                                                                                      >  1:706369/1‑146 (MQ=255)                           tCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGAAAT‑AAA‑ATTTTCTCCTTTCTGAGt                                                                                                                                >  1:1079534/1‑150 (MQ=255)                               aaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTcc                                                                                                                           <  1:671105/151‑1 (MQ=255)                                         ggTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACt                                                                                                                 <  2:1078348/151‑1 (MQ=255)                                           tAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGGTTAAGAAAAA‑AAA‑ATTT‑CTTCTTTCTGAGTGGTcccc                                                                                                                         >  2:569832/1‑139 (MQ=255)                                                aTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGAAAA‑ACGAATTTTTCCCTTTCA‑AGGGCTCCCCCTCTCaacatttt                                                                                                            >  1:852128/1‑142 (MQ=255)                                                 ttCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGACAA‑AAG‑ATTTTCTCCTTTTTGAGGGGTCCGCCTTTCCCtatttt                                                                                                            >  1:537878/1‑146 (MQ=255)                                                  tCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGAAAA‑ACG‑AATTTCTCCTTTTCGAGGGCTCCGCTTCTCCCTATTTtct                                                                                                          >  1:315805/1‑149 (MQ=255)                                                        aTTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGGTTAAGACAA‑AAA‑AATTTCTCCTTTTCTAGGGGTCCCCTTCTCCCTCTTAtctctccc                                                                                                     >  2:902730/1‑144 (MQ=255)                                                         ttCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAGGAAAA‑AAC‑ATTTTCTCCTTTTTCAGTGGGCCGCTTCTCACcatttttt                                                                                                          >  1:725703/1‑138 (MQ=255)                                                             gAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑t                                                                                           <  2:369618/151‑1 (MQ=255)                                                               aaaaTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGACAT‑AAA‑ATTTTCTCCTTTCTTGGTGCTCCGCCTCTCCCTATTTTTTCACGCAGAA‑Ttt‑tt                                                                                          >  2:538730/1‑150 (MQ=255)                                                                  aTACAATTCTCTGCTGCAAGATGAATATTGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGACAA‑ACG‑AATTTCTCCTTTCATAATGCTCCCCCTCTTACTATTTTCTCACGCAGTA‑Ttt‑tt                                                                                          >  1:1192225/1‑147 (MQ=255)                                                                  aTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑Agg                                                                                     <  2:281540/151‑1 (MQ=255)                                                                     cAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGaa                                                                                  <  1:1128324/151‑1 (MQ=255)                                                                           tctGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAGGAAAA‑ACG‑ATTTTCTCCCTTTTGGGGGCTCCCCCTCTCCCTTTTTTTTCACCCAGGATTTT‑TTA‑AGGGaa                                                                                  >  1:861354/1‑146 (MQ=255)                                                                           tctGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAAGAAAA‑ACG‑AATTTCTCCCTTCCGAGGGGTCCCCTTCTCACTATTTTCTCCCCCAGTA‑TTT‑TTT‑AAGGAAAGaaa                                                                             >  1:238336/1‑148 (MQ=255)                                                                              gctgcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAagga                                                                         <  2:290125/151‑1 (MQ=255)                                                                               ctgcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAAGAGAA‑AAA‑AATTTCTCCTTTTTGCGGGGTCCCCCTCTCCCCTTTTTTTCCCGCAGGA‑TTCTTTA‑AGGGAA‑GATaaa                                                                           >  1:1159563/1‑147 (MQ=255)                                                                                tgcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAAGACAA‑AAA‑ATTTTCTCCTTTTCTGGGGGTCCCCCTCTCACTTTTATCTCACCACCCA‑TTT‑TTA‑ATGGGACGATaaagggg                                                                       >  1:620366/1‑146 (MQ=255)                                                                                tgcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGAAAT‑AAA‑ATTTTTTCCTTTTAGAGTGATCCGCCTTTCACTATTTTCTCTCACCAAA‑TTT‑TTT‑TTGGAACAATAaggagg                                                                       >  2:859879/1‑151 (MQ=255)                                                                                tgcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGAAAA‑ACG‑ATTTTTTTCTTTTAGAGTGGTCCGCTTCTCACTATTATTTTCCCCAGTA‑TTT‑TTT‑AGGGAAAGaaaaagggg                                                                       >  2:1013177/1‑146 (MQ=255)                                                                                  cAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAAGAAAT‑AAG‑AATTTCTCCTTTTATAGTGGTCCCCTTCTCCCTATTTTTTCTCGCAGTATTTT‑TTA‑AGGGAA‑GATAaagagg                                                                       >  1:846026/1‑149 (MQ=255)                                                                                    aGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAAGAAAA‑AAG‑AGTTTCTCCTTTTTGAGTGGTCCGCCTCTCCCTATTTTTTTACCCAGCATTTT‑TTA‑AGGGAAAGA‑AAAGAGG‑aaa                                                                   >  1:851507/1‑149 (MQ=255)                                                                                        gAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGACAA‑AAG‑ATTTTCTCCTTTTTGAGAGCTCCCCTTCTCACTATTATTTCACCCAGGA‑TTT‑TTA‑AAGGAAAGAAAAAGAGG‑aa                                                                    >  1:1120965/1‑145 (MQ=255)                                                                                         aataatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTACGACAA‑AAG‑CTTTTTTTTTTTTATAGAGCTCCGCCTCTCACTATTATCTTACGACGTA‑TTC‑TTT‑AGGGGACGATAAAGAGG‑AACCaa                                                                >  1:875930/1‑147 (MQ=255)                                                                                          ataatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAAGAAAA‑AAG‑AATTTCTCCTTTTTTAGGGGTTCGCTTCTCCCTATTATTTTACGCAGGA‑TTT‑TTA‑AAGGGAAGATAAAGAGG‑AACCATGaa                                                             >  1:163827/1‑150 (MQ=255)                                                                                             atGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACAtt                                                         <  2:635408/151‑1 (MQ=255)                                                                                               gTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTAc                                                       <  2:307953/151‑1 (MQ=255)                                                                                                 ttATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAGGACAT‑AAG‑ATTTTCTCCTTTTCGAGTGGTCCCCTTCTCACCATTTTTTCACCCCGTA‑TTT‑TTT‑AAGGAAAGAAAAAGAGG‑GAACAAGAAAATTTccc                                                     >  1:918671/1‑151 (MQ=255)                                                                                                    tCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGAAAA‑AAG‑ATTTTCTCCTTTTTGAGTGCTCCCCCTCTCCCCATTATTTCACGCAGTA‑TTT‑TTT‑AGGGGACGAAAAAGAGG‑AAACATGAACATTAcccc                                                    >  1:26476/1‑149 (MQ=255)                                                                                                      tACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGAAAT‑AAG‑ATTTTCTCCTTTCCGAGTGGTCCGCTTCTCCCTTTTATTTCACGCCGTA‑TTT‑TTT‑AAGGGACGAAAAAGAGG‑GAACATTAACATTACCCCCttt                                                >  1:136664/1‑151 (MQ=255)                                                                                                       aCAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTAccctg                                                   <  2:172149/147‑3 (MQ=255)                                                                                                       aCAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCg                                                   <  1:685838/147‑1 (MQ=255)                                                                                                          gCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGAAAA‑AAG‑ATTTTCTCCTTTTAGAGTGCTCCGCCTCTCCCTATTTTCTCACGCCGTA‑TTC‑TTT‑AGGGAACGATAAAGAGG‑AACCATGAACATTAcccccctt                                                >  1:426281/1‑143 (MQ=255)                                                                                                           cATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGAc                                           <  1:439169/151‑1 (MQ=255)                                                                                                           cATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGAc                                           <  2:1024447/151‑1 (MQ=255)                                                                                                              ttCCTTAAAAGATATGTCAGGCTTGCGGAGCGGC‑GGCTTCAGAAAA‑AAG‑ATTTCCTCCCTTTAGATTGTTCCTCTTTTCACTAATATCTCCCGCAGTT‑TTC‑TTATAGGGAA‑GTGAAGGAGGAAACAATGAACATTTCCCCGTTTccactcc                                         >  1:111026/1‑146 (MQ=255)                                                                                                                 cTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGAAAT‑AAA‑ATATTCTCCTTTCTTTGGGGTCCGCCTCTCACCATTTTTTTACGCAGGATTTT‑TTA‑AGGGAA‑GATAAAGAGG‑AAACATGAAAAATAACCCCTTTTCGACGCctt                                      >  2:884739/1‑150 (MQ=255)                                                                                                                  ttAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAAGAAAA‑AAG‑AATTTCTCCTTTTCGAGTGCTCCCCCTTTCACTATTTTTTCACCCCGTA‑TTCTTTA‑AGGGAA‑GATAAAGAGG‑AAACATGAACATTACCCCCCTTCCCACGCttt                                      >  1:203976/1‑149 (MQ=255)                                                                                                                   tAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGAAAT‑ACA‑AATTTCTTCTTTCTTAATGGTCCGCTTCTCACTATTTTTTCCCCCCGTA‑TTT‑TTT‑ATGGGAAGAAAAGGAGG‑AACCATGAACAATAACCCCCTTCCGCCGCttttg                                    >  2:664511/1‑148 (MQ=255)                                                                                                                    aaaaGATATGTCAGGCTTGCGGAGTGGC‑GGTTAATGAAAA‑AAA‑TATTACCCCCTTCAGAGTGCTCCGCCTTTCCCTAATATCTTACAACCTG‑TTT‑TTT‑AGGGGAAGAAAAAGGGG‑AACAACGAACATTATCCCGttttcccggccttaacc                                  >  1:904163/1‑137 (MQ=255)                                                                                                                     aaaGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAAGAAATAAAA‑ATTTTTTCCTTT‑TTAGTGGTCCGCCTCTCTCTATTTTTTCACGCCGTATTTT‑TTA‑AGGGAA‑GGTAAAGAGG‑GAACATGAAAATTAACCCGCTTTCCACGCTTTTGcc                                  >  2:209623/1‑150 (MQ=255)                                                                                                                     aaaGATATGTCAGGCTGGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAGCGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCg                                 <  2:130795/151‑1 (MQ=255)                                                                                                                       aGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAAGAAAA‑AAG‑ATTTTCTCCCTTTAGAGGGGTTCCCCTCTCACTATTTTTTCACGCCGGATTTT‑TTA‑AGGGAA‑GATAAAGAGG‑AAACATGAAAATTACCCCCGTTTCGACGCCTTTGcccg                                >  1:1088590/1‑148 (MQ=255)                                                                                                                       aGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAAGACAA‑AAG‑ATTTTCTCCTTTTTGAGGGCTCCGCCTCTCACTCTTTTTTCACCCCGGA‑TTT‑TTT‑AAGGAACGATAAAGAGG‑GAACATGAAAAATACCCCCCTTCCGGCGCCTTTGcccg                                >  1:1108112/1‑148 (MQ=255)                                                                                                                             gTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACtat                         <  1:1179409/151‑1 (MQ=255)                                                                                                                                  gCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATg                    <  1:982179/151‑1 (MQ=255)                                                                                                                                  gCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATg                    <  1:1008783/151‑1 (MQ=255)                                                                                                                                       cGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATa               <  2:464309/151‑1 (MQ=255)                                                                                                                                             ggC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGtt         <  1:883282/151‑1 (MQ=255)                                                                                                                                             ggC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGt          <  1:1035729/150‑1 (MQ=255)                                                                                                                                              gC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTa        <  2:915116/151‑1 (MQ=255)                                                                                                                                               c‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTAt       <  1:610269/151‑1 (MQ=255)                                                                                                                                                  gTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTATCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGACGCCATAAATGTTATCg     >  1:747951/1‑151 (MQ=255)                                                                                                                                                   ttAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCgg    >  1:573083/1‑151 (MQ=255)                                                                                                                                                     aaGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGAc  >  2:92487/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                           CATTTAGCGCCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGAC  >  NC_000913/1487001‑1487298

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