BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,804	T→G	L192L (CTT→CTG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,804	0	T	G	98.4%	188.5 / NA	62	L192L (CTT→CTG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand): ref base T (1/0); major base G (39/21); minor base A (1/0); total (41/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGAT  >  NC_000913/4618555‑4619050

tCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCC‑ACTGATCACCGCCTCTATTCTGGCGAAGAAACTTGc                                                                                                                                                                                                                                                          >  1:334433/1‑251 (MQ=255)

    atgctctcTGGTCGCGGCCCCGGTCATACTGGCGGTACGCGCGACAAACTGGAATCCACCCCGGGTTTCGAAATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATAAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACGGGCTCCGGGGGATAAACGTTTCTCCG‑CGCCCCGTGATATTCCCGAAACCGTGGAC‑TCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGa                                                                                                                                                                                                                                                       <  1:265777/247‑1 (MQ=255)

       atCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGCTGCCAACCGCTTCCGA‑GAAATTATTAAAGACGTCGGCGTGGCGCTTATCGGTCAGACCAGCTCCCTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑ATCATCCCGCTGATCAACGCCTCTTTTCTGGCGAAGAAAATGGCGGAAgg                                                                                                                                                                                                                                                    >  1:287387/1‑251 (MQ=255)

            gggTCGCGCCCTCGTTCTTCCTGGCTGTACGCCCGACAAACGGGAACCCATCCCGGGCTTCAACATTTCCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGGGGCGATTATCGGTCAGACCAGCTCACTGCCTCCGGCGGATAAACGTTTCTCCG‑CGCCCCGTGATATTCCCGCAACCGTGGAC‑CCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTgg                                                                                                                                                                                                                                               <  1:129507/250‑1 (MQ=255)

                  cgGCCTCGGTCATACTGGCGGTACGCTCGACAAAATGGAATCCATCCCAGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGGTAAACGTTTTTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGGAGGTCTGGgcgcgc                                                                                                                                                                                                                                         >  1:336807/1‑251 (MQ=255)

                     ccTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCCGGCTTCGACATTTTCCC‑GGATGACAACCGGTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGGTTATCGGTTCGACCAGCTCCCTGGCTCCGGCGGGGAAAAGTTTCTAAG‑CGACCCGGGATTTTAACCCAACCGTGGAC‑TCCCTCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGGCGGTCTGGgcgcgcggg                                                                                                                                                                                                                                      >  2:289626/1‑248 (MQ=255)

                      cTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGAAATTACCGCAACCATGGAC‑TCCATCCCGCTGATCGCCGGCTCAATTCTGGCGAAGAAACTGGCGGAAGGGCTTGACGCGC‑Tggg                                                                                                                                                                                                                                    >  1:378869/1‑250 (MQ=255)

                               cctggcggaccgtccaaCAATCTGAAATCCCTCCCTGTCTACGAC‑TTTTCCCTGTATGACAACCGTTCCCGC‑AAAATTATTAAAGACGCTGGCGGGGCTATTATCGGTCAAACCGGCTCATGGGCTCCGGCGGAAAACCTTTTCTACG‑CGCCCCGTTATATTACCCCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGt                                                                                                                                                                                                                           <  2:390100/235‑1 (MQ=255)

                                          gCTCGACAAACTGGAATCCATCCCTGGCTTAGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGAAATTACCGCAACCGGGGAC‑TCCATCCCGCTGATCACCGCATCTATTCTGGCGAAAAAACTGGCGGAAGGTCTGGACGCGC‑TGGGGATGGGCGTGGAAGAGGGt                                                                                                                                                                                                                 >  1:370487/1‑250 (MQ=255)

                                                  aaCTGAAACCCATCCCTGGCTTCGACATTTCCCG‑GGATGACACCCGTTCCCGC‑AAAATTATAAAAGCCGGCGGCGTGGCTATTATCGGCCAGCCCAGCTCACTGGCTCCGGCGGATAAACTTTTCTCCG‑CGCCCCGTGATATTCCCGCACCCGTGGCC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGcgcg                                                                                                                                                                                                        <  1:289626/251‑1 (MQ=255)

                                                    cTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCCGCTCACTGGCTCCGGCGGGTAAACGTTTCTACG‑CGACCCGTGGTATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGACGAAAATGGAGGAAGGACTGGACGGGCGTGG‑GGTGGACGTGGAAGAGGGGAGCGGGGCGtt                                                                                                                                                                                                      >  2:398098/1‑251 (MQ=255)

                                                         atccatccCTGACTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAACTTATTAAAGACGTCGGCGTGGCGATTATCGGTAAGACCCGCTCACTGGCTCCGGCGGATAAACGGTTCTACG‑CGAACCGTGATATTACCGCAACCGGGGAA‑GTCATCCCGGTGATACCCGCCTCTATTCTGGCGAAGAACCTGGCGGAGGGTCTGGCCGCGC‑TGGTGATGGGCGTGAAAGgggggcggggagcgtgtaggg                                                                                                                                                                                                 >  2:262608/1‑234 (MQ=255)

                                                              tccCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCTGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGGGGAC‑TCCATCCCGCTGATAACCGCCTCTTTTCTTGCGAAGAAACTGGCGGAAGGGCTGGACGCGC‑GGGCGGTGGAAGGGGAAGTGGGTAacggggggttttttcaggcc                                                                                                                                                                                            >  2:660192/1‑235 (MQ=255)

                                                                  gggttccGACATTTCCCC‑GAATGACAACCGTTTCCGC‑AAAATTATTAAAGCCGCCGGCGTGGCGATTATCGGCCAGACCAGCCCACTGGCCCCGGCGGATAACCGTTTCTACG‑CGCCCCGTGATATTCCCGCACCCTGGGCC‑CCCACCCCGCTGACCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACg                                                                                                                                                                                       <  1:124021/245‑1 (MQ=255)

                                                                     cTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAAc                                                                                                                                                                                    >  1:611348/1‑251 (MQ=255)

                                                                          aCATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCCCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGCAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGCGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAActatct                                                                                                                                                                               >  1:157354/1‑251 (MQ=255)

                                                                           cATTTTCCCCGGATAACAACC‑TTTCCCCAGAAATAATTAAAGCGGACGGGGTGGCTATTATGGGTCAGCCCTGCTCACTGGCGCGGGCGGA‑AAACCGTTCTTCGGCGCCCCGTGATATTCCCGCACTCG‑GGACCTCCATCCTGCTGTCCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAActctct                                                                                                                                                                               <  2:256210/251‑1 (MQ=255)

                                                                            aTTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTg                                                                                                                                                                              >  1:103083/1‑250 (MQ=255)

                                                                                 ccc‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑Tg                                                                                                                                                                                                                                      >  1:594691/1‑190 (MQ=255)

                                                                                              ccGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGGCGCGC‑TGGCGATGGACGTGAAAGTGGGTAGCGGCGCGGTG‑TTGCCGGCCTACGAGATCTCTGAAACCCCTGCCGGA‑GCgg                                                                                                                                                           >  2:209354/1‑250 (MQ=255)

                                                                                               cGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCCACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGACGGTCTGGTCGCGC‑TGGGGATGGACGTGAAAGGTGGTAGCGGCGGGATT‑AATCCGACCTACGAACTCCCTCAAGCCCTTGCCggagggat                                                                                                                                                           >  2:204290/1‑244 (MQ=255)

                                                                                                             ttattaAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGATCACTTGCTCCGGCGGATAAACGGTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCACCTCTATTCTGGCGAAGAAACTGGCGGCAGGGCTGGCCGCGC‑TGGTGATAGACGTGACAGTGGGTAGCCGCGCGGTT‑ATTCCGACCTACGAACTCCCTGTAGCCCTGGCCGGA‑ACGGTTGTTGGCGTCGCt                                                                                                                                             >  2:549798/1‑251 (MQ=255)

                                                                                                                  aaaGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGGCCTACGAACTCCCTGAAGCCCTTGCCGCA‑GCGGTTGTTGGCGTGGCCAACgg                                                                                                                                        >  2:577704/1‑251 (MQ=255)

                                                                                                                   aaGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGAGGAAGGTCTGGACGCGC‑TGGAGATGGACGTGAAAGTGGGTAGCGGCGGGTTT‑ATGCCGACCTACGGACTCTCTGAATCCCTTTCCGGA‑ACGGTTGTTGGCGTGGCTATCGgc                                                                                                                                       >  2:150452/1‑251 (MQ=255)

                                                                                                                     gACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGGGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTTAT‑CCGGCCCACGAACTCTCTGAAGCCCTTTCCGAA‑GCGATTGTTGGCGTGGCTAACggtcc                                                                                                                                     >  2:225109/1‑248 (MQ=255)

                                                                                                                      aCGTCGGCGGGGCGGTTATCGGTCAGACCAGCTCACTGGCTCAGGCGGATAAACGTTTCTACG‑AGACCCGTGATATTACCGCAACCGTGGCA‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGAGTGGTCGCGC‑TGGTGATGGTCGTGAAAGTGGGGAGCGGCGCGTTT‑CTGGAGACGCACGAAATCTCTGGAGCCCTTGCCGAAAGC‑ATTGTTGTCGTggtcacatgcgcg                                                                                                                                    >  2:6286/1‑240 (MQ=255)

                                                                                                                       cggcggcggGGCGTTTATCGGTCAGCCCACCTACCTGGCTCCGGGGGAAAACCTTTTCTACG‑CGACCCGTAATATTACCGCAACCGGGGAC‑TCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTg                                                                                                                                   <  1:385968/248‑1 (MQ=255)

                                                                                                                             cGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGAGGGC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAGGAAACTAGCGGAGAGTCTGGACGCGC‑TGATGAGGGACGGGAAAGGGAGTAGCGGCGAGTTT‑AGGCCGGCCTAGGAACTCTCTGAAGCCCTTTCCGAA‑GCGATTGTTGGCGTGGCTAACGGAGCTGGCGTgg                                                                                                                             >  1:1241/1‑250 (MQ=255)

                                                                                                                                    aTTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGctccac                                                                                                                     >  1:676389/1‑251 (MQ=255)

                                                                                                                                      tATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGGTTGTTGGCGTGGCTAAAGGCGCTGGGGTG‑CGCACCCCcg                                                                                                                   >  2:214461/1‑251 (MQ=255)

                                                                                                                                      tATCGGACAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGAACCGTGATATTACCGCAACCGGGGCC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGAGGAAGGTCTGGAAGCGC‑TGGTGATGGACGTGAAAGTGGGACGCGGCGAGTTT‑CTGCCGACCTCCGAACGATGTGACGCCCTGACTGAG‑GCGATTGTTAGCGTGGCTaatggagtgggagggcgaccacccg                                                                                                                    >  2:26380/1‑228 (MQ=255)

                                                                                                                                       aTCGGTCAGACCAGCTCCCTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCAGACCTACGAACTCTATGAAGCCCTTGCCGGA‑GCGGTTGTTGGCGTGGCTAACGGCGCTGGCggtcggcccaccgc                                                                                                                   >  1:99586/1‑238 (MQ=255)

                                                                                                                                          ggTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCGTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGGG‑GGCACCCCCgcgcc                                                                                                               >  2:37382/1‑250 (MQ=255)

                                                                                                                                             cAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCgcggct                                                                                                            >  1:189242/1‑251 (MQ=255)

                                                                                                                                              aGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTc                                                                                                           >  1:346346/1‑251 (MQ=255)

                                                                                                                                                aCCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCAc                                                                                                         >  1:204937/1‑251 (MQ=255)

                                                                                                                                                                 ggcggaaaaaCGTTTCTACG‑CGACCGGTGATATTACCGCAACGGTGGAC‑TCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTAc                                                                                        <  1:225109/244‑1 (MQ=255)

                                                                                                                                                                  gcgaaTAACCGTTTCTCCG‑CGCCCCGTGATATTCCCGCAACCGTGGCC‑TCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACt                                                                                       <  1:214461/247‑1 (MQ=255)

                                                                                                                                                                                 tACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACa                                                                                                    >  1:163387/1‑223 (MQ=255)

                                                                                                                                                                                             tGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCgggtg                                                             >  1:48980/1‑248 (MQ=255)

                                                                                                                                                                                                tatTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGAAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCAAGTGCAGGTAACGCgggtta                                                            >  1:27347/1‑245 (MQ=255)

                                                                                                                                                                                                 atTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGGTGAAGt                                                         >  1:675936/1‑251 (MQ=255)

                                                                                                                                                                                                           ccGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGAGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGGAGTTCGTGAAGCg                                               >  1:181289/1‑251 (MQ=255)

                                                                                                                                                                                                              gtggcccccttCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGt                                             <  1:209354/241‑1 (MQ=255)

                                                                                                                                                                                                                gcc‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCa                                          <  1:548493/249‑1 (MQ=255)

                                                                                                                                                                                                                        ccccGCTGACCTCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTg                                   <  1:493675/250‑1 (MQ=255)

                                                                                                                                                                                                                         cccgttggtccccGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGa                                  <  1:398098/240‑1 (MQ=255)

                                                                                                                                                                                                                          ccGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCAGAAGGTCTGGCCGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGcc                                                                                                                                                                                               >  2:275763/1‑95 (MQ=255)

                                                                                                                                                                                                                           ccgtGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGGGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACGCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTGCCTGACg                                <  2:611348/248‑1 (MQ=255)

                                                                                                                                                                                                                             ttgtcctccgtctctATTCTGGCGTAGAAACTGGCGGAAGGGCTGGACGCGC‑TGGTGATGGACGTGAAAGTTGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATAAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACgg                               <  2:103083/239‑1 (MQ=255)

                                                                                                                                                                                                                               tatcgccgccgCTATTCTGGCGAAGAAACTGGCTGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACAGGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTGCGTGAAGCGGTGCAGTTCCTGATGGGt                             <  2:294152/245‑1 (MQ=255)

                                                                                                                                                                                                                               aaccACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTTAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTg                            <  1:180315/248‑1 (MQ=255)

                                                                                                                                                                                                                                  cACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGGGCAGTTCCTGACgggggagt                         >  1:117209/1‑246 (MQ=255)

                                                                                                                                                                                                                                  acccGCCTCTTTTCTGGCGAAGAACCTGGCGGATGGACTTGACTCGT‑TAGTGATTGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCGACGTACTCTCGGAAGCCTTTTCCGAT‑TCGTTTGTTGGCGTGTCTAACGGCGCTGGC‑TGTCGCACCTCCGCTCTGATCACCGACATGAATCATGTACCTGCCTCCAGTGCAGGGAACGCGGTTGAAGTTCGTTACGCGGTGCAGTTCCTGTCGGTTGagt                         <  2:370487/249‑3 (MQ=255)

                                                                                                                                                                                                                                      gCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAATCGGTGCAGTTCCTGACGTGTGAGTATCg                     <  1:141819/251‑1 (MQ=255)

                                                                                                                                                                                                                                       ccTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGt                    >  1:373146/1‑251 (MQ=255)

                                                                                                                                                                                                                                                tGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGGATCGTAACCCGCGt           >  1:133928/1‑251 (MQ=255)

                                                                                                                                                                                                                                                  gCGAAGAAACTGGTGGAAGGTCGGGACGCGC‑TGGTGATGGACGGTAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCGCTGAAGCCCGTGCCGAA‑GCGATTTTTGGCGAGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTGCGGTAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCt         <  2:346346/251‑1 (MQ=255)

                                                                                                                                                                                                                                                   cGAAGAAACTGGCGGAAGGACTGGCCGAGC‑TGGTTATGGAAGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGCCCTACGCACTCTCTGCAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCCCCACCGCGCTGCTCACCGACAGGAATCCGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGCAGTTCGTGCAGCGGTGCAGTTCCTGACGGGTGGTTATCGTAACCCGCGtatg        >  2:415265/1‑248 (MQ=255)

                                                                                                                                                                                                                                                     aagaaACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGACAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCCCCGACATGCATCAAGTACTGGCCTCCCGTGCCGGTAACGCGGTTTAAGTTTGTGAAGCGGGGCAGTTCCTGACgggggggttttgtaacccgcggctggt      >  2:623405/1‑227 (MQ=255)

                                                                                                                                                                                                                                                         aaCTGGCGTAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGGCCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGAt  <  1:426169/251‑1 (MQ=255)

                                                                                                                                                                                                                                                         aaCTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGAt  >  1:416079/1‑251 (MQ=255)

TCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTG‑CGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGAT  >  NC_000913/4618555‑4619050

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