Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,350	0	C	T	54.4%	10.7 / 58.8	57	S64F (TCT→TTT)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (7/19);  new base T (31/0);  total (38/19)
Fisher's exact test for biased strand distribution p-value = 1.03e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-08

CGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATT  >  NC_000913/4604210‑4604500                                                                                                                                               |                                                                                                                                                           cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCTGTCAGGAGTCATTTCAGGATGCCATGTCGCGGACCTCGTTTGCGGCGGT‑CCTT‑TTTTTTTTTTTTTCCaccgga                                                                                                                                                  >  2:121164/1‑146 (MQ=255)     aTGAGTAGTATCGGTATTGAGAGCTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACAAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCATTGttttttttttctcctagcgagtg                                                                                                                                                >  1:451532/1‑140 (MQ=255)       gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCTT‑TTTTTTTTTTTTGCCCTGAGAAggggg                                                                                                                                            >  2:312880/1‑147 (MQ=255)       gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTTCCCTGAGAAggggg                                                                                                                                            >  2:93334/1‑147 (MQ=255)       gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTGCCATGAGAAGTGAg                                                                                                                                            >  1:32672/1‑151 (MQ=255)          tagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTCTGCCCTGAGAAGTGAgcg                                                                                                                                          >  1:302659/1‑150 (MQ=255)          tagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                                                                                           <  1:53122/149‑1 (MQ=255)            gtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTTCCATGGGAAGGGGGGGCag                                                                                                                                       >  2:384200/1‑151 (MQ=255)             taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTTCCATGAGGAGGTAGCGGagg                                                                                                                                      >  1:436486/1‑150 (MQ=255)                 ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCTT‑TTTTTTTTTTTTTCCCTGAGAAGGGAGCGCCGAGAgg                                                                                                                                  >  2:376537/1‑151 (MQ=255)                  gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTTCCCTgaggagggggggccggggggg                                                                                                                                 >  2:211855/1‑136 (MQ=255)                   tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                                                 <  1:682955/150‑1 (MQ=255)                          agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGATCGCAGAGAGGGATTATCt                                                                                                                          <  1:714359/150‑1 (MQ=255)                          agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCt                                                                                                                          <  1:247334/150‑1 (MQ=255)                           gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCTT‑TTTTCTTTTTCTGCCATGAGAAAGGAGCGCAGGTAGGGGTTATCTTg                                                                                                                        >  1:604573/1‑151 (MQ=255)                             ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGcc                                                                                                                      <  2:373840/151‑1 (MQ=255)                                               ggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTTCCATGAGAAGTGAGCGCCGAGAGGGGTTATTTTGCCTGA‑CTGAACTGGC‑GG‑ttt                                                                                                 >  2:264192/1‑150 (MQ=255)                                               ggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGGGAGGGCAGAGAGGGATTATCTTGCCTGA‑CTGAAATGGC‑GA‑TTa                                                                                                 >  2:458262/1‑151 (MQ=255)                                                 tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑g                                                                                              <  2:260265/151‑1 (MQ=255)                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCTT‑TTTTTTTTTTTTAGCATCAGAATTGGGCGTAGAGGGGGATTATTTTGCCTTATCTGAACggggggtaaaatt                                                                                               >  1:50793/1‑141 (MQ=255)                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑Gttt                                                                                           >  1:46073/1‑151 (MQ=255)                                                      cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGAGGGATTATTTTTCCTGA‑CTGAAATGGG‑GA‑TTAA‑Gtttt                                                                                          >  2:162098/1‑149 (MQ=255)                                                       ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CTTT‑TTTTTTTTTTTTGCCATGGGAAGGGAGCGCAGAGAGGGGTTATTTTTCCTGA‑CTTAACTGGC‑GC‑TTTA‑GTTTcc                                                                                         >  2:65891/1‑150 (MQ=255)                                                        tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTTCCATGAGAAATGAGCGCAGAGAGGGGTTATCTTTCCTGA‑ATGAAATGGG‑GG‑TTAA‑GTTTCcg                                                                                        >  2:181503/1‑150 (MQ=255)                                                             gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCTT‑TTTTTTTTTTTTGCCCTGAGGAGGGGGCTGCGAGGGGGGGTATCTTTTCTTA‑CTGGAATGGC‑GG‑TTAA‑GTTTTCGCGtcg                                                                                   >  1:495472/1‑148 (MQ=255)                                                             gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCTT‑TTTTTTTTTTTTGCCCTGAGAAGGGGGCGCAGAGAGGGATTATCTTTCCTTA‑CTGGACTTGG‑GATTTAA‑GTTCCcgggaa                                                                                    >  1:460543/1‑146 (MQ=255)                                                             gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTGGCCTGGGAAAGGAGGGCAGAGAGGGATTTTTTTTCCCTG‑CTGAACTGGCGGT‑TTAA‑GTTT‑CGCGTAc                                                                                   >  1:652386/1‑150 (MQ=255)                                                             gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTGCCATGAGAAGTGAGCGCCGAGAGGGGTTATCTTGCCTTG‑CTGAACTGGG‑GA‑TTTA‑GTTTTCGCGTAc                                                                                   >  1:164316/1‑150 (MQ=255)                                                                ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCTT‑TTTTTTTTTTTTGCCCTGGGAAGGGAGCGCCGAGAGGGGTTATTTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCCCGTACCCgg                                                                               >  2:792714/1‑151 (MQ=255)                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑TCTT‑TTTTTTTTTTTCGCCATTACGAGAGCCGGCAGGGGGGGGTTTTCTTTCCTGG‑CTCAAATGGC‑GC‑TTAT‑GTCTCCGCGTACCCGGCGt                                                                            >  1:143090/1‑151 (MQ=255)                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCTT‑TTTTTTTTTTCTTCCATGAGAAGGGAGCGCAGAGGGGGGTTATCTTGCCTGA‑CCGAACTTGCGGT‑TTAA‑GTTT‑CGCGGACCCGGCGt                                                                            >  2:33638/1‑151 (MQ=255)                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTGTCCTGAGAAGTGAGCGCAGGGAGGGGTTATCTTGCCTTA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCTCGTACCCGGGGt                                                                            >  1:218066/1‑151 (MQ=255)                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTTTGCCATTAGAAAGGAGCGCAGGGAGGGATTATCTTTCCTGA‑CTGAACTGGC‑GA‑TTAA‑ATTTCCGCGTACCCGGCgg                                                                            >  2:600055/1‑150 (MQ=255)                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGC‑CATT‑TTTTTTTTTTTTGTCCTGGGGAGAGAGCGCAGAGGGGGCATATTTTGCCTGA‑CTGCACTGGG‑GA‑TTTA‑GGTTCCGCGTACCCgtggt                                                                            >  1:465206/1‑147 (MQ=255)                                                                    aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTGCCCTGAGAAGGGAGAGCAGGGAGGGGTTTTTTTGCCTGA‑CTGAAATGGC‑GG‑TTAA‑ATTTTCGCGGTCCCGGCGtt                                                                           >  2:525736/1‑151 (MQ=255)                                                                    aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTCTGCCATGAGGAGGGGGCGCCGAGGGGGGTTTTCTTTCCCGA‑CTGAAATGGC‑GG‑TTAA‑AGTTTCGCGTACCCGGCGtt                                                                           >  2:31330/1‑151 (MQ=255)                                                                    aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTTTGCCATGAGAAGTGGGCGCAGAGAGGGGTTATCTTGCCTGA‑CTGAACTGGC‑GG‑TTAA‑GTTTTCCCGTACCCCGCGtt                                                                           >  2:141966/1‑151 (MQ=255)                                                                        ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGGGGGGGTATCTTGCCTGA‑CTTAACTGGC‑GA‑TTAA‑ATTTCCGCGTACCCCGCGtttgtt                                                                       >  1:320894/1‑148 (MQ=255)                                                                         taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGTT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGtt                                                                      <  1:721803/151‑1 (MQ=255)                                                                           ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGC‑TCTT‑TTTTTTTTTTTTGCCTTGAGAAGAGCGGGCAAAGGGGGGTTATCTTGCCTGT‑TTGTACTGGC‑GA‑TTAAGGTTTCCCCGTCCCCGCCGttattttt                                                                     >  2:105163/1‑145 (MQ=255)                                                                               tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTg                                                                  <  1:93361/149‑1 (MQ=255)                                                                                    aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatg                                                            <  2:210465/150‑1 (MQ=255)                                                                                              ggATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTCTGCCCTGAGAAAGGAGGACAGAGAGGGGTTATTTTGTCTGA‑CTGAAATGGC‑GG‑ATAA‑ATTTCCGCGTACCCGGCGGTTAGTTTTTGCGGGTGCTGATATTgta                                                 >  1:765133/1‑149 (MQ=255)                                                                                               gATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGa                                                  <  1:428231/149‑1 (MQ=255)                                                                                               gATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGa                                                  <  1:294935/149‑1 (MQ=255)                                                                                                aTGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTTTTTTTCTGTCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTTA‑CTGTACTGGC‑GG‑TTAA‑GTTTTCGCGTACCCCGCGGTTTGTTTTTGCGGATGATGGTATTTAAg                                                >  1:87062/1‑150 (MQ=255)                                                                                                aTGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTTTCTTGCCTGA‑CTGAAACGGG‑GA‑TTAA‑ATTTTCCCGTTCCCCGCGTTTAGTTTTTGGGGGTGGTGATATTTAAgt                                               >  2:581265/1‑150 (MQ=255)                                                                                                aTGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                                <  2:24640/150‑1 (MQ=255)                                                                                                aTGCCATGTCGCGGATCTCGTTTGCGGCGGG‑CATT‑TTTTTTTTTTTTGCCACGGGGAAGAAGCGGCGGGGGGGGTTATCTTTTCTGA‑ATGAAACGGT‑GA‑TTTA‑GTTTCCGGGGACCCCGCGGTTTGTTATTTGGGAGGATGAtttttgag                                                >  1:464904/1‑146 (MQ=255)                                                                                                   ccATGTCGCGGATCTCGTTTGCGGCGGT‑CCTT‑TTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGGGAGGGGTTATTTTGCCTTA‑CTGAAATGGG‑GG‑TTTA‑GTTTTCGCGGAACCGGGGGTTAGTTATTGCGGGTTATGGTAttttaagttg                                            >  1:221452/1‑143 (MQ=255)                                                                                                     aTGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                          <  2:296516/151‑1 (MQ=255)                                                                                                         cgcgGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGGGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAt                                      <  2:226793/151‑1 (MQ=255)                                                                                                               tctcGTTTGCGGCGGT‑CATT‑TTTTTTTTTTTTTCCATGAGAAGTGAGGGCAGAGAGGGGTTAtt                                                                                                                           >  2:681270/1‑63 (MQ=255)                                                                                                               tctcGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATc                                                                                                                           <  1:680964/64‑1 (MQ=255)                                                                                                               tctcGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTc                                <  1:422837/151‑1 (MQ=255)                                                                                                               tctcGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTc                                <  2:655659/151‑1 (MQ=255)                                                                                                                    tttGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGt                           <  1:507702/151‑1 (MQ=255)                                                                                                                    tttGCGGCGGT‑CATT‑TTTTCTTTTGCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGt                           <  1:342302/151‑1 (MQ=255)                                                                                                                                       ttctttttctGCCATGAGAAGTGAGCTCAGAGAGGGATTATCGTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc          <  1:336329/151‑1 (MQ=255)                                                                                                                                       ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc          <  1:69743/151‑1 (MQ=255)                                                                                                                                       ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc          <  1:734428/151‑1 (MQ=255)                                                                                                                                            tttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgtgt     <  2:132245/151‑1 (MQ=255)                                                                                                                                               ctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAtt  <  2:675357/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                           CGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATT  >  NC_000913/4604210‑4604500

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