BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,878,947	0	C	T	23.8%	116.0 / 13.2	63	E276K (GAA→AAA)	cas1	multifunctional nuclease Cas1
Reads supporting (aligned to +/- strand): ref base C (25/23); new base T (15/0); total (40/23)
Fisher's exact test for biased strand distribution p-value = 4.46e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCG  >  NC_000913/2878702‑2879195

gCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGACTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGA‑TTTTTccc                                                                                                                                                                                                                                                                   >  2:207658/1‑249 (MQ=255)

               tGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTcc                                                                                                                                                                                                                                                    <  2:447851/250‑1 (MQ=255)

                     cAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTTTA‑TTTTTCCAGCgggaacgccgcccactata                                                                                                                                                                                                                                              >  1:118086/1‑233 (MQ=255)

                     cAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCtata                                                                                                                                                                                                                                              >  1:144078/1‑250 (MQ=255)

                       gttcatttaccGAATTATTGCGGATACATCAACTACATATACCCCTGCACGTACCTCCAACAACCAGTTGGCTAATCTGCCTCGTATGCGCGGCGGACCATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACGCCGATGGCCTGCATCTCCCAGTGCAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGGGGGGCCGGCGGTTGTAAGTTCTCC‑GCGGCAAGAACggtctctgaatgc                                                                                                                                                                                                                                           >  2:42925/11‑239 (MQ=255)

                         tCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCgg                                                                                                                                                                                                                                         >  1:11891/1‑251 (MQ=255)

                         tCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTGGA‑TTTTTGCAGCGGCAAGCACGACCTCTATaaaggg                                                                                                                                                                                                                                         >  1:196188/1‑247 (MQ=255)

                                   aaTTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTG‑TTTTTCCAGCGGGAAGCCAGTCctctcttagaggaaaaaatttg                                                                                                                                                                                                                               >  2:541175/1‑233 (MQ=255)

                                    aTTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTT‑TTTTTTCAGCGGCCAACAACACCTCTCTATACGGAATCA‑ATTTgg                                                                                                                                                                                                                             >  2:338292/1‑251 (MQ=255)

                                    aTTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTgg                                                                                                                                                                                                                             <  1:162687/251‑1 (MQ=255)

                                      ttttGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTTTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAACGGAATCC‑ATTTTGCt                                                                                                                                                                                                                           >  1:119702/1‑251 (MQ=255)

                                                    caacacaatatACCCCTGCACGGACCTCCAACAAACAGATGGTTAATTTGGCTCGTAAGCGCGGAGGGTCATTTTCAGTGACCACGACCAACATACTCATTTCAGCGACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAAGGGCAACAGGCGTTGCATCTTCAGGTGTGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTacgact                                                                                                                                                                                                            <  2:564359/244‑1 (MQ=255)

                                                     aCCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGGCCTCTATAAACGGAATCA‑ATTTTGCTAA‑TGTTTTACTACTg                                                                                                                                                                                                           >  1:218512/1‑251 (MQ=255)

                                                        tACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAATCACGTCCTCTATAAACGGAAACA‑ATTTTGGTAA‑TGTTTTACTACTGCgg                                                                                                                                                                                                        >  1:507692/1‑250 (MQ=255)

                                                             ataCCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑aaaaa                                                                                                                                                                                                   <  2:373456/250‑1 (MQ=255)

                                                             ataCCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGGGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTTTCCAGCGGCAAGCACGAGCTATATAAATGGAATCA‑ATTTGGCCAA‑ACTTTTACAACTGcccaaaa                                                                                                                                                                                                    >  1:117120/1‑244 (MQ=255)

                                                                      cACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCCGCCGCCAGCACGTCCTCCTTAAACGGGATCA‑ATTTGGCTTA‑TGTTTTACTACTTCG‑AAAAataatcctccc                                                                                                                                                                                         >  2:529947/1‑243 (MQ=255)

                                                                           aCCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGcc                                                                                                                                                                                    <  1:102545/250‑1 (MQ=255)

                                                                           aCCTCAAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTTTA‑TTTTTCCCTCGGGAAACACGGCCTCTATAATCGGAAACC‑AATTGGCTCA‑AGTTTTATTACTGCG‑CAAAAAATCCCTGCA‑CGcc                                                                                                                                                                                    >  1:398269/1‑250 (MQ=255)

                                                                               ccAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTG‑TTTTTCCAACGGGCAGCCCGTCCTCCTTAAACGGGATCA‑ATTTGGGTTA‑TGTTTTTTTACTTCG‑GAAAAAATACCTGCA‑CGCCCAAc                                                                                                                                                                                >  2:590478/1‑250 (MQ=255)

                                                                               ccAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTTTCCAGCGGCAAGAACGTCCTCTTTAAACGGAATAA‑ATTTGGGTAAATTTTTTA‑TACTGCG‑CAAAATATCCCTGccacacaaac                                                                                                                                                                                 >  2:463000/1‑241 (MQ=255)

                                                                                     acCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCCGCGGCAAGCACGTCCTCTATAAACGGAATCA‑ATTTTGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑Cttt                                                                                                                                                                       >  1:304716/1‑250 (MQ=255)

                                                                                       cAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTGTA‑TTTCTCCAACGGCAAGCAAGCCCTCTATAAGCGGGATCA‑ACTTGGCTAA‑TGGTTTTCTACTTCG‑GAAAATAACCCCTCA‑CGCCCAACGG‑A‑CCTccc                                                                                                                                                                     >  1:390686/1‑251 (MQ=255)

                                                                                          aTGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTTTTCCGCGGGCAGCACGGCCCCTATAAACGGAGACA‑ATTTTGTTAA‑TGTTTTTCTACTGCGTAAAAAAATCC‑TGCA‑CGCCAAACGGGA‑CTCCCCgg                                                                                                                                                                   >  2:303049/1‑251 (MQ=255)

                                                                                                aaTCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCCCGTCCTCTATAAACGGGATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATACCTGCA‑CGCCAAAAGG‑A‑CTTCCCGGTCCGGcc                                                                                                                                                            >  1:553567/1‑250 (MQ=255)

                                                                                                       cTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTTTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATACCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑gg                                                                                                                                                  >  1:228506/1‑251 (MQ=255)

                                                                                                         cGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑gg                                                                                                                                                  <  2:119702/249‑1 (MQ=255)

                                                                                                            aaGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGA‑TTTTTCCAGCGGCAAGCACGGCCTCTATAAGAGGAATCA‑AATTGGGTAA‑AGTTTTACTAATGCG‑AAAAAAATCCCTGCA‑CGCCAAACGG‑A‑CATCCCGGCCCGGCCC‑ACC‑A‑GGGTTa                                                                                                                                              >  1:516819/1‑250 (MQ=255)

                                                                                                                 cgGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGTTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGA‑TTTTTCCACCGGCCAGCACGTCCTCCATAAGCGGGAGTA‑ATTTTGCTCA‑TTTTTTTCTTCCGCG‑AAAAATATCCCTGCAACCCCAAAGGG‑A‑‑TTCCCCGTCCgtcccccccggggtacgcggc                                                                                                                                           >  1:325537/1‑231 (MQ=255)

                                                                                                                  ggAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTTTCCCGCGGCAAACACGTCCTCTATAAACGGAGTTA‑ATTTGGCTAA‑TGTTTTACTAATGCG‑AAAAATATTCCTGCA‑CGCCAAACGG‑A‑CTTCCCCGTCCGGCTC‑ACCCA‑GGGT‑ACGAcgcg                                                                                                                                       >  1:84736/1‑251 (MQ=255)

                                                                                                                   gAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCa                                                                                                                                                                                                                                                                   <  2:353711/135‑1 (MQ=255)

                                                                                                                                      cACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGCTGTT‑ATTCTTCCGCCGGAAGCACGACCTCTATAAACGGGATTA‑ATTTTGCTAA‑TGGTTTTCTACTGAG‑AAAAAAATCCCCGCA‑CGCCCAACCG‑A‑CCTCCCGGTCCGGCCC‑ACC‑CCGGGTT‑CGACGCGCTAT‑CTCCAAAACATTC‑ggg                                                                                                                 >  1:211787/1‑250 (MQ=255)

                                                                                                                                       aCGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTa                                                                                                                <  1:510702/251‑1 (MQ=255)

                                                                                                                                         gACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTGCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTaca                                                                                                              <  2:211787/251‑1 (MQ=255)

                                                                                                                                             aaCATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGATGGA‑TTTCTCCAGCGGCAAGCACGTCCTCTTTAAGCGGAAAAACATTTGGC‑AA‑TTTTTTAATACCGCG‑GAAAATTTCCCTACC‑CGCCAAACCG‑G‑CTTCCCCGCCCGTCCC‑ACT‑C‑GGGGTACGACGCGCTAT‑CTTAAAAGCTTTC‑CGTcca‑aca                                                                                                          >  1:521901/1‑250 (MQ=255)

                                                                                                                                              aCATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCCGCGGGAAGCACGTCCTCTCTAAACGGGATCA‑AATTTGCTAA‑TGGTTTCTTACTGCG‑TAAAATATCCCTGCA‑CGCCAAACGG‑AACTTCC‑GGTCCGGCCC‑CCC‑C‑GGGGTACGACGCCCCAT‑CTCCACAGCTTTC‑CGTgca‑acc                                                                                                          >  2:461306/1‑248 (MQ=255)

                                                                                                                                                     cATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑aa                                                                                                  <  1:541175/249‑1 (MQ=255)

                                                                                                                                                                       ggCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCGCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAata                                                                              <  1:355184/251‑1 (MQ=255)

                                                                                                                                                                            gCATATCCCAGTGAAACAGGAAGCGGAATGGCAACAGTCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGGCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑GCC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGt                                                                          <  2:347436/250‑1 (MQ=255)

                                                                                                                                                                                             aGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCTGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCt                                                        <  2:437258/251‑1 (MQ=255)

                                                                                                                                                                                              ggAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAAAGG‑A‑CTTCCCGGTCCGGCCC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑AAAGTGTC‑AAAATTAATGATGTCTGCAATATCGTAAAAAAAAGAAAAAGGCtt                                                       >  1:106238/1‑251 (MQ=255)

                                                                                                                                                                                               gAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCttt                                                      <  1:554716/251‑1 (MQ=255)

                                                                                                                                                                                                 aGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CGCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGaa                                                               <  1:458885/240‑1 (MQ=255)

                                                                                                                                                                                                  gCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCAAGTCCTCTATAAGCGGAATCA‑ATTTTGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTT‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTccc                                                   >  1:418286/1‑250 (MQ=255)

                                                                                                                                                                                                      aaTGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGA‑TTTCTTCCGCGGCCAGGACGCACTCTTTAAGCGGGAACA‑ATTTGGCTAA‑TGTTTTTATACTGCG‑GAAAATATCCCTGCA‑ACCCCAACGG‑A‑CTTCCCCGTCCGGCTC‑ACC‑A‑GGGGTACGACGCGCTAT‑CTCAAAAACATTCTGGT‑CA‑ACAGGGTC‑AAATTTTATGATGTCTGCAATATCGTAAAAAAAGGAAAAAGGCTTTTCTGTAt                                               >  1:351966/1‑251 (MQ=255)

                                                                                                                                                                                                      aaTGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTGTA‑TTTTACCAGCGGCAAGAACGTCCTCCTTAAACGGGATCA‑ATTTTGTTAA‑TTTTTTTCTTCTCCG‑AAAACTATCCCTGCA‑CGCCAAACCG‑A‑CTTCCCCGTCCGGCCCCACC‑A‑GG‑ATACTACGCCCTATACTCACAAACTT‑C‑GG‑ACACACAGTGTCCAAATTTAAT‑ATTATTGCAATATCGTAAACAAAAGAAAACGACTTTCCTGtct                                               >  2:461332/1‑249 (MQ=255)

                                                                                                                                                                                                      aaTGGCAACAGGCAGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTAt                                               <  2:287619/251‑1 (MQ=255)

                                                                                                                                                                                                         ggCAACAGGCTGTGCCTCTTCATGTGGAGCCGGCGGTTATA‑TTTCTCCAGCTGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATAACCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATTATGTCTGCAATTTCGTAAACAAAGGAAGGAGGCTTTCCTTTATGc                                             <  2:477002/250‑1 (MQ=255)

                                                                                                                                                                                                         ggCAACAGGCTGTGCATCTTCAGGTGGGGCCGTCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGc                                             <  2:351966/250‑1 (MQ=255)

                                                                                                                                                                                                          gCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCCCGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGGGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAAGAAAGAGGCTTTCCTGTATGccc                                           >  1:209370/1‑249 (MQ=255)

                                                                                                                                                                                                                          cTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCCCTATAAGCGGAAACA‑ATTTGGGTAA‑TGTTTTACTAATGCG‑AAAAATATTCCTGCC‑CGCCAAACCG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCCAAAGCATTC‑GGGACA‑ACAGGGGC‑AAATTTAATGATGTCTGCAATATCGTAAAAAAAGGAAAGAGGCTTTCCTGTATGGACAAACCCAATAGCTGGt                           >  1:492856/1‑251 (MQ=255)

                                                                                                                                                                                                                                   gggCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGc                                             <  1:161536/224‑1 (MQ=255)

                                                                                                                                                                                                                                      ccGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCt               <  1:578849/251‑1 (MQ=255)

                                                                                                                                                                                                                                        gtcgtttGTA‑TTTCTCCAGCGGCAAGCACGGCCTCTATAAGCGGTATCA‑ATTTGGCTGA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCTCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGc             <  2:577933/246‑1 (MQ=255)

                                                                                                                                                                                                                                           ggTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAg          <  1:338292/251‑1 (MQ=255)

                                                                                                                                                                                                                                            gTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGt         >  2:262077/1‑251 (MQ=255)

                                                                                                                                                                                                                                             ttGTA‑TTTCTCCCGCGGCAAGCACGTACTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTAATACTGCG‑AAAAATATCCCTGAA‑CGCCACACGG‑A‑CTTCACGGTCCGGCTC‑ACC‑A‑GTGTTAAGACGCGCTAT‑CTCAACAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGAATGCACAAACCCAATCGCTGGTGCATAACAAGCTGCAAGt         >  2:474344/1‑250 (MQ=255)

                                                                                                                                                                                                                                               gTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATc      >  2:183534/1‑251 (MQ=255)

                                                                                                                                                                                                                                               gTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGAGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCC‑CGCCACACGG‑A‑CTTCCCCGTCCGGCTC‑ACC‑C‑GCGTTACGACGCTCTAT‑CTCACACGCTTTC‑GTTACA‑ACAGTGTC‑AACTTTTATGTTGTCTGCAATATCGTTATCAAACTATAGTGGCTTTCCTGTATGCACTAACCCATTATCTGGTGCATAACCAGCTGCAAGTATc      >  1:379430/1‑251 (MQ=255)

                                                                                                                                                                                                                                                 a‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgc    >  2:197622/1‑251 (MQ=255)

                                                                                                                                                                                                                                                 a‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgc    >  2:572032/1‑251 (MQ=255)

                                                                                                                                                                                                                                                 a‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGCCTCA‑TGTTTTACTACTGCG‑AACACTATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCACTATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATCGCTGGTGCATAACCAGCTGCAAGTATcgc    >  2:406994/1‑251 (MQ=255)

                                                                                                                                                                                                                                                    ttCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgcc   <  1:511958/250‑1 (MQ=255)

                                                                                                                                                                                                                                                     tctcCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgccg  <  2:106238/250‑1 (MQ=255)

GCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTA‑TTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCA‑ATTTGGCTAA‑TGTTTTACTACTGCG‑AAAAATATCCCTGCA‑CGCCAAACGG‑A‑CTTCCCGGTCCGGCTC‑ACC‑A‑GGGTTACGACGCGCTAT‑CTCAAAAGCTTTC‑GGTACA‑ACAGTGTC‑AAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCG  >  NC_000913/2878702‑2879195

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