Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,862	0	G	A	61.5%	16.6 / 65.5	65	intergenic (+31/‑72)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (4/21);  new base A (40/0);  total (44/21)
Fisher's exact test for biased strand distribution p-value = 2.08e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-07

ACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC  >  NC_000913/2386712‑2386998                                                                                                                                                        |                                                                                                                                                    acCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAg                                                                                                                                                     <  2:1094926/151‑1 (MQ=255)       tttgATCGTATTCAACAAGTCCACTTGCCCTGTGATACCGCAGACAGATTGAATTATCAGAAACTGGATAAGGTAATGGAAGATGGTATGACTTTGGTATACAACAATAAATAATGAAAATGTGTCGAGGCGGT‑TCCAGAAACAGCGACac                                                                                                                                               <  2:834300/151‑2 (MQ=255)           aTCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGc‑a‑c                                                                                                                                         <  2:138777/151‑1 (MQ=255)              gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAAAAAAAACGACATGC‑A‑Catt                                                                                                                                      >  1:821739/1‑151 (MQ=255)              gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCCAAAAAAAGACaaggccaat                                                                                                                                        >  1:494716/1‑143 (MQ=255)                 ataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GTCAGAAAAAAAGAGAAAC‑A‑Catt‑ttt                                                                                                                                  >  2:881581/1‑150 (MQ=255)                   aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TTCGGAAAAAAAGACCTGCAA‑Catt‑tttt                                                                                                                                 >  1:535870/1‑149 (MQ=255)                   aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑Catt‑ttatt                                                                                                                                <  1:1099373/151‑1 (MQ=255)                    aaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑Catt‑ttattt                                                                                                                               <  2:175513/151‑1 (MQ=255)                     aCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TTCCGAAAAAACGAAAATG‑A‑Catt‑ttttttt                                                                                                                              >  1:506592/1‑151 (MQ=255)                      cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TTCCGAAAAAACGGAAAGC‑A‑Catt‑tttttttt                                                                                                                             >  2:79994/1‑150 (MQ=255)                      cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TGCAGCAAAAAAGACATGC‑A‑Caat‑tttttttt                                                                                                                             >  1:10820/1‑150 (MQ=255)                      cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCGCGAAAAAAGAACATGC‑A‑Caat‑tttttttt                                                                                                                             >  2:1085427/1‑150 (MQ=255)                      cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAAAGGCAAGG‑A‑Caat‑tttttttt                                                                                                                             >  1:67789/1‑150 (MQ=255)                        aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTccc                                                                                                                           <  2:692773/151‑1 (MQ=255)                            cACTAGCCCGGTGATACCGCAGACAGATTGAA‑TATAAGAAAAAGGATCAGGTAATCGATGGTGTTATAACTTTGTTAGACAACATTACATAATAAAAAGTTCTCGAGGCGCT‑TCCAGGAAAAGCGCCATGC‑A‑TATT‑TTAATTGGCCTCATc                                                                                                                      <  2:386328/151‑1 (MQ=255)                                  cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCCAAAAAACAACATGC‑C‑CATT‑TTTTTTTTCCCCCTCTtttt                                                                                                                 >  1:607228/1‑151 (MQ=255)                                      gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAAAAATAAAATTTTTCGGGGGCGGGCCCCAAAAAAAGAGA‑ATAC‑A‑AATT‑TTTTTTTtccctttttttttttaa                                                                                                             >  1:449572/1‑137 (MQ=255)                                       tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑a                                                                                                           <  2:400212/151‑1 (MQ=255)                                       tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑a                                                                                                           <  2:725486/151‑1 (MQ=255)                                           aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TTCAGAAAAAAAGACAAGC‑A‑CATT‑TTTTTTTTTCTCATCTTTTTTTAG‑ATGC‑g                                                                                                      >  2:733790/1‑151 (MQ=255)                                           aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAAAGACACGC‑A‑CAAT‑TTTTTTTTCCTCCTCTGtttttaaaaggg                                                                                                        >  1:142143/1‑143 (MQ=255)                                           aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCACAAAAAACAACAATC‑ACCATT‑TTTTTTTTCC‑CCTCTGtttttttcagc                                                                                                         >  2:285707/1‑143 (MQ=255)                                                 gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAAAGACATGC‑C‑CATT‑TTTTTTTTCCTCTTCTGtttttaaaatggaaaaat                                                                                                  >  1:1075203/1‑137 (MQ=255)                                                     gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAACGAAAAGG‑A‑CATT‑TTTTTTTTCCTCACCTGTTTTTAG‑AAGG‑G‑aaaaaatgaa                                                                                           >  1:971512/1‑144 (MQ=255)                                                        tGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTc                                                                                    <  1:39959/151‑1 (MQ=255)                                                                   aaaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACggagg                                                                           <  1:881955/148‑1 (MQ=255)                                                                        ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATGAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAAGACAAGC‑A‑CAAT‑TTATTTTTTCTCCTCTGTTTTTAA‑ATGC‑G‑GAAC‑A‑T‑G‑AAAATT‑AAACGAGGaaaaaaaa                                                                   >  2:457832/1‑151 (MQ=255)                                                                        ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAAAGACAAGG‑C‑CAAT‑TTTTTTTTCCTCCTCTTTTTTTAG‑ATGCGG‑AAAA‑A‑G‑A‑AAATTA‑AAGGGGGGaaaaaaa                                                                    >  1:529184/1‑149 (MQ=255)                                                                         gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GGCAAAAAAAACCAAAAGC‑C‑CAAT‑TTTTTTTTTCTCATCTGTTTTTAGAAGGG‑G‑AAAA‑A‑T‑G‑AAATT‑‑AAAGGAGGGAaaataaa                                                                  >  2:798359/1‑151 (MQ=255)                                                                              ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAACGACAAGC‑A‑CATT‑TTTTTTTTCCCCCTCTTTTTTTTA‑ATGC‑G‑AAAA‑A‑A‑T‑AGATTT‑AACCGAGGGAAAAAAAAAGGGa                                                             >  1:730843/1‑151 (MQ=255)                                                                              ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCAGGAAAAAAGGCAAGC‑A‑CATT‑TTTTTTTTTCTCCTCTTTTTTTAA‑AAGG‑G‑GAAC‑A‑A‑A‑AAAATT‑AACGGGGGGAAAAAAAATgggg                                                             >  1:370490/1‑150 (MQ=255)                                                                              ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCAAAAAAAAAGACAAGC‑A‑CATT‑TTTTTTTTCCTCCTCTTTTTTTAA‑ATGC‑G‑GAAA‑A‑A‑G‑AAAATC‑AACGGGGGGaaaaaaaaaagga                                                             >  1:113103/1‑146 (MQ=255)                                                                              ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGGAAATAATAAAAATTTCTCGAGGCGGG‑TCCGGGAAAAACCAAAAAC‑A‑CCAT‑TTATTTTTTCTCTTCTGTTTTTAA‑ATTCGG‑AAAC‑A‑‑‑C‑AAATTT‑TAAAGAGGGAAAATAAATgggg                                                             >  2:670129/1‑150 (MQ=255)                                                                                    ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAACGA‑AAGC‑A‑AATT‑TTATTTTTTCTTCTCTGGTTTTAG‑AAGG‑G‑AAAA‑A‑TAG‑AAATTA‑AAGCGGGGAAAATAAAATGGGatggtt                                                        >  2:290570/1‑146 (MQ=255)                                                                                    ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCTAAAAAAAAACATGC‑A‑AATTATTATTTT‑CCTCCTCTGTTTTTAA‑ATGC‑GAAAAC‑A‑T‑A‑AAATCA‑AAAGGGGGAAAAAAAAAAGGGagtaat                                                        >  2:814978/1‑146 (MQ=255)                                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCGGAAAAAAAGGCATGC‑A‑CCTT‑TTTTTTTTCCCCTTCTGTTTTTAG‑ATGG‑G‑TAAA‑A‑T‑GAAAATTC‑AAGGGGGGAAAAA‑AAATGGGGATGaattggta                                                  >  1:58548/1‑144 (MQ=255)                                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAAAGACATGA‑A‑CAAT‑TTTTTTTTCCCCATCTGTTTTTAG‑ATGC‑G‑GAAA‑A‑A‑G‑AAATTT‑AACCGAGGGAAAATAAATTGGAATGAGTGGGTa                                                  >  1:966249/1‑151 (MQ=255)                                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAACGACATGC‑A‑CATT‑TTTTTTTTCCTCCTCTGTTTTTAG‑ATGG‑G‑AAAA‑A‑T‑G‑AAAATT‑AAAGGGGGAAAAAAAAATGGGGAAGATTgggga                                                  >  1:964180/1‑148 (MQ=255)                                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAAGGCAATA‑C‑CATT‑TTTTTTTTCCTCCTCTGTTTTTAA‑AAGC‑G‑AAAA‑A‑A‑A‑AAAATT‑AAAGGAGGAAAAAAAAATAGGGAAGATTgggtt                                                  >  1:1121599/1‑148 (MQ=255)                                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAAGAAATGC‑A‑CATT‑TTTTTTTTCCTCATCTGTTTTTAA‑ATGC‑G‑AAAA‑A‑A‑G‑AAAATT‑AAAGGAGGGAAAAAAAAAGGGGAAGATTGGGTa                                                  >  1:1103162/1‑151 (MQ=255)                                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAAAACAAGC‑A‑AATT‑TTATTTTTTCTCATCTGTTTTTAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTT‑AAAGGAGGGAAAATAAATGGGAATGAATgggtt                                                  >  2:1081407/1‑148 (MQ=255)                                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCGAAAAAAACGGCAAGG‑A‑CAAT‑TTTTTTTTCCTCATCTGTTTTTTG‑AAGC‑G‑CAAC‑A‑T‑T‑AAATTC‑AACTGGGGGAAAAAAAATGGGAATGATTgggtt                                                  >  1:961643/1‑148 (MQ=255)                                                                                          tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TTCCGGAAAAACGACATGC‑C‑AATT‑TTTTTTTTTCTCCTCTGTTTTTAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAAATA‑AAAGGAGGAAAAAAAAATAGGAATGGTTGGCTAt                                                 >  2:112470/1‑151 (MQ=255)                                                                                           gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GGCACAAAAAAAAACATGC‑A‑CCAA‑TTTTTTTCCCTCATCTTTTTTTAA‑AAAC‑G‑AAAC‑A‑A‑C‑AAAATT‑AAAAGAGGGCAAATAAATTGGAATGATTGGTTAtt                                                >  2:184781/1‑151 (MQ=255)                                                                                                    ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAACGACAAGC‑A‑CATT‑TTTTTTTTCCCCATCTGTTTTTAA‑ATGC‑G‑AAACAA‑T‑A‑AAATTA‑AAC‑GAGGAAAAAAAAAAGGGGATGATTTGCTTTTTTTCAGaaa                                       >  1:932329/1‑151 (MQ=255)                                                                                                      gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAACCCATTC‑A‑CCTT‑TTTTTTTTCCTCATCTGTTTTTAA‑ATGC‑G‑AAAA‑A‑T‑G‑AAAATT‑AACCGGGGGAAAATAAAAGGGGATGTTTTGCTTTTTTTTAGaaaat                                     >  1:116980/1‑149 (MQ=255)                                                                                                      gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAACGACACAC‑A‑AACT‑TTATTTTTCCCCATCTGTTTTTAG‑ATGC‑G‑TAAA‑A‑A‑A‑AAATTT‑AAATGAGGTAAAAAAAAAGGGGAAGAATGGGTGTTTTTCAGaaaat                                     >  1:952147/1‑149 (MQ=255)                                                                                                      gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCAGAAAAAACGAAATGC‑A‑CAAT‑TTTTTTTTCCCCATCTGTTTTAAA‑ATGG‑G‑AAAA‑A‑T‑GAAAATTA‑AACGAGGGAAAAAAAAATGGGAA‑GGTTGGGTATTTTTTAGaaaat                                     >  1:788413/1‑149 (MQ=255)                                                                                                       ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAACGGCAAGC‑A‑CATT‑TTTTTTTTCCTCATCTTTTTTAAG‑ATGC‑G‑AAAA‑A‑T‑T‑AAATTC‑AACGGGGGGAAAAAAAATGGGGATGATTTGCTATTTTTCAGAAAATg                                    >  1:668887/1‑151 (MQ=255)                                                                                                            caacaGTAAATAATAAAAATTTCTCGAGGCGGG‑GTCACAAAAAAAGAAATGA‑A‑CATT‑TTTTTTTTTCTCATCTGTTTTTAG‑ATGC‑G‑AAAA‑A‑T‑A‑AAAATCTAACGG‑GGAAAAAAAAAAGGGAAAGATTGTGTATTTTTTAGAAAATTACTCa                               >  2:658278/1‑151 (MQ=255)                                                                                                              acaGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAAAGACATGC‑A‑CATT‑TTTTTTTTCCTCCTTTTTTTTTAA‑ATGG‑G‑AAAA‑ATT‑A‑AAATTA‑AAGGGGGGAAAAAAAAATGGGAA‑GATTTGGTATTTTTTAGAAAATTACttaaa                             >  1:1150928/1‑147 (MQ=255)                                                                                                                  tAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaaa                           <  1:600813/149‑1 (MQ=255)                                                                                                                          aaaaaTTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCa                 <  2:788063/151‑1 (MQ=255)                                                                                                                            aaaTTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAAc               <  1:113521/151‑1 (MQ=255)                                                                                                                            aaaTTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAAc               <  1:986581/151‑1 (MQ=255)                                                                                                                                        gCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc     <  1:343474/149‑1 (MQ=255)                                                                                                                                        gCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc     <  2:67765/149‑1 (MQ=255)                                                                                                                                        gCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc     <  1:931922/149‑1 (MQ=255)                                                                                                                                        gCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc     <  2:1102662/149‑1 (MQ=255)                                                                                                                                        gCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc     <  1:1108341/149‑1 (MQ=255)                                                                                                                                         cGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAc  <  2:657426/151‑1 (MQ=255)                                                                                                                                           gg‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa   >  2:788529/1‑148 (MQ=255)                                                                                                                                           gg‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa   >  2:1070839/1‑148 (MQ=255)                                                                                                                                                  gAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGaa                                                            >  2:578661/1‑85 (MQ=255)                                                                                                                                                  gAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGaa                                                            <  1:578949/85‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                    ACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGC‑A‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC  >  NC_000913/2386712‑2386998

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