Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,487,163	0	C	T	50.0%	58.2 / 39.5	68	intergenic (+200/‑72)	hrpA/ydcF	putative ATP‑dependent helicase/DUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base C (4/29);  new base T (33/0);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.83e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

ATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGC  >  NC_000913/1487021‑1487310                                                                                                                                                  |                                                                                                                                                        attaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCttt                                                                                                                                                  <  1:165857/150‑1 (MQ=255) attaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCttt                                                                                                                                                  <  1:185871/150‑1 (MQ=255) attaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCttt                                                                                                                                                  <  1:429980/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:1008783/151‑1 (MQ=255)                                                                                                              gCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATg                                <  1:982179/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)                                                                                                                                      cAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtgg         <  1:1004942/151‑1 (MQ=255)                                                                                                                                           cG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggctgg     <  1:54151/151‑1 (MQ=255)                                                                                                                                              aTTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGgcg   <  1:48351/151‑1 (MQ=255)                                                                                                                                                 tcctccTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGgcgc  >  1:1187546/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                        ATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGC  >  NC_000913/1487021‑1487310

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