BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,745,397	0	A	C	32.8%	51.9 / 13.0	63	L198F (TTA→TTC)	yiaN	2,3‑diketo‑L‑gulonate:Na(+) symporter ‑ membrane subunit
Reads supporting (aligned to +/- strand): ref base A (16/25); new base C (0/20); total (18/45)
Fisher's exact test for biased strand distribution p-value = 5.75e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCA  >  NC_000913/3745148‑3745631

cgccgccgGTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTAc                                                                                                                                                                                                                                               <  2:183305/251‑1 (MQ=255)

   cgccgTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑AGGGCG‑GCACG‑CCGATGCTTACATCA‑TGGGGGC‑AGGACAGGCGCCTTAA‑TTTAcc‑cc                                                                                                                                                                                                                                           >  1:874328/1‑251 (MQ=255)

           ccgttctACTGGGGCGGATGATGCGCAGCGCCAACTATCCGGTCACCCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTTTTCCCCCTTCCATCCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTTTTTATGGCGGGCATT‑GCCCCGGGCATGATGA‑TGGGCG‑CAACG‑TTGATGCTTCCCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTCCC‑CCGCCAGCaa                                                                                                                                                                                                                                   <  1:93391/247‑1 (MQ=255)

           ccGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGACAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑TCCCCCGGCATGATAG‑TGGGGG‑CACCA‑ATGGTGCTTACCTGG‑ATGAGGG‑ACGCACGTCGCCCTAA‑TTTACC‑CCCCCAGCa                                                                                                                                                                                                                                    >  1:372702/1‑250 (MQ=255)

             gCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATACATCAGTAAGCTGTTTATGGCGGGCATT‑TCCCCCGGCAGTAGGA‑TGGGCG‑GCACC‑CCGATTCTTACCTGG‑GGGTGGGGAGGCCAG‑CGCCCTAA‑TTTAcaccccccacaaac                                                                                                                                                                                                                                  >  2:909420/1‑238 (MQ=255)

                  aCTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGGTGG‑TGGGCG‑CAAAG‑CTGGTTCTTTCCTGG‑GGGGGGG‑CGGCCAGTCGGCTTTA‑TTTACA‑CCGCCAACAAAAAGCa                                                                                                                                                                                                                             >  1:170630/1‑250 (MQ=255)

                     ggTGCCGATGTTGCGCAGCGCCAATTATCCGGTCCACCGGGCGGCGGGGCTGATTTCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATG‑GCCCCTGGCATGAGGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAT‑TTTCCC‑CCGCCAGCAAAAAGCAACGa                                                                                                                                                                                                                         <  1:641727/251‑1 (MQ=255)

                         ccGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCCTGATTG‑TGGGCGTCAAC‑‑ATGATTCTTTCCAGT‑TGGTGGG‑CGGCCCGTGGGCTTAA‑TTTATC‑CCCCCAGCaaaaaaaaaccatgc                                                                                                                                                                                                                      >  2:859729/1‑240 (MQ=255)

                         ccGATGAGGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCAAATTATTCCCCCTTCCTTTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCATTAAGCTGTTTATGGCGGGCTTT‑CCCCCGGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑GGGGGGC‑AGGCCATTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCa                                                                                                                                                                                                                     <  1:272095/251‑1 (MQ=255)

                                      gcgcCACCTTTCCGGCCAACCGGGGGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTTTTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGTTTATCCATCAGTAAGCTGTTTAGGGCGGGCTTT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTCCCTGG‑TGGTGGC‑AGGCCAGTCGCCTAAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCAc                                                                                                                                                                                                      <  1:677231/251‑1 (MQ=255)

                                             acaccctgccacccGGGCGGGGGGGTTGTTTGCCTGTGGCGGCTTAATCGCCCCAATTATCCCCCATTCCACCCCTTTTTTTTTCTTCGGCGTTTCCATCGGATTATCCATAAGGAAGCTGTTGTTGGCGGTCTTT‑GCCCCTGCCATGAGTA‑TGGGCG‑AACCG‑CGGATGCTAACCGGG‑TGGTGGC‑AGCCCATTCTCCTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCATTCTTTTg                                                                                                                                                                                               <  2:744420/239‑1 (MQ=255)

                                                cccGGTCACCCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGtct                                                                                                                                                                                            <  1:104655/250‑1 (MQ=255)

                                                  cGGTCTACCGGGGGGCGGTGCTGTTTCCCTCTGTCGGCTTTACCGCGCCATTTTTCCCCCCTTCCATTCCTTTTTTTATCTCCGGCTTTCCCCGGGGATTATCCACCTGTACTCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTCCCCGG‑TGGTGGC‑AGCCCATACCCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTcc                                                                                                                                                                                          <  2:392609/251‑1 (MQ=255)

                                                     ccacccGGGGGGTGGGGGTGATTCCCTCGGGCAGCATTCTCGCGCCATTTATCCCCCCCTCCTTCCCTTTTATTACCGTCGGCGTTTCCAGGGGATTTTCCATCATAAAGCTTTTTATGGCGGGCATT‑CCCCCGGGAAGTATGA‑TGGCCG‑CAACG‑TGGATGCTTACCTGG‑TGGGGGC‑AGCCCAGTCCCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCGCCGGt                                                                                                                                                                                       <  2:271444/247‑1 (MQ=255)

                                                     caaaCCGGGGGGGGGGGCTTATTTCCTCGGGCGGCTTTATCGCGCCAATTTTTCCCCCTCCCATTCCGTTTATTTCCTCCGGCGTTTCCAGCGGATTACCCATCAGTAAGCTGTTTATGGCGGGCATG‑GCCCTGGGCAGGATGA‑GGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCCCTTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGt                                                                                                                                                                                       <  2:798166/249‑1 (MQ=255)

                                                      cAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATAAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTCCCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGt                                                                                                                                                                                       <  1:762822/250‑1 (MQ=255)

                                                        aCCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCACCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATc                                                                                                                                                                                    <  1:503844/251‑1 (MQ=255)

                                                           gggCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCTTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATG‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTCCCTGG‑TGGTGGC‑AGGCCAGTCGCTTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTgg                                                                                                                                                                                 <  1:825835/251‑1 (MQ=255)

                                                                gcggGGGTGATTGCCTCTGTGGGCATTATCGCGCCATTTATTCCCCTTTCCATTCCGTTTTTAATCTTCGGCGTTCCCAGCGGATTATCAATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑AAACG‑CTGATGCTTACCTGT‑TGTGGGC‑AGGCTAGTCGCCTTAT‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCt                                                                                                                                                                            <  1:704413/251‑1 (MQ=255)

                                                                gcggGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATCCCTTTTATTATCTTCGTCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCt                                                                                                                                                                            <  1:576096/251‑1 (MQ=255)

                                                                  ggggCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATTA‑TGGGCG‑CAACG‑CTGGTGCTTACCTGG‑GGGGGGC‑AGGCCAGTCGCCTTAA‑TTTTCC‑CCGCCAGCAAAAAGGAAAGAGGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGGATATGGgcgcagt                                                                                                                                                                          >  1:604259/1‑248 (MQ=255)

                                                                   gggCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCCGGCCTGATGAGCGGGC‑‑CAACT‑CGGTTGCTTACCTGT‑GGGGGGC‑AGGCCAACCGCCTTTA‑ATTACA‑CCGCCCGAAAACAACAAACATACAAGAAAACT‑G‑GCCCTCTTTTCTTTCCGTTATATGGGCGCCGtt                                                                                                                                                                         >  1:891006/1‑251 (MQ=255)

                                                                           tGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑TGCCCTGGCATGGTGG‑AGGGGG‑CCACC‑CTGATGCTTTCCTTT‑TGGGGGC‑ATGCCAGTCGCCTTTA‑TTTACC‑CCGCCCGCAAAAAGCAACGATGCAGGGAATAT‑GTGC‑CTCTTTTTTCTCCTGTATATTGTGGCTGTTTCTTcc                                                                                                                                                                   >  1:502613/1‑249 (MQ=255)

                                                                                     ggcATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑TCCCCTGGCATGCTGT‑TGGGCG‑CAACG‑ATGATGCTTACCTGG‑GGGGGGG‑GGGCCAGTTGCCTTTA‑TTTTCC‑CCGCCCGCAAAAAGCAACGCTGCAGGGAATATCG‑GC‑CTCTTTTTTCTCCGGTATCTGGGCGCTGtttttttccgccctcctttt                                                                                                                                                       >  2:3936/1‑234 (MQ=255)

                                                                                     ggcATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑TACCCCGGCATGATGA‑TGGGCG‑CCACG‑CTGGTGCTTTCCTGG‑TGGGGGC‑AGGCCCGTCGCCTTAA‑TTTAAC‑CCCCCACCAAAAAGCAACGATGGACGAAATAT‑G‑GCCCTCTTTTGTCTCCCGGATCTAGGCGATGTTTCTTCCACTTATCattat                                                                                                                                                       >  2:777009/1‑251 (MQ=255)

                                                                                          tATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGGTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGGG‑CAACG‑CTGATGCTTACCTTG‑TGGTGGC‑CGGCCAGTCGACTTTA‑TTTACC‑CCGCCAGCAAAAAGCAACGAGGCAGGGAATCT‑G‑GCACTCTTTTTTCTCCGGTATCTTGGCGCTTTTTCTTCCCGTCATCCTTATTgggg                                                                                                                                                  >  1:73993/1‑249 (MQ=255)

                                                                                             cgcgCCAATTATTCCCCCTCCCATTACGTTTATTATCGTCGGCGTTTCCAGCGGATTATCCACCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGAGGA‑GGGGCG‑CAACG‑CTGATGCTTCCCTGG‑GGGTGGC‑AGGCCATTCGCCTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTGTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCggg                                                                                                                                                <  1:410005/250‑1 (MQ=255)

                                                                                             cccgcaaatttctCCCCCCTCCATTCCGTTTATTACCTCCGGCGTTTCCCGCGGAGTATCCAGCAGTAAGCTTTTTTGGGCGGCCTTT‑GCCCCTGGCATGATGA‑TGGGCG‑TACCG‑CTGATGCTTCCCGGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGt                                                                                                                                               <  2:96653/239‑1 (MQ=255)

                                                                                                 caaaTTTTTCCCCCTTCCATCCCTTTTATTATCTTCGGCGTTTCCAGCGGATTACCCATCAGTAAGCTGTTTAGGGGGGGCATT‑GCCCCTGGCATTATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGGGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCg                                                                                                                                           <  1:45199/249‑1 (MQ=255)

                                                                                                          cccctTTCCATTCCTTTTATTATCTTCGGCGTTCCCAGCGGATTTTCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑GGGGCC‑CAACG‑CTGAGTCTTCCCTGG‑TGGTGGC‑AGGCC‑GTTCCCTTAAATTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTAGCTCCGGcc                                                                                                                                   <  1:166039/250‑1 (MQ=255)

                                                                                                           ccctttctttcccGTTTATTTTTTTGGGGGTTTCCAGCGGATTATCCATCAGTAGGCTGTTTATGGCGTGCATG‑CCCCCTGGCATGATGA‑TGGGCG‑CACCG‑CTGATCCTTACTTGG‑TGTTGGC‑AGGCCAGCCGCTTTAT‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTg                                                                                                                                 <  1:660923/240‑1 (MQ=255)

                                                                                                           ccccTTCCATCCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTAGGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CACCG‑CTGAGGCTTCCCGGG‑GGGGGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTg                                                                                                                                 <  1:85075/251‑1 (MQ=255)

                                                                                                                 ccATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTGCCCC‑GGCGAGATGG‑TGGGGG‑CCAAG‑CTGATGCTTAACTGG‑TGGTGGC‑CGGCCAGTCGCCTTAA‑TTTACC‑CCCCCAGCAAAAAGCAACGCTGCAGGGAATCT‑G‑GCACTCTTTTGTCTCCGGGATCTGGGCGCTGTTTCTTCCAGGCATCATTATTGGCGGGGTTccgtccccgccttttttc                                                                                                                           >  1:889042/1‑234 (MQ=255)

                                                                                                                  cATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTAccc                                                                                                                          <  1:798988/251‑1 (MQ=255)

                                                                                                                   aTTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCCGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTAACTGG‑TGGGGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAACAAAAAGCAACGATGCAGGGAATCT‑G‑GCACTCTTTTGTCTCCGGTTTCTGGGCGCTGTTTCTTCCCGTCATCATTATTGGCggggttttttcccggccgttttacac                                                                                                                         >  1:746641/1‑228 (MQ=255)

                                                                                                                         tttATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCg                                                                                                                    <  1:412431/250‑1 (MQ=255)

                                                                                                                           tattatCTTCGGCGTTCCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATTA‑GGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGGGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGaa                                                                                                                  <  1:612792/250‑1 (MQ=255)

                                                                                                                              tatCTTCGGCGTTGCCGGCGAATTATCCACCAGGAAGCTGTTGATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑AAACG‑CTGAGGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTAACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGcc                                                                                                               <  1:498820/250‑1 (MQ=255)

                                                                                                                              tatCTCCGGCGTTTCCAGCGGATTACCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGcc                                                                                                               <  1:3936/250‑1 (MQ=255)

                                                                                                                                   tCGGCGTTTCCAGCGGATTTTCCATCATGAAGCGGTTTAGGGCGGGCATG‑GCCCCTGGCATTAGGA‑TGGGCG‑CACCG‑CTGTGGTTTCCCTGG‑GGGGGGC‑GGCCCTGTCGCTTTAA‑TTGCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCGCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGc                                                                                                          <  2:724311/250‑1 (MQ=255)

                                                                                                                                    gggCGTTTCCAGGGGATTATCCATCAGAAAGCTGTTTATGGCGGGCATG‑GCCCCGGGCATGATGA‑TGGGCG‑CACCG‑CTGATGCTTACCGGG‑GGGTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCg                                                                                                        <  1:331854/250‑1 (MQ=255)

                                                                                                                                      gCGTTTCCAGCGGATTACCCACCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCAGGATGA‑TGGGCG‑CACCG‑CTGATCCTTACCTGG‑TGGGGGC‑AGGCCAGCCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATGATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGtt                                                                                                      <  1:487226/251‑1 (MQ=255)

                                                                                                                                                   aTTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCGGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATg                                                                                         <  1:259320/251‑1 (MQ=255)

                                                                                                                                                        ccttctgtaggttGTTGATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑TGGAGGCTT‑CCTGGGTGGTGCC‑AGGCCAGTCGCTTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTg                                                                                    <  2:470443/239‑1 (MQ=255)

                                                                                                                                                          aTCAGTAAGCTGTTGATGGGGGGCATT‑GCCCCGGGCTTGATGA‑TTGGCG‑CAACG‑CTGATGTTTACCTGG‑TGTTGGC‑AGGCCAGTCGCTTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCTCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGtt                                                                                  <  1:691145/251‑1 (MQ=255)

                                                                                                                                                             aggaagttGTTTTTGGCGGGCTTT‑TCCGCTGGCATTAGGG‑TGGGCG‑‑AGCGTCTGATGCTTCCCTGG‑GGGGTGC‑AGGCCTGTCGCCTTGA‑TT‑ACCTCCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTCCCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGt                                                                               <  2:254376/244‑1 (MQ=255)

                                                                                                                                                                         gtggggggctgt‑GCCCCGGGCAAGATGA‑TGGGCG‑CAGCG‑CTGAGGCTTCCCTGG‑TGGGGGC‑GGGCCGGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCAGCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTa                                                                   <  1:230951/239‑1 (MQ=255)

                                                                                                                                                                          gggCGGGCATT‑TCCTCGGGCAGGATGA‑TGGGCG‑CAACC‑CTGATGCTTACCGGG‑TGGGGGC‑AGGCCAGTCGCCTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTAtt                                                                 <  1:771523/250‑1 (MQ=255)

                                                                                                                                                                                  aTT‑GCCCCTGGCAGGAGTA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGGGGC‑AGGCCATTCGCCTTAA‑TTGCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTATTTACCGTGa                                                         <  1:259651/251‑1 (MQ=255)

                                                                                                                                                                                      gCCCCTGGCATGATGC‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCattatt                                                                                                                                                      >  2:389219/1‑156 (MQ=255)

                                                                                                                                                                                      gCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCattatt                                                                                                                                                      <  1:389219/156‑1 (MQ=255)

                                                                                                                                                                                         ccTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCTCTCCGGCCTGGTTTCCCCCACCGACGCCGGGGCCGTTGCTGCTTTTTTTTCGCTGTTTGTTGCCACCAGTT‑TTTACCGTGAAATGAc                                                   >  1:580594/1‑251 (MQ=255)

                                                                                                                                                                                             ggtgtatga‑tgGGCG‑CAACG‑CGGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTATTTACCGTGAAATGACCTTc                                               <  1:449002/246‑1 (MQ=255)

                                                                                                                                                                                                     a‑tgGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCGGCCTGTTTTCCACCACCGCCGCCGGGGGCGGTGTTGCCTTTTCTTCGCTGTTTGTCGTCCC‑ACTTCTTTTCCGTGGAATGACCtttcccccacc                                       >  2:12445/1‑242 (MQ=255)

                                                                                                                                                                                                          gcg‑cAACG‑CTGATGCTTCCCTGG‑TGGTGGC‑AGCCCAGTCGCCTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTATTTACCGTGAAATGACCTTCGCCACACTTTgg                                   <  1:272525/251‑1 (MQ=255)

                                                                                                                                                                                                             ccAACG‑CTGTTGCTTCCCGGG‑TGTTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCa                                 <  1:350237/250‑1 (MQ=255)

                                                                                                                                                                                                                       gttgcttgcctgtgggtggC‑AGGCCATTCGCCTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCa                         <  1:585399/237‑1 (MQ=255)

                                                                                                                                                                                                                       aTGCTTCCCGGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTCCC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCAt                        <  1:567501/251‑1 (MQ=255)

                                                                                                                                                                                                                         gCTTACCTGG‑TGGTGGC‑AGGCAAGTCTCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTg                      <  1:796874/251‑1 (MQ=255)

                                                                                                                                                                                                                                 gg‑tggtggC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTTCTCCCGCCTGTTTTCCCCCACCCAAGCCGGGGGCGTTGCTGCTTTTTTTGCGCTGTTTTTCGCCCC‑AGTTTTTTTCCCTGAAATGACCCTCGCCCCACTTTTGCATGTACTCATTTGCgcgggg              >  1:843381/1‑249 (MQ=255)

                                                                                                                                                                                                                                 gg‑tgggggC‑AGTCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCgcgg                <  1:325793/249‑1 (MQ=255)

                                                                                                                                                                                                                                           gcagcccatgcgctttattgtACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAacca       <  1:141273/231‑1 (MQ=255)

                                                                                                                                                                                                                                               ccAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCCGGCTGTTTTCCCCTACCCTAGCCGGGGCCGTTGCTGCTTTTTATGCGCTGTTTGTTGCCCC‑AGGTATTTACCGTGAAATGACCTTCGCCACACTTTGGGCTGTACTCATTTGCGCGGCGAAAACCCCCTcc  >  1:609948/1‑250 (MQ=255)

CGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATT‑GCCCCTGGCATGATGA‑TGGGCG‑CAACG‑CTGATGCTTACCTGG‑TGGTGGC‑AGGCCAGTCGCCTTAA‑TTTACC‑CCGCCAGCAAAAAGCAACGATGCAGGAAATCT‑G‑GCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCAC‑AGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCA  >  NC_000913/3745148‑3745631

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