Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,151,691	C→T	20.0%	intergenic (+45/+20)	yegJ → / ← yegK	uncharacterized protein/ser/thr phosphatase‑related protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,691	0	C	T	20.0%	64.9 / 3.5	43	intergenic (+45/+20)	yegJ/yegK	uncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base C (17/17);  new base T (4/5);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACC‑ATTGATCCTTTATTGGTCA‑CCT‑ATCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATG  >  NC_000913/2151496‑2151925                                                                                                                                                                                                              |                                                                                                                                                                                                                                               tGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAAC‑ATTTATTCCTTTTTTGTCACCCA‑ATCCAAAA‑AAGAAAGCGAC‑GAGC‑TAAAAGGAAAAATAATTTGCTTA‑TCAACC‑G‑GTTTATTTTTTTTTTTAAATTGGGGGTTTGA‑TTTTTTTT‑TT‑TAA‑CACCCCGTAAAAAACACACTTTCCCcacacccagcg                                                                                                                                                                                                    >  2:477713/1‑231 (MQ=255)               aaaaCAAAAAAGAAGATATTTCAGCCTGGATGTATACGCGGGGGGATAAAATTTACGGCGTTTATCCC‑ATTGATCCTTTATTGGTCC‑CCT‑ACCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGAAAAATTAGTTCGCTAA‑TCACCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGTTTGAT‑TTGTTGTC‑TT‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcat                                                                                                                                                                                   <  2:436470/239‑1 (MQ=255)                     aaaaaaaaGATATTTCAGCTGGGATTTATGCGGGGGGGAAAAAAATTACCGGGGTTTATACC‑ATTACTCTTTTTTGGTTCC‑CCT‑AACAAAAAGAAAAAGCCGAG‑GGGA‑AAAGAAAAAAATTTGTTCGCT‑A‑TCCGCCTT‑TTATA‑TTTTTTGTTAAAACGGGGGTTTAAT‑TTGTTGTT‑TC‑AAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCg                                                                                                                                                                             <  1:1077744/237‑1 (MQ=255)                                    cagatggtttGTATCCGCGGGGGGATAAAATTTGCGGCGTTTATACC‑ATTTATCTTTTATTGTTCC‑CCT‑ATCAAAAAGAAGAAGCCGAC‑GAGA‑TAAGAGAAAATTTAGTCCGTAAC‑TCAGCC‑T‑TTATT‑TTTTTTTTTTAACCGGGCGTTTGAT‑TTTTTTTC‑TA‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAccg                                                                                                                                                              <  1:477713/230‑1 (MQ=255)                                                   cgcgggGGAAAAAAATTAAGGGGGTTTATACA‑ATTGACCTTTTATGGTCCC‑CCT‑ACCAAAAAAAAGAAGCCGAG‑GAGC‑TAAGAG‑AAAATTTATTCCCTAATTCAGCC‑T‑TTATT‑TTTTTTTTTAAAACGGGGGGTTAAT‑TTTTTGTC‑TC‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaaaa                                                                                                                                               <  2:653900/239‑1 (MQ=255)                                                                   ttACGG‑GGTTTTACCAATTGTCCTTTTATTGGCCC‑CCT‑A‑CCAAAAAAAAGAGCCCAA‑GAACATAAAACAAAATTTAGTCCGTAAA‑TCAGCT‑T‑TTTTA‑TTTTTTTTTAAACCGGGGGGTTAAT‑TTTTTTTC‑TC‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTc                                                                                                                               <  2:614736/239‑1 (MQ=255)                                                                   tgaggggggtacaaCC‑ATTGATCTTTTA‑TGTTCC‑CCTAAACAAAAAGAAG‑AGCCGGCAGAGA‑AAAGA‑AAAAATTTGGTCCATAA‑ACAGGC‑GTTTATT‑TTTTTTTTTAAACCGGGCGTTTGAT‑TTGTGGTT‑TT‑CAT‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTc                                                                                                                               <  1:361039/226‑1 (MQ=255)                                                                                      tGATCCTTTATTGGTCA‑CCT‑ATCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAA‑C‑CGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTTATCATCCGTACGCCCGGTTACCCCCGGGCTGGACAAAAACTGCTTTAACAATTCAg                                                                                                                             >  2:463056/1‑223 (MQ=255)                                                                                      tGATCCTTTATTGGTCA‑CCT‑ATCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAA‑C‑CGCCATGAGGAACTCAATCTGCCCACAGAGCCAGTGCCCACGGTTTATCATCATCCGTACGGCCGTTTCACGCCGGGCTGGACAAAAACTGCTTTAACAAATCAGGTAATAAAACGAggtt                                                                                                             >  2:793499/1‑236 (MQ=255)                                                                                      tGATCCTTTATTGGTCA‑CCT‑ATCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAA‑C‑CGCCATGAGAAACTCAATCTTCCCACAGCGCCCGTGCCAACGGTTTTTCATCCTCCGTCCGCTCGTTTACCGCCGGGCTGGACAAAAACTGCCTTAAAAATTCAGGGAAAAAAACGAgttg                                                                                                             >  2:895467/1‑239 (MQ=255)                                                                                       gATCCTTTATTGGTCA‑CCT‑ATCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTGTTTTC‑TT‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGGGCCAACGGTTTATCATCATCCGTACGCTCGCTTACCCCCGGGCTGGACAAAAAACGCTTTAACAATTCAGGGAATAAATCGAgttgt                                                                                                            >  1:1070438/1‑239 (MQ=255)                                                                                          ccTTTATTGGTCA‑CCT‑ATCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑ATGTTGTT‑TT‑CAA‑C‑CGCCAAGAGAAACTCAATCTGCCCACAGCGCCAGTGCCAACGGTTTTTCATCATCCGTACGCTCGTTTAACGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttgtt                                                                                                           >  2:1048379/1‑237 (MQ=255)                                                                                             ttATTGGTCA‑CCT‑ATCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTTTTTTC‑TT‑CCA‑C‑CGCCATTAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACCGTTTTTCCTTCTCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAAACGCTTTAACAATTTAGGGAAAAAATAGAattttttccgc                                                                                                      >  1:604140/1‑234 (MQ=255)                                                                                             ttATTGGTCA‑CCT‑ATCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGGTGTTCCTGc                                                                                                      >  2:801298/1‑239 (MQ=255)                                                                                               ttttgggtcccttaccaaaaaaaAGAGGCCGAA‑GACA‑AAAGACAAAAATTGGTTCGATAA‑CCGCCC‑T‑TTATA‑TTTATTTTTAAACCGGGGGGTTAAA‑TTTTTTTC‑TC‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGc                                                                                                      <  2:692305/222‑1 (MQ=255)                                                                                                      a‑CCT‑ATCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTGTTTTT‑TT‑CAA‑C‑CGCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCACAGGTTTTTCATCATCCGCACGCTCGTTTACCGCCGGGCTGGGCCAAAACCGGTTTTACCATTCAGGTAATTAATCGAGTTGGTCCTGCGTttctgc                                                                                              >  1:1158122/1‑238 (MQ=255)                                                                                                      a‑CCT‑ATCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAA‑C‑CGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgctgc                                                                                              >  2:492121/1‑238 (MQ=255)                                                                                                            aTCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGTA‑TTTTTTCC‑TT‑CAA‑C‑CGCCATTAGAAACTCACTCTGCCCCCAGCGCCACGTCCAACCGTTTTTTATCAACCGTACGCCCGCTTTCCGCCCGGCTGGACCAAAAACCCTTTTACAAATCAGGTAATAAATCGAGTTGTTCCTCCGTTGCTgtcgccc                                                                                         >  1:502177/1‑238 (MQ=255)                                                                                                                       gaagCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTTAA‑TTGTTGAC‑TT‑TAA‑C‑CGCCATGAGGAAATTACTCTGCCCACAGGGTCAGTGCCAACGTTTTATCATCATCCGTTCGCTCGTTTACCGCCGGGGTGGACAAAAAATGCTTTTACAATTCATGGACTAAATCGAGTTTTTCCTGCGGTGCTTGCTCCCGTCCATTGata                                                                              >  1:252931/1‑237 (MQ=255)                                                                                                                        aagCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTTTTTTC‑TT‑TAA‑A‑CGCCATGAGAAACTCACTCTGCCCACAACGCCAGTGCCAACGTTTTATTATCATCCGCAAGCCCGTTTACCCCCGGGCTGGACACAAAATGCATTATAAATTTAGGGAATAAAACGAGTTGTTTCTGCGTTGCTGCCGCCAGTCCAtttcaaa                                                                             >  1:808528/1‑234 (MQ=255)                                                                                                                         agCCGAC‑GAGT‑AAAAAAAAAATTTAGTCCTTTAA‑CCACCC‑T‑TTATT‑TTTTTTTTTTAATCGGGC‑GTTTAA‑TTTTTTTC‑TTCCAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaaa                                                                            <  2:604140/239‑1 (MQ=255)                                                                                                                                acgataaaagaaaaaatTTAGTTCCCTAA‑CCACCC‑G‑TTATA‑TTTTTTTTTAAAACGGGCGGTTAAT‑TTGTTGTT‑TA‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACggg                                                                        <  2:808528/225‑1 (MQ=255)                                                                                                                                 ggcaagaaaacacaaaTTAGCTCTCTCA‑TTCCCC‑T‑TTATA‑TTT‑TTGTTGTAACTGGCGGTTTAAGTTGTTGTC‑TA‑CCA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCACCCGTTCGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACgggg                                                                       <  1:895467/225‑1 (MQ=255)                                                                                                                                      tAAAAGAAAATTTTGTTCG‑TTA‑TCATCC‑G‑TTTTTTTTTATTGTTCAGACGGGCGTTTGAT‑TTGTTGTC‑TC‑CAC‑C‑CGCCATGAGAAACTCACGCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGTCAGTCCATTGAAAAACGGGGTaaaga                                                                 <  1:401509/239‑1 (MQ=255)                                                                                                                                               aaTTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAA‑C‑CGCCATTAGAAACTCACTCTGCCCACAGCGCCAGTGCCAAAGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCCGGCCATTTAAAAAAGGGGTAAAGAAAGGca                                                           >  2:858271/1‑236 (MQ=255)                                                                                                                                                    gTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTTTTGTT‑TT‑CCA‑C‑CCCCATGAGAAACTCACTCTTCCCCCAACCCCCGTGCCCACGTTTTTTCCTCATCCGTACGCTCGTTTTCCCCCCGGCCGGACAAAAAATTCCTTTACAATTTAGGTAATTAAACGAGGTGTTCCTTCGTTGCTTCCGCCCATCCCTTTAAAAAAAGGGGAAAGAAAAGCACATgaggg                                                   >  1:1073900/1‑238 (MQ=255)                                                                                                                                                     ttCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑TCA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCACCGTTTTATCATCATCCGTACTCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCCAGTTGTTTCTGCGTTTCTGCCGCCCGTCCCTTGGAAAAAGGGGGTAAGACAGGCACCTgagggg                                                  >  1:863741/1‑237 (MQ=255)                                                                                                                                                     tattcctctaacccccttttATT‑TTTTTTGTTAAAGGTGG‑GGTTTGA‑TTGTTTTT‑TTCCAA‑C‑CGCCATGAGAACCTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                     <  1:801298/221‑1 (MQ=255)                                                                                                                                                         cctaaataGCC‑T‑TTTTT‑TTTTTTGTTTAAACCGGGGGTTTAA‑TTTTTTTC‑TC‑CAACC‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGaa                                                <  1:343433/232‑1 (MQ=255)                                                                                                                                                              aaccagccttatta‑ttttttGTTAAACTTGGTGTTTAAT‑TTGTTTCC‑TT‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTAt                                            <  2:805712/227‑1 (MQ=255)                                                                                                                                                                 gCC‑T‑TTATA‑TTTTTTGTTAAAACTGGGGGTTGAA‑TTGTTGTC‑TA‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                                       <  2:1113816/239‑1 (MQ=255)                                                                                                                                                                   tctttttta‑ttttttGTTAAACCTGGTGTGTGAT‑TTTTTTTT‑TA‑CAC‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                                       <  2:140672/232‑1 (MQ=255)                                                                                                                                                                    cctttATA‑TTTTTTTTTTAAAGTGGCGGTTAAT‑TTTTTGTC‑TC‑AAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCa                                      <  2:944003/236‑1 (MQ=255)                                                                                                                                                                            tttttttGTTTAACGTGGGGTTTGAT‑TTGTTTTC‑TC‑AAC‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTaa                               <  2:863741/238‑1 (MQ=255)                                                                                                                                                                               tattGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAA‑C‑CGCCATGAGAAAATCAATCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCCTACGCCCGTTTACCGCCGGGCTGGGCAAAAAATGCTTTAACAATTTAGGTAATAAATCGAGGTTTTCCTGCGTTGCTGCCGCCCGTCCATTGAAAAACGGGGTAAAGAAAAGCACATGAGGGGAATTATCCAGCATATTTaaag                             >  1:514476/1‑236 (MQ=255)                                                                                                                                                                                attGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑TAA‑C‑CGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTATCATCATTCGTAAGCTCGTTTACCCCCGGGCTGGACAAAAACTGCTTTAACAATTTAGGAAATAAATAGAGTTGTTCCCGCGTTGCTGCCGCCAGTTCATTCAAAAAAGGGGTAAAGAAAGGCACATGAGGAGGACTATTCCGCATCTTTAACGc                            >  2:60681/1‑238 (MQ=255)                                                                                                                                                                                      tAAACTGGCGGTTGAA‑ATGTTGTC‑TT‑CAA‑C‑CGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTTAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAAAGCCAGTc                       >  2:442125/1‑237 (MQ=255)                                                                                                                                                                                          cggcgggTTAAA‑TTTTTCTC‑TC‑CAC‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCAGCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGa                 <  2:1158122/234‑1 (MQ=255)                                                                                                                                                                                           ggggtgTTTAA‑TTTTTTTCTTT‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGa                 <  2:302663/234‑1 (MQ=255)                                                                                                                                                                                              gtggtgaattttttgtC‑TT‑AAC‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATa               <  2:197106/226‑1 (MQ=255)                                                                                                                                                                                              cGGGTGAA‑TTTTTGTC‑TT‑CAA‑A‑CGCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCCACGTTTTATCATCATCCCTACGCTCGCTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTTCTGCGTTGCTGCCGCCCGTCCCTTGAAAAACGGGGTGAAACAAAGCACATAGGGGGAATTCTCCAGCATATTTAACGCCAATTGCTgggaacc             >  1:470340/1‑234 (MQ=255)                                                                                                                                                                                                          ttgtC‑TA‑CAC‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAtg  <  2:514476/239‑1 (MQ=255)                                                                                                                                                                                                            gtC‑TT‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATCCCATCCGTAAAt   >  1:406874/1‑236 (MQ=255)                                                                                                                                                                                                              |                                                                                                                                                                                                                                               TGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACC‑ATTGATCCTTTATTGGTCA‑CCT‑ATCCAAAAGAAGAAGCCGAC‑GAGC‑TAAGAGCAAAATTAGTTCGCTAA‑TCAGCC‑G‑TTATA‑TTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAA‑C‑CGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATG  >  NC_000913/2151496‑2151925

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