Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,344	0	C	T	57.9%	‑1.8 / 54.1	57	S62F (TCT→TTT)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (5/19);  new base T (33/0);  total (38/19)
Fisher's exact test for biased strand distribution p-value = 6.67e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.02e-06

ATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  NC_000913/4604199‑4604492                                                                                                                                                    |                                                                                                                                                         aTAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATt‑ttttcttttt                                                                                                                                                     >  1:675052/1‑151 (MQ=255)            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)                                                                                                                                                    |                                                                                                                                                         ATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATT‑TTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGA‑CTGAACTGGC‑GA‑TTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  NC_000913/4604199‑4604492

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