Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,487,157	0	C	A	42.0%	103.1 / 32.5	69	intergenic (+194/‑78)	hrpA/ydcF	putative ATP‑dependent helicase/DUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base C (12/28);  new base A (29/0);  total (41/28)
Fisher's exact test for biased strand distribution p-value = 3.47e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGG  >  NC_000913/1487008‑1487307                                                                                                                                                        |                                                                                                                                                            cgcCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACg                                                                                                                                                            <  2:1175887/151‑1 (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:185871/150‑1 (MQ=255)              attaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCttt                                                                                                                                               <  1:429980/150‑1 (MQ=255)              attaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCttt                                                                                                                                               <  1:165857/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                                                    <  2:1024447/151‑1 (MQ=255)                                                                                                    cATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGAc                                                    <  1:439169/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)                                                                                                                                                        |                                                                                                                                                            CGCCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTAAGGACAT‑ACG‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTA‑AGGGAACGATAAGGAGG‑AACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGG  >  NC_000913/1487008‑1487307

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