BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,781,145	C→G	17.2%	intergenic (+128/‑70)	lldD → / → trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,781,145	0	C	G	17.2%	242.7 / 11.3	99	intergenic (+128/‑70)	lldD/trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand): ref base C (51/29); new base G (2/16); total (53/46)
Fisher's exact test for biased strand distribution p-value = 4.94e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

CCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGA  >  NC_000913/3781047‑3781243

cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GTGCCg                                                                                                     >  2:2323163/1‑100 (MQ=255)

 ccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GTGGAGA‑GGGCCgg                                                                                                    >  2:1706201/1‑100 (MQ=255)

  cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCggg                                                                                                   >  1:1149488/1‑100 (MQ=255)

  cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCggg                                                                                                   >  1:2857467/1‑100 (MQ=255)

      gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGa                                                                                              >  1:1397977/1‑100 (MQ=255)

       ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAg                                                                                             >  1:2263699/1‑100 (MQ=255)

       ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAg                                                                                             <  2:1397982/100‑1 (MQ=255)

        gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAg‑g                                                                                           <  2:654108/100‑1 (MQ=255)

         agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAg‑gg                                                                                          <  1:195232/100‑1 (MQ=255)

          gagGGGACCGTTCGGCGCTGTATGTACTCCCTCTCTCTGAAACGACACCGCACTCTTTTTTTCTCCCGCGCCCC‑TAC‑GGGGAGA‑GGGCCGGGG‑TGAg‑gg                                                                                          <  1:1480414/99‑1 (MQ=255)

            ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACGGCACTCTTTTTTTCTCCCGCGCCCC‑TGG‑TGGGCGA‑GGGCCGGGG‑CGAg‑gggga                                                                                       >  2:2186393/1‑98 (MQ=255)

               gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAg                                                                                    <  1:2619333/100‑1 (MQ=255)

                aCCGTGCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAgg                                                                                   <  1:2725191/100‑1 (MQ=255)

                 ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGc                                                                                  >  1:1456035/1‑100 (MQ=255)

                 ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGc                                                                                  >  1:40750/1‑100 (MQ=255)

                 ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑CCC‑GGGGAGA‑GGGCGGGGG‑GAAG‑GAGAAAAGGc                                                                                  >  2:321256/1‑100 (MQ=255)

                  cGTTCGGCGCTGTATGTCCTCCCTCCCTCTGAAACGACGCCGCACTCTTTTTTTCTCCCCCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGcc                                                                                 <  1:1582287/100‑1 (MQ=38)

                  cGTTCGGCGCTGTATGTCCTCCCGCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGcc                                                                                 <  1:1131941/100‑1 (MQ=255)

                  cGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACGGCACTCTTTTTGTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGcc                                                                                 <  1:1231060/100‑1 (MQ=255)

                   gTTCGGCGCTGTATGTCCTCCCTCCCTCTGAAACGACACGGCACTCTTTTTTTCTCCCGCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGcc                                                                                 <  1:2435896/99‑1 (MQ=38)

                       ggCGGTGTATGTCCTCCGCCACTGGGAACGGACGCGGCGCTGTTTGTTGCCCCCGCGCCCC‑GCC‑GGGGAGG‑GGGCGGGGG‑GGAG‑GGGAAAAGGCCGCAc                                                                             <  1:2542424/99‑1 (MQ=255)

                             ttgtttaCGCCCTCACCGTGAACCGACCCCGCCGTCTTTGTTTCTCCCTCCCCCC‑CCC‑GGGGAGA‑GGGCCGGGG‑GGAG‑GGGAAAAGGCCGCACTGCTCAc                                                                      <  1:2234317/95‑1 (MQ=255)

                             gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCCCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GTGAAAAGGCCGCACTGCTCAc                                                                      <  2:2099693/100‑1 (MQ=255)

                              tatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCAcc                                                                     >  2:1414302/1‑100 (MQ=255)

                              tatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCAc                                                                      >  1:1508013/1‑99 (MQ=255)

                              tatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑GGAG‑GGGAAAAGGCCGCACTGCTCAcc                                                                     >  2:755740/1‑100 (MQ=255)

                                 gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTTTTTTTTTTTTT‑TTT‑GGGGTGA‑GGGCCGGGG‑TGCG‑GGCAAAAAGACTCACTTTTCACCAAt                                                                  >  2:1124199/1‑100 (MQ=255)

                                 gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAAt                                                                  >  1:1551652/1‑100 (MQ=255)

                                 gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAAt                                                                  >  1:2778634/1‑100 (MQ=255)

                                  taCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAAtt                                                                 >  2:468099/1‑100 (MQ=255)

                                  taCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAAtt                                                                 >  2:2912692/1‑100 (MQ=255)

                                    cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAAtttt                                                               >  2:2631601/1‑100 (MQ=255)

                                    cTCCCCCCCTTTAAAACGACCCCTCCCTTTTTTTTTTCCCCCCCCCCC‑CCG‑GGGGGGA‑GGGCCGGGG‑GGAG‑GGGAAAAGGCCGCACTGCTCACCAAttt                                                                <  2:647622/99‑1 (MQ=255)

                                     tCCTTCCCTCCGCAAGGACCCCGGTCTCTTTTTTTCCCCCTCGCCCC‑TCC‑AGGAAGA‑GAGCTGGGT‑TAAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTc                                                              <  1:2712395/100‑1 (MQ=255)

                                     tCCCTCACTCTGAAACGATACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTc                                                              >  1:57557/1‑100 (MQ=255)

                                     tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATtctc                                                              >  2:494166/1‑100 (MQ=255)

                                      cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTTTTTCTCTTC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑TGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                             >  1:137907/1‑100 (MQ=255)

                                      cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                             >  2:1329630/1‑100 (MQ=255)

                                       ccTCACTCTGAAACGACACCGCACTCTTTTTTTTTTTTTTTCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGTT‑TGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                            >  1:1713387/1‑100 (MQ=255)

                                       ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                            >  1:38292/1‑100 (MQ=255)

                                            ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GTGTCGGGG‑TGCG‑GGGAAAAGGCCCCACTGCTCACCAATTTTCTTCTAt                                                        >  1:1155695/1‑99 (MQ=38)

                                            cgctGAAACGACCCCGCGTTCTTTTTTTCCCCCCCGCCCC‑CCG‑GGGGAGA‑GGGCGGGGG‑TGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                       <  1:2740623/98‑1 (MQ=255)

                                              ctGAACCCACCCCGCACTTTTTTTTTCCCCCCCGCCCC‑CCG‑GGGGGGG‑GGGCCGGGG‑TGGG‑GGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                      <  2:562943/99‑1 (MQ=255)

                                              ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                     >  2:127314/1‑100 (MQ=255)

                                              ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCA‑TCC‑GGGGAGA‑GGGCCGGGG‑TGCG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                     >  1:1559538/1‑100 (MQ=255)

                                              ataaaaCGACCCCGCCCTCTTTTTTTCTCCC‑CCCCCC‑CCC‑GGTGCAC‑AGTCCGGGG‑CGGGCGGCACAAGGCCGCACTGCTCCCCAATTTTCTGCTATcc                                                      <  1:2414339/96‑1 (MQ=255)

                                               tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                    >  2:1951849/1‑100 (MQ=255)

                                               gaaaaCCACCCCTCCCTCTTTTTTTCCCCCTGCCCCC‑CCG‑GGGGGGG‑GGGCGGGGG‑GGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                    <  2:567853/98‑1 (MQ=255)

                                                gAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAACAGGCCGCGCTGCTCACCAATTTTCTGCTATCCTGc                                                   >  2:542932/1‑100 (MQ=255)

                                                 aaccGACCCCCCCCTTTTTTTTCCCCCCGCGCCCC‑CCG‑GGGGAGG‑GGGCGGGGG‑TGGG‑GGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                      <  2:496865/93‑1 (MQ=255)

                                                 aaaCGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCGGGGG‑TGCG‑GGGCAAAGGCCGCACTGCGCACCAATTGTTTGCTATCCTGc                                                   >  1:2795470/1‑99 (MQ=255)

                                                  aaCGCCACCGCACTCTTTTTTTCTCCCGCGCCCC‑GCG‑GGGGAGG‑GGGCGGGGG‑GGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                                  <  2:150342/99‑1 (MQ=255)

                                                  aaCGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                 >  2:549721/1‑100 (MQ=255)

                                                  aaCGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGACCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                 >  2:129223/1‑100 (MQ=255)

                                                    cGACACCGCCCTTTTTTTCTCTCCCCCCCCCC‑CCC‑GGGGGGG‑GGGCCGGGG‑GGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                               <  1:2366995/100‑1 (MQ=255)

                                                        ccggcgctctTTTTTCCTCCCTCCCCCG‑CCC‑GGGAGGG‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                 <  1:446421/88‑1 (MQ=255)

                                                         gatctctctTTTTTTCCCCCTCGCCCC‑TCG‑GGGGAGA‑GGGCGGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                              <  1:2533912/91‑1 (MQ=37)

                                                         ccGCCCTCTTTTTTCCTCCCCCGCCCC‑CCC‑GGGGGGG‑GGCCGGGGG‑TAGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                               <  1:1601606/95‑1 (MQ=255)

                                                             aCTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGGAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                      >  1:793853/1‑100 (MQ=255)

                                                              cttttttttcCCCCCGCCCCCC‑CCG‑GGGGAGG‑GGGCGGGGG‑GGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                     <  2:2613656/97‑1 (MQ=255)

                                                              ctctTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCgagatcggaa                                     >  2:2533927/1‑91 (MQ=255)

                                                              ctctTTTTTCCTCCCTCG‑CCC‑TCC‑GGGGGGG‑GGGCCGGGGTTGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                     <  2:2613736/100‑1 (MQ=255)

                                                               cctTTTTTTGCCCCTCCCCCC‑CCC‑GGGGATG‑GGCCCGGGG‑GGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                     <  2:2832007/98‑1 (MQ=255)

                                                                ctTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                   >  2:1348683/1‑100 (MQ=255)

                                                                ctTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                   >  1:2687772/1‑100 (MQ=255)

                                                                 tttttttC‑CCCTCCCCCC‑CCCGGGGGAGG‑GGGCGGGGG‑GGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                  <  2:2574602/100‑1 (MQ=255)

                                                                 ttttttcCCCCCTCGCCCC‑CCC‑GGGGAGA‑GGGCCGGGG‑TGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                  <  1:2323148/100‑1 (MQ=39)

                                                                 tttttgtCCCCCTCGCCGC‑CCC‑GGGGGGA‑GCGCCGGGG‑TGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                      <  2:354536/96‑1 (MQ=255)

                                                                  tttttcCTCCCCCTCCCC‑CGC‑GGGGGGG‑GGGCGGGGT‑TGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg                                       <  1:1926261/94‑1 (MQ=255)

                                                                    ttttCTTCCTCGCCCC‑TGC‑GGTGATA‑GGGCCTGTG‑TTAT‑GTAAAAAGACCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                               <  1:692763/100‑1 (MQ=255)

                                                                    ttttCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                               >  2:119891/1‑100 (MQ=255)

                                                                    ttttCCCCCTCCCCCCTTCC‑GGGGGGG‑GGCGCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                <  1:2879766/100‑1 (MQ=255)

                                                                     tttCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTGTCTGCTATACTGCCCCCGCACTAAGGGGGCAGt                               >  1:2417306/1‑99 (MQ=255)

                                                                     tttCCCCCCCGCCCC‑GCG‑GGGGAGG‑GCGCGGGTG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa                              <  2:2687788/100‑1 (MQ=255)

                                                                        ccccccCGCCCC‑CCG‑GGGGAGA‑GGGCGGGGG‑TGGG‑GGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                           <  2:2703137/98‑1 (MQ=255)

                                                                         ccccccccccc‑TCG‑GGGGAGA‑GGGCCGGGG‑TGGG‑GGCAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                           <  2:2581359/92‑1 (MQ=255)

                                                                          ccgtCGTCCC‑TCG‑GGGGAGA‑GGCCCGGGG‑GGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                           <  1:390509/95‑1 (MQ=255)

                                                                          cccctccccc‑TCC‑CGGGGGACGGCCCGGGA‑TGGG‑GGTGAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                             <  2:1748773/91‑1 (MQ=255)

                                                                          cccccccccc‑GCG‑GGGGAGG‑GGGCGGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                          <  2:1551658/93‑1 (MQ=255)

                                                                          cccTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTGAGGGGGCAGTATGCTa                         >  1:1700339/1‑100 (MQ=255)

                                                                          cccTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                         >  1:1598642/1‑100 (MQ=255)

                                                                          cccTCCCCCC‑CCG‑GGGGGGG‑GGGCGGGGG‑GGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                           <  2:208853/98‑1 (MQ=255)

                                                                           ccTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                        >  2:1242560/1‑100 (MQ=255)

                                                                            cTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                       >  2:1850714/1‑100 (MQ=255)

                                                                            cTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                       >  1:1806323/1‑100 (MQ=255)

                                                                                cccc‑TCG‑GGGGGGG‑GGGCCGGGG‑GGAG‑GGGAACAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                           <  1:1833390/92‑1 (MQ=255)

                                                                                cccc‑TCC‑GGGGAGA‑GGGCCGGGG‑AGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                   >  2:235590/1‑100 (MQ=255)

                                                                                 ccc‑ccc‑GGGGGGA‑GGGCGGGGG‑GGGG‑GGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt                    <  1:1335379/98‑1 (MQ=255)

                                                                                   cccccgggggggg‑gggCCGGGG‑GGGG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                  <  2:1124678/94‑1 (MQ=255)

                                                                                   ccccc‑GGGGAGG‑GGCCCGGGG‑TGAG‑CGCAAAAGGCCGCACTGCTCGCCAATTTTCTTCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                  <  1:2864211/96‑1 (MQ=255)

                                                                                   c‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa                >  2:2012862/1‑100 (MQ=255)

                                                                                      cc‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt              >  1:2269162/1‑100 (MQ=255)

                                                                                       c‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACtt             >  1:2066597/1‑100 (MQ=255)

                                                                                            gtga‑gggaggGCG‑TGTG‑GGTAAAAGGCCGCACGGCTCACCAATTTTTAGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAAATCGTACttt            <  1:981832/95‑1 (MQ=255)

                                                                                                 gggCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc     >  2:1911857/1‑100 (MQ=255)

                                                                                                 ggaccGGGG‑TGGG‑GGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc          <  1:2186380/92‑1 (MQ=255)

                                                                                                  ggCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCa    >  2:1753153/1‑100 (MQ=255)

                                                                                                    ccGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAATATGCTAAACATCGTACTTTACGAACCAg   >  1:2058701/1‑99 (MQ=255)

                                                                                                    ccGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATACTGCCCCCGCACTGAGGGGGCAGTATGCTAACCATCGTACTTTACGAACCAGa  >  1:1087095/1‑100 (MQ=255)

                                                                                                    ccGGGG‑TGAG‑GGGAAAAGGCCGCACTACTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGa  >  2:1776446/1‑100 (MQ=255)

CCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCC‑GGGGAGA‑GGGCCGGGG‑TGAG‑GGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGA  >  NC_000913/3781047‑3781243

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