Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,883	0	C	T	43.9%	45.7 / 51.9	68	intergenic (+52/‑51)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base C (10/27);  new base T (29/0);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.49e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACG  >  NC_000913/2386735‑2387032                                                                                                                                                        |                                                                                                                                                              aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTccc                                                                                                                                                             <  2:692773/151‑1 (MQ=255)     cACTAGCCCGGTGATACCGCAGACAGATTGAA‑TATAAGAAAAAGGATCAGGTAATCGATGGTGTTATAACTTTGTTAGACAACATTACATAATAAAAAGTTCTCGAGGCGCT‑TCCAGGAAAAGCGCCATGCA‑TATT‑TTAATTGGCCTCATc                                                                                                                                                        <  2:386328/151‑1 (MQ=255)           cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCCAAAAAACAACATGCC‑CATT‑TTTTTTTTCCCCCTCTtttt                                                                                                                                                   >  1:607228/1‑151 (MQ=255)               gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAAAAATAAAATTTTTCGGGGGCGGGCCCCAAAAAAAGAGA‑ATACA‑AATT‑TTTTTTTtccctttttttttttaa                                                                                                                                               >  1:449572/1‑137 (MQ=255)                tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑a                                                                                                                                             <  2:725486/151‑1 (MQ=255)                tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑a                                                                                                                                             <  2:400212/151‑1 (MQ=255)                    aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TTCAGAAAAAAAGACAAGCA‑CATT‑TTTTTTTTTCTCATCTTTTTTTAG‑ATGC‑g                                                                                                                                        >  2:733790/1‑151 (MQ=255)                    aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAAAGACACGCA‑CAAT‑TTTTTTTTCCTCCTCTGtttttaaaaggg                                                                                                                                          >  1:142143/1‑143 (MQ=255)                    aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCACAAAAAACAACAATCACCATT‑TTTTTTTTCC‑CCTCTGtttttttcagc                                                                                                                                           >  2:285707/1‑143 (MQ=255)                          gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAAAGACATGCC‑CATT‑TTTTTTTTCCTCTTCTGtttttaaaatggaaaaat                                                                                                                                    >  1:1075203/1‑137 (MQ=255)                              gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAACGAAAAGGA‑CATT‑TTTTTTTTCCTCACCTGTTTTTAG‑AAGG‑G‑aaaaaatgaa                                                                                                                             >  1:971512/1‑144 (MQ=255)                                 tGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTc                                                                                                                      <  1:39959/151‑1 (MQ=255)                                            aaaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACggagg                                                                                                             <  1:881955/148‑1 (MQ=255)                                                 ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATGAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAAGACAAGCA‑CAAT‑TTATTTTTTCTCCTCTGTTTTTAA‑ATGC‑G‑GAAC‑A‑T‑G‑AAAATT‑AAACGAGGaaaaaaaa                                                                                                     >  2:457832/1‑151 (MQ=255)                                                 ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAAAGACAAGGC‑CAAT‑TTTTTTTTCCTCCTCTTTTTTTAG‑ATGCGG‑AAAA‑A‑G‑A‑AAATTA‑AAGGGGGGaaaaaaa                                                                                                      >  1:529184/1‑149 (MQ=255)                                                  gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GGCAAAAAAAACCAAAAGCC‑CAAT‑TTTTTTTTTCTCATCTGTTTTTAGAAGGG‑G‑AAAA‑A‑T‑G‑AAATT‑‑AAAGGAGGGAaaataaa                                                                                                    >  2:798359/1‑151 (MQ=255)                                                       ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAACGACAAGCA‑CATT‑TTTTTTTTCCCCCTCTTTTTTTTA‑ATGC‑G‑AAAA‑A‑A‑T‑AGATTT‑AACCGAGGGAAAAAAAAAGGGa                                                                                               >  1:730843/1‑151 (MQ=255)                                                       ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCAGGAAAAAAGGCAAGCA‑CATT‑TTTTTTTTTCTCCTCTTTTTTTAA‑AAGG‑G‑GAAC‑A‑A‑A‑AAAATT‑AACGGGGGGAAAAAAAATgggg                                                                                               >  1:370490/1‑150 (MQ=255)                                                       ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCAAAAAAAAAGACAAGCA‑CATT‑TTTTTTTTCCTCCTCTTTTTTTAA‑ATGC‑G‑GAAA‑A‑A‑G‑AAAATC‑AACGGGGGGaaaaaaaaaagga                                                                                               >  1:113103/1‑146 (MQ=255)                                                       ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGGAAATAATAAAAATTTCTCGAGGCGGG‑TCCGGGAAAAACCAAAAACA‑CCAT‑TTATTTTTTCTCTTCTGTTTTTAA‑ATTCGG‑AAAC‑A‑‑‑C‑AAATTT‑TAAAGAGGGAAAATAAATgggg                                                                                               >  2:670129/1‑150 (MQ=255)                                                             ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAACGA‑AAGCA‑AATT‑TTATTTTTTCTTCTCTGGTTTTAG‑AAGG‑G‑AAAA‑A‑TAG‑AAATTA‑AAGCGGGGAAAATAAAATGGGatggtt                                                                                          >  2:290570/1‑146 (MQ=255)                                                             ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCTAAAAAAAAACATGCA‑AATTATTATTTT‑CCTCCTCTGTTTTTAA‑ATGC‑GAAAAC‑A‑T‑A‑AAATCA‑AAAGGGGGAAAAAAAAAAGGGagtaat                                                                                          >  2:814978/1‑146 (MQ=255)                                                                  gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCGGAAAAAAAGGCATGCA‑CCTT‑TTTTTTTTCCCCTTCTGTTTTTAG‑ATGG‑G‑TAAA‑A‑T‑GAAAATTC‑AAGGGGGGAAAAA‑AAATGGGGATGaattggta                                                                                    >  1:58548/1‑144 (MQ=255)                                                                  gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAAAGACATGAA‑CAAT‑TTTTTTTTCCCCATCTGTTTTTAG‑ATGC‑G‑GAAA‑A‑A‑G‑AAATTT‑AACCGAGGGAAAATAAATTGGAATGAGTGGGTa                                                                                    >  1:966249/1‑151 (MQ=255)                                                                  gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAACGACATGCA‑CATT‑TTTTTTTTCCTCCTCTGTTTTTAG‑ATGG‑G‑AAAA‑A‑T‑G‑AAAATT‑AAAGGGGGAAAAAAAAATGGGGAAGATTgggga                                                                                    >  1:964180/1‑148 (MQ=255)                                                                  gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAAGGCAATAC‑CATT‑TTTTTTTTCCTCCTCTGTTTTTAA‑AAGC‑G‑AAAA‑A‑A‑A‑AAAATT‑AAAGGAGGAAAAAAAAATAGGGAAGATTgggtt                                                                                    >  1:1121599/1‑148 (MQ=255)                                                                  gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAAGAAATGCA‑CATT‑TTTTTTTTCCTCATCTGTTTTTAA‑ATGC‑G‑AAAA‑A‑A‑G‑AAAATT‑AAAGGAGGGAAAAAAAAAGGGGAAGATTGGGTa                                                                                    >  1:1103162/1‑151 (MQ=255)                                                                  gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAAAACAAGCA‑AATT‑TTATTTTTTCTCATCTGTTTTTAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTT‑AAAGGAGGGAAAATAAATGGGAATGAATgggtt                                                                                    >  2:1081407/1‑148 (MQ=255)                                                                  gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCGAAAAAAACGGCAAGGA‑CAAT‑TTTTTTTTCCTCATCTGTTTTTTG‑AAGC‑G‑CAAC‑A‑T‑T‑AAATTC‑AACTGGGGGAAAAAAAATGGGAATGATTgggtt                                                                                    >  1:961643/1‑148 (MQ=255)                                                                   tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TTCCGGAAAAACGACATGCC‑AATT‑TTTTTTTTTCTCCTCTGTTTTTAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAAATA‑AAAGGAGGAAAAAAAAATAGGAATGGTTGGCTAt                                                                                   >  2:112470/1‑151 (MQ=255)                                                                    gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GGCACAAAAAAAAACATGCA‑CCAA‑TTTTTTTCCCTCATCTTTTTTTAA‑AAAC‑G‑AAAC‑A‑A‑C‑AAAATT‑AAAAGAGGGCAAATAAATTGGAATGATTGGTTAtt                                                                                  >  2:184781/1‑151 (MQ=255)                                                                             ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAACGACAAGCA‑CATT‑TTTTTTTTCCCCATCTGTTTTTAA‑ATGC‑G‑AAACAA‑T‑A‑AAATTA‑AAC‑GAGGAAAAAAAAAAGGGGATGATTTGCTTTTTTTCAGaaa                                                                         >  1:932329/1‑151 (MQ=255)                                                                               gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAAACCCATTCA‑CCTT‑TTTTTTTTCCTCATCTGTTTTTAA‑ATGC‑G‑AAAA‑A‑T‑G‑AAAATT‑AACCGGGGGAAAATAAAAGGGGATGTTTTGCTTTTTTTTAGaaaat                                                                       >  1:116980/1‑149 (MQ=255)                                                                               gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAACGACACACA‑AACT‑TTATTTTTCCCCATCTGTTTTTAG‑ATGC‑G‑TAAA‑A‑A‑A‑AAATTT‑AAATGAGGTAAAAAAAAAGGGGAAGAATGGGTGTTTTTCAGaaaat                                                                       >  1:952147/1‑149 (MQ=255)                                                                               gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCAGAAAAAACGAAATGCA‑CAAT‑TTTTTTTTCCCCATCTGTTTTAAA‑ATGG‑G‑AAAA‑A‑T‑GAAAATTA‑AACGAGGGAAAAAAAAATGGGAA‑GGTTGGGTATTTTTTAGaaaat                                                                       >  1:788413/1‑149 (MQ=255)                                                                                ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑GCCCGAAAAAACGGCAAGCA‑CATT‑TTTTTTTTCCTCATCTTTTTTAAG‑ATGC‑G‑AAAA‑A‑T‑T‑AAATTC‑AACGGGGGGAAAAAAAATGGGGATGATTTGCTATTTTTCAGAAAATg                                                                      >  1:668887/1‑151 (MQ=255)                                                                                     caacaGTAAATAATAAAAATTTCTCGAGGCGGG‑GTCACAAAAAAAGAAATGAA‑CATT‑TTTTTTTTTCTCATCTGTTTTTAG‑ATGC‑G‑AAAA‑A‑T‑A‑AAAATCTAACGG‑GGAAAAAAAAAAGGGAAAGATTGTGTATTTTTTAGAAAATTACTCa                                                                 >  2:658278/1‑151 (MQ=255)                                                                                       acaGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCCGAAAAAAAGACATGCA‑CATT‑TTTTTTTTCCTCCTTTTTTTTTAA‑ATGG‑G‑AAAA‑ATT‑A‑AAATTA‑AAGGGGGGAAAAAAAAATGGGAA‑GATTTGGTATTTTTTAGAAAATTACttaaa                                                               >  1:1150928/1‑147 (MQ=255)                                                                                           tAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaaa                                                             <  1:600813/149‑1 (MQ=255)                                                                                                   aaaaaTTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCa                                                   <  2:788063/151‑1 (MQ=255)                                                                                                     aaaTTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAAc                                                 <  1:113521/151‑1 (MQ=255)                                                                                                     aaaTTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAAc                                                 <  1:986581/151‑1 (MQ=255)                                                                                                                 gCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                                       <  2:67765/149‑1 (MQ=255)                                                                                                                 gCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                                       <  2:1102662/149‑1 (MQ=255)                                                                                                                 gCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                                       <  1:931922/149‑1 (MQ=255)                                                                                                                 gCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                                       <  1:343474/149‑1 (MQ=255)                                                                                                                 gCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                                       <  1:1108341/149‑1 (MQ=255)                                                                                                                  cGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAc                                    <  2:657426/151‑1 (MQ=255)                                                                                                                    gg‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa                                     >  2:1070839/1‑148 (MQ=255)                                                                                                                    gg‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa                                     >  2:788529/1‑148 (MQ=255)                                                                                                                           gAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGaa                                                                                              <  1:578949/85‑1 (MQ=255)                                                                                                                           gAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGaa                                                                                              >  2:578661/1‑85 (MQ=255)                                                                                                                                  cGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAACCCTTTGAt                     <  2:189262/151‑1 (MQ=255)                                                                                                                                  cGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt                     >  2:978438/1‑151 (MQ=255)                                                                                                                                    acttGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATgg                   >  2:1034832/4‑151 (MQ=255)                                                                                                                                     cATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGTTTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGa                  >  1:247642/1‑151 (MQ=255)                                                                                                                                            caTT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  1:222683/151‑1 (MQ=255)                                                                                                                                            caTT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            >  2:832143/1‑151 (MQ=255)                                                                                                                                            caTT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  1:55913/151‑1 (MQ=255)                                                                                                                                            caTT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  2:78118/151‑1 (MQ=255)                                                                                                                                            caTT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            >  2:289402/1‑151 (MQ=255)                                                                                                                                            caTT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  2:675035/151‑1 (MQ=255)                                                                                                                                              tt‑ttattttCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAGTGAAAAACCTTTGATGGATAATATAc          <  2:532880/151‑1 (MQ=255)                                                                                                                                              tt‑ttattttCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc          <  2:248649/151‑1 (MQ=255)                                                                                                                                              tt‑ttattttCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc          <  1:80018/151‑1 (MQ=255)                                                                                                                                              tt‑ttattttCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc          <  1:444139/151‑1 (MQ=255)                                                                                                                                              tt‑ttattttCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc          <  2:425540/151‑1 (MQ=255)                                                                                                                                              tt‑ttattttCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc          <  2:457786/151‑1 (MQ=255)                                                                                                                                                    atttCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATGC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGac     <  2:35612/150‑1 (MQ=255)                                                                                                                                                       tCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACg  >  2:1122643/1‑151 (MQ=255)                                                                                                                                                        cccTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGAAGAAATTGACTCAGAAAAGATTAAGCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACg  >  2:902469/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                              ATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGG‑TCCAGAAAAAGCGACATGCA‑CATT‑TTATTTTCCCTCATCTGTTTTAAG‑ATGC‑G‑AAAC‑A‑T‑G‑AAATTC‑AACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACG  >  NC_000913/2386735‑2387032

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