Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,862	0	G	A	57.0%	‑4.3 / 98.0	81	intergenic (+31/‑72)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (5/29);  new base A (45/0);  total (52/29)
Fisher's exact test for biased strand distribution p-value = 8.36e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-12

GTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAA  >  NC_000913/2386725‑2387003                                                                                                                                           |                                                                                                                                                            gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TACAGCCAAAAAC‑ACA‑TGA‑A‑Catt                                                                                                                                            >  1:1037017/1‑151 (MQ=255) gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCAGAAAAAACG‑ACA‑TGC‑A‑Catt                                                                                                                                            >  2:185766/1‑151 (MQ=255)      aaaCAATNCCACGTGCCCGGTGATACCGCNNNNNNATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑Cattttatt                                                                                                                                       <  1:266514/151‑1 (MQ=255)      aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TGCAGAACAAAAG‑ACA‑TAC‑A‑Aactttttt                                                                                                                                       >  1:800707/1‑151 (MQ=255)      aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAAG‑ACA‑AGA‑A‑Cattttttt                                                                                                                                       >  1:657638/1‑151 (MQ=255)      aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCAGAAAAAACGCACA‑‑AC‑A‑Cattttatt                                                                                                                                       >  1:252928/1‑151 (MQ=255)      aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGC‑GCCAGGAAAACCG‑ACA‑AAC‑A‑Cattttttt                                                                                                                                       >  1:1253829/1‑151 (MQ=255)        aCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑Cattttatttt                                                                                                                                     <  2:96983/151‑1 (MQ=255)         cAATACCACTTGCCCGGTGATACCGCGGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TGCAGCAAAAGAC‑ATA‑TCC‑C‑Catttttttttt                                                                                                                                    >  1:1333403/1‑150 (MQ=255)         cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCAAAAAAACG‑AAA‑AGC‑A‑Caattttttttt                                                                                                                                    >  2:60060/1‑150 (MQ=255)            tACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCt                                                                                                                                 <  2:1461320/151‑1 (MQ=255)                aCTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATc                                                                                                                             <  2:677188/151‑1 (MQ=255)                aCTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATc                                                                                                                             <  2:34730/151‑1 (MQ=255)                aCTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATc                                                                                                                             <  1:1247818/151‑1 (MQ=255)                     cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAAAAAAAAGG‑GCA‑TGC‑ACCATTTTTTTTTccccctcttttt                                                                                                                         >  1:1026472/1‑142 (MQ=255)                          tGATACCGCAGACATATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGa                                                                                                                 <  1:1102550/151‑1 (MQ=255)                          tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGa                                                                                                                 <  1:116485/151‑1 (MQ=255)                          tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGa                                                                                                                 <  1:1214471/151‑1 (MQ=255)                          tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGa                                                                                                                 <  2:1026980/151‑1 (MQ=255)                          tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGa                                                                                                                 <  1:929809/151‑1 (MQ=255)                          tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGa                                                                                                                 <  2:801162/151‑1 (MQ=255)                                  cagacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCAGAAAAAAAG‑ACA‑TGC‑A‑AATTTTTTTTTTCCTCATC‑TG‑tttttaaaagggaaaa                                                                                                         >  2:1241075/1‑139 (MQ=255)                                      cagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAAAAAAAAAA‑AAA‑AAC‑A‑CATTTTATTTTTTCTTATC‑TC‑TTTTacaaatctaaacaaga                                                                                                     >  1:378768/1‑136 (MQ=255)                                        gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TGCCAAAAAAACG‑ACA‑AGC‑A‑AATTTTTTTTTTCCTCCAC‑TG‑TTTTTAGaattggaaaaacaaa                                                                                                   >  2:533087/1‑137 (MQ=255)                                         aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑At                                                                                                    <  2:456506/146‑1 (MQ=255)                                              aTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAc                                                                                         <  2:636725/151‑1 (MQ=255)                                                     gAAAATGGCTCAGGAAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TTCAGAGAAAACG‑AAA‑AGC‑A‑CATTTTATTTTTCCTCCTT‑TT‑TTTTAAAATT‑C‑GAAAA‑ATT‑aaaattaaaaggggaa                                                                                  >  2:1345775/1‑138 (MQ=255)                                                           ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCAAAAAAACG‑GAA‑TGC‑A‑CATTTTTTTTTTTCTCATC‑TG‑TTTTTAAATG‑C‑GAAAC‑ATGAAAATTTAAG‑GGG‑GG‑aaaaaaa                                                                          >  2:967416/1‑149 (MQ=255)                                                           ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCAGAAAAAAAG‑ACA‑AGC‑A‑CATTTTTTTTTTCCTCCTC‑TG‑TTTTTAAAAG‑C‑GAAAA‑ATG‑AAAATTAAC‑GGA‑GG‑GAaaaaaa                                                                         >  2:168639/1‑151 (MQ=255)                                                            gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCAGAAAAAACG‑ACA‑AGG‑A‑CATTTTATTTTTCCCCATC‑TT‑TTTTTAAAAG‑G‑GAAACAAAA‑AAATTAAAC‑GGG‑GG‑aaaaaaaa                                                                         >  2:1075507/1‑151 (MQ=255)                                                              tCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATg                                                                      <  1:888235/151‑1 (MQ=255)                                                                aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGGAAAAAAC‑‑AAA‑TGT‑A‑CCTTTTATTTTTCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAA‑ATG‑AAAATCAAC‑GGA‑GG‑AAAAATAAATagg                                                                    >  1:1279399/1‑151 (MQ=255)                                                                 ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAACG‑ACA‑TGC‑A‑AATTTTTTTTTTCCTCCTC‑TC‑TTTTTAAATG‑CGGAAAA‑AGA‑AAATTAACC‑GGG‑GG‑AAAAAAAAAgggg                                                                    >  2:194253/1‑151 (MQ=255)                                                                 ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAAAA‑ACA‑AGC‑C‑AAATTTATTTTTCCTCCTT‑CG‑TTTTTAAAAG‑C‑GGAAA‑AAG‑AAAATCAAC‑TGT‑GG‑AAAAATAAATgggg                                                                   >  1:130190/1‑150 (MQ=255)                                                                 ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAAG‑ACA‑TGC‑A‑CATTTTTTTTTTCCTCATC‑TG‑TTTTAAGAAG‑T‑GAAAA‑AAG‑AAAATTAAA‑GGG‑GG‑AAAAAAAAAAGGGa                                                                   >  2:900526/1‑151 (MQ=255)                                                                 ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAAAG‑ACA‑AGC‑A‑CATTTTTTTTTTTCTCATC‑TG‑TTTTTAGATG‑C‑GGAAC‑AAG‑AAAATCAAA‑CGA‑GG‑GAAAAAAAATtggg                                                                   >  1:1373674/1‑150 (MQ=255)                                                                  gTAGTGGGTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAAAG‑AAA‑AGG‑A‑CCTTTTTTTTTTCCTCATT‑TG‑TTTTTAAAAG‑C‑GGAAC‑AAG‑AAAATTAAA‑GGA‑GG‑AAAAAAAAATgggga                                                                  >  2:480455/1‑149 (MQ=255)                                                                  gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TTCACGAAAAAAG‑ACA‑TGG‑A‑AATTTTTTTTTTCCTCATC‑TT‑TTTTTAAATG‑C‑GGAAA‑ATG‑AAATTTAAA‑GGA‑GG‑AAAAATAAATTGGaa                                                                  >  1:362804/1‑151 (MQ=255)                                                                  gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TTCAAAAAAAAAG‑AAA‑AGG‑A‑CAATTTTTTTTTCCTCCTCTTT‑TTTTTAAAAG‑G‑AAAAA‑AGA‑AAATTAAAG‑GGG‑GG‑AAAAAAAAATGGGa                                                                   >  1:237116/1‑151 (MQ=255)                                                                  gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAACG‑ACA‑AGC‑A‑CATTTTTTTTTTCCTCCTC‑TT‑TTTTTAAAAG‑C‑GAAAA‑AAGAAAATTAAAC‑GGG‑GG‑AAAAAAAAAGGGGa                                                                   >  2:1246112/1‑151 (MQ=255)                                                                  gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAAAG‑ACA‑TGC‑A‑CATTTTTTTTTTCCTCCTC‑TT‑TTTTTAGATT‑C‑GAAAA‑AAG‑AAAATTAAA‑GGG‑GG‑GAAAAAAAATgggga                                                                  >  2:1337880/1‑149 (MQ=255)                                                                  gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAACG‑GCA‑AGC‑C‑CATTTTTTTTTTTCCCATC‑TT‑TTTTTAAAAG‑G‑GGAAA‑ATT‑AAAATTAAC‑GGAGGG‑AAAAAAAAATGGGa                                                                   >  2:203931/1‑151 (MQ=255)                                                                   tAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAAt                                                                 <  1:1320159/151‑1 (MQ=255)                                                                     gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TTCCGAAAAAAAG‑ACA‑TGC‑A‑CATTTTATTTTTTCTCATC‑TG‑TTTTTAGAAG‑C‑GGAAA‑AAG‑AAAATTAAA‑GGA‑GG‑GAAAAAAAAAGGGAATG‑a                                                              >  1:887938/1‑151 (MQ=255)                                                                      tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑At                                                             <  2:498903/151‑1 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGTGGCGGG‑TGCCGGCAAAATG‑ACA‑TGC‑A‑CAATTTTTTTTCCCTTCTC‑TG‑TTTTAAAAAA‑C‑GAAAC‑ATA‑AAATTCAAC‑TGA‑GG‑AAAAATAAATTGGAATG‑Atttgg                                                         >  1:1466025/1‑146 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGT‑TTCCGAAAAAAAG‑AAA‑AAA‑A‑AATTTTATTTTTTCTCAAC‑TG‑TTTTTAAAAG‑C‑GAAAA‑AAG‑AAAATTAAC‑AGA‑GGCAAAAAAAAATGGG‑ATT‑ATTGGCtt                                                       >  1:775479/1‑150 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGCCCCCCCAAATCCG‑AC‑‑TTC‑A‑CACTTTTTTTTTTCTCCAC‑TG‑TTTTAATATGCC‑GAACA‑ACA‑AAAATAAAC‑GCG‑GG‑ACAAA‑AAATTGGAATG‑Atttggc                                                        >  2:1042192/1‑146 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TTCAGAAAAAAAG‑AAA‑TGA‑A‑AAATTTTTTTTTTCTCCAC‑TG‑GTTTTAGAAA‑CGGAAAA‑AGA‑AAATTCAAC‑‑GA‑GG‑GAAAAAAAATGGGGATT‑AATGGGTa                                                       >  1:313837/1‑151 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAACG‑ACA‑TGC‑A‑CATTTTTTTTTTCCCCCTT‑TT‑TTTTTAGATG‑G‑GGAAC‑AAG‑AAATTCAAC‑CGG‑GG‑AAAAAAAAATGGGGATG‑ATTgggtt                                                       >  2:802059/1‑148 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAACG‑GAA‑TGC‑A‑CAATTTTTTTTTTCCCATC‑TT‑TTTTTAAATG‑C‑GGAAC‑AAG‑AAAATCAAA‑GGG‑GG‑GAAAAAAAATAGGGATG‑ATTTGCTa                                                       >  2:1204909/1‑151 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAAC‑ACA‑AGC‑C‑CAATTTATTTTTCCCCACC‑TGCTTTTATAATT‑C‑AAAAA‑AT‑‑AAATTTAAA‑GGA‑GG‑GAAAATAAAAGTGGATG‑agttggta                                                       >  1:143391/1‑144 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GTCAGAAAAAAAG‑ACACCGC‑A‑AATTTTTTTTTCCTTCA‑C‑TG‑TTTTTAAAAG‑A‑GGAACAATA‑AAAATAAAG‑GGG‑GA‑AAAAAAAAATGGGATTG‑TTTgcaa                                                        >  1:491491/1‑148 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GTCAGAAAAAAAG‑ACA‑AGA‑A‑AATTTTATTTTTTCTCATC‑TG‑TTTTAAGATG‑A‑GAAAC‑AAG‑AAAATCAAC‑GGA‑AG‑AAAAAAAAAAGGGAAAG‑ATTggtgt                                                       >  1:1160420/1‑148 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GTCAGAAAAAAAG‑AAA‑AGA‑A‑AATTTTATTTTTTCTCATC‑TG‑TTTTTAAATG‑C‑GAAAC‑AAG‑AAATTTAAA‑GGA‑GG‑GAAAATAAATGGGGATG‑ATTTGCTa                                                       >  1:466271/1‑151 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GGCCGAAAAAAAG‑ACA‑AGC‑A‑AATTTTTTTTTTCCTCCTC‑TG‑TTTTTAAATG‑C‑GGAAA‑ATG‑AAAATTAAA‑GGA‑GG‑GAAAAAAAATGGGAAAG‑ATTGGCTa                                                       >  2:1163170/1‑151 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAAAG‑ACA‑AGC‑A‑CATTTTTTTTTTCCTCATT‑TG‑TTTTTAGATG‑C‑GAAAA‑AAG‑AAATTTAAA‑GGA‑GG‑AAAAAAAAATGGGGATG‑ATTGGCTa                                                       >  2:1221259/1‑151 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑AACCGCAAAAAAG‑ACA‑AAG‑A‑AATTTTATTTTTTCTCATC‑TG‑TTTTTTGATG‑A‑GAAAA‑ATG‑AAAATCAAA‑GGA‑GG‑AAAAAAAAAAAGGAATG‑AATgggtt                                                       >  1:165409/1‑148 (MQ=255)                                                                              gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAACG‑ACA‑TGC‑A‑CATTTTTTTTTTCCTCCTC‑TG‑TTTTTAGATG‑C‑GGAAA‑AAG‑AAAATTAAC‑GGG‑GG‑GAAAAAAAATGGGAAAG‑ATTGGGTAtt                                                     >  2:1481007/1‑151 (MQ=255)                                                                                       ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCAAAAAAAACG‑ACA‑AGC‑A‑AAATTTTTTTTTCCTCATCTTT‑TTTTAAAATG‑G‑AAAAA‑AGA‑AAATTAAAG‑GGG‑GA‑AAAAAAAAATGGGATTA‑TTTGGTtttttttcgaaa                                             >  2:91022/1‑145 (MQ=255)                                                                                         gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAAAG‑ACA‑AGC‑C‑CATTTTTTTTTTCCCCCTT‑TT‑TTTTTAAATG‑G‑GGAAA‑ATG‑AAAATTAACAGGG‑GG‑AAAAAAAATTGGGAAGA‑TTTGGTTTTTTTGCGAAAAt                                           >  2:79813/1‑151 (MQ=255)                                                                                         gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑CTCCGAAAAAACG‑ACA‑TGCAA‑CATTTTTTTTTTCCCCTTT‑TT‑TTTTAA‑AAA‑C‑GCAAA‑ACG‑AAAATCCAA‑CCC‑GG‑GAAAATAAATTGGGATG‑ATTTGCCATTTTGTAGaaaat                                          >  2:1513543/1‑149 (MQ=255)                                                                                            ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAAAG‑ACA‑ATC‑A‑CATTTTTTTTTTTCTCATC‑CG‑TTTTTTAATG‑C‑GGAAA‑CTT‑AAATTAAAA‑GGA‑GG‑GAAAAAAAATAGGGATGAATTGGTttttttctcaaaatgaa                                        >  1:679940/1‑140 (MQ=255)                                                                                            ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGAc                                       <  1:14504/151‑1 (MQ=255)                                                                                            ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAAAG‑ACA‑AGC‑A‑CATTTTATTTTTTCTCATT‑TT‑TTTTTAGATG‑G‑GAAAA‑AAGAAAATTAAAA‑GGG‑GG‑AAAAA‑AAATTGGGATG‑ATTTGGTCTTTTTCAGAAATTGaa                                       >  1:601695/1‑150 (MQ=255)                                                                                              acaacaTTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GA‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTc                                     <  2:1124686/151‑1 (MQ=255)                                                                                               caacaGTAAATAATAAAAATTTCTCGAGGCGGG‑GGCAGAAAAAACG‑ACA‑TTC‑A‑CATTTTTTTTTTTCTCATC‑TT‑TTTTATAATG‑C‑GAAAC‑AAG‑AAAATCAAC‑GGAGGG‑AAAAAAAAATGGG‑ATG‑ATTTGCGTTTTTGTTGaaaaatactca                                    >  1:647561/1‑143 (MQ=255)                                                                                                        ataataAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGAtt                           <  1:1223251/151‑1 (MQ=255)                                                                                                        ataataAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGAtt                           <  2:1001213/151‑1 (MQ=255)                                                                                                           ataaAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAc                        <  1:1268944/151‑1 (MQ=255)                                                                                                            taaAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGAAG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAcc                       <  2:527750/151‑1 (MQ=255)                                                                                                               aaaTTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAAc                    <  1:654636/151‑1 (MQ=255)                                                                                                               aaaTTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAAc                    <  1:516473/151‑1 (MQ=255)                                                                                                                           gCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc          <  1:461588/149‑1 (MQ=255)                                                                                                                           gCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc          <  1:14065/149‑1 (MQ=255)                                                                                                                              gg‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa        >  2:206388/1‑148 (MQ=255)                                                                                                                              gg‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa        >  2:123611/1‑148 (MQ=255)                                                                                                                              gg‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa        >  1:1341187/1‑148 (MQ=255)                                                                                                                               g‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTg     <  2:1471948/150‑1 (MQ=255)                                                                                                                               g‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGa    >  1:415640/1‑151 (MQ=255)                                                                                                                                  ccAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaa  <  1:426238/151‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                            GTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCG‑ACA‑TGC‑A‑CATTTTATTTTCCCTCATC‑TG‑TTTTAAGATG‑C‑GAAAC‑ATG‑AAATTCAAC‑GGA‑GG‑AAAAATAAATGGGAATG‑ATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAA  >  NC_000913/2386725‑2387003

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