BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,838,044	+G	coding (498/1458 nt)	ascF →	PTS system enzyme IIBC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,838,043	1	.	G	99.5%	768.4 / NA	212	A166G (GCG→GGG)	ascF	PTS system enzyme IIBC
Reads supporting (aligned to +/- strand): ref base . (1/0); major base G (99/111); minor base T (0/1); total (100/112)
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.99e-01

Only 100 of 236 total aligned reads displayed.

GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑C‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTG  >  NC_000913/2837951‑2838131

gaTGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg‑c                                                                                           <  1:5744840/93‑1 (MQ=255)

gaTGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg‑c                                                                                           <  1:907338/93‑1 (MQ=255)

 aTGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg‑c‑g                                                                                         >  1:3674301/1‑93 (MQ=255)

  tGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg‑c‑gg                                                                                        <  1:2069231/93‑2 (MQ=255)

    aGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGCTGGTcg‑c‑ggca                                                                                      >  1:34941/1‑90 (MQ=255)

    aGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg‑c‑ggca                                                                                      <  2:1857095/93‑4 (MQ=255)

     gcggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg‑c‑ggcat                                                                                     <  1:1029836/93‑5 (MQ=255)

     gcggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg‑c‑ggcat                                                                                     <  1:6615021/93‑5 (MQ=255)

     gcggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg‑c‑ggcat                                                                                     >  2:4225952/1‑89 (MQ=255)

       ggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCt                                                                                    <  1:2457630/93‑1 (MQ=255)

       ggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCt                                                                                    >  2:3278144/1‑93 (MQ=255)

        gcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGAGATTGGTGACGGTGCTTTCTTTTTCCTGCCGCTGATGGTCG‑CGGCATCTg                                                                                   >  2:1561843/1‑93 (MQ=255)

         cgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGc                                                                                  <  1:3230090/93‑1 (MQ=255)

         cgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGc                                                                                  <  2:250264/93‑1 (MQ=255)

         cgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGc                                                                                  <  1:2562922/93‑1 (MQ=255)

         cgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTATCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGc                                                                                  <  1:36567/93‑1 (MQ=255)

          gTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGcc                                                                                 <  2:4663620/93‑1 (MQ=255)

          gTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGcc                                                                                 >  2:5656891/1‑93 (MQ=255)

            gCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCcgc                                                                               >  1:6317697/1‑93 (MQ=255)

              tGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCcgcca                                                                             >  1:189351/1‑91 (MQ=255)

               gACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCcgccat                                                                            >  2:4065073/1‑90 (MQ=255)

               caCAAAACGATCGCCGACTTTACCTATTCTGAATGTGATTGGTGACGGTCCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCcgccat                                                                            <  2:539584/92‑4 (MQ=255)

                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCa                                                                         <  1:1242281/93‑1 (MQ=255)

                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCa                                                                         >  2:1148556/1‑93 (MQ=255)

                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCAGATGGTCG‑CGGCATCTGCCGC‑CATCa                                                                         <  2:1656900/93‑1 (MQ=255)

                   aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGAG‑GCATCTGCCGC‑CATCaaa                                                                       >  1:771452/1‑93 (MQ=255)

                    aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAAt                                                                      >  1:1445986/1‑93 (MQ=255)

                     aGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAAtt                                                                     <  2:1191355/93‑1 (MQ=255)

                      ggATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAAttt                                                                    >  2:5015562/1‑93 (MQ=255)

                       gATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTa                                                                   <  1:5708410/93‑1 (MQ=255)

                        aTCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTaa                                                                  >  2:6520158/1‑93 (MQ=255)

                         tCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTaaa                                                                 <  1:5154213/93‑1 (MQ=255)

                            ccGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAAcc                                                              >  1:5300070/1‑93 (MQ=255)

                            ccGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAAcc                                                              <  2:1992510/93‑1 (MQ=255)

                             cGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCa                                                             <  1:4830055/93‑1 (MQ=255)

                              gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCaa                                                            >  1:837204/1‑93 (MQ=255)

                               aCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAat                                                           >  1:5447644/1‑93 (MQ=255)

                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                          <  1:5907209/93‑1 (MQ=255)

                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                          >  2:5111875/1‑93 (MQ=255)

                                 cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAatat                                                         <  2:882737/93‑1 (MQ=255)

                                 cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAatat                                                         <  1:4912955/93‑1 (MQ=255)

                                  ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                        <  2:2439622/93‑1 (MQ=255)

                                  ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                        >  2:48991/1‑93 (MQ=255)

                                   tAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGt                                                       <  2:598040/93‑1 (MQ=255)

                                    aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTc                                                      <  2:571930/93‑1 (MQ=255)

                                    aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTc                                                      <  2:3085683/93‑1 (MQ=255)

                                    aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTc                                                      <  2:5270663/93‑1 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                     >  1:1121061/1‑93 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                     >  1:3195829/1‑93 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                     <  2:2224174/93‑1 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                     >  1:3929596/1‑93 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                     >  1:4375389/1‑93 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                     >  2:5125534/1‑93 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                     >  2:6212038/1‑93 (MQ=255)

                                       cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                                   <  2:5482050/93‑1 (MQ=255)

                                       cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                                   <  1:4840608/93‑1 (MQ=255)

                                       cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                                   <  1:2903588/93‑1 (MQ=255)

                                        aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                                  >  2:1086887/1‑93 (MQ=255)

                                        aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                                  >  1:164991/1‑93 (MQ=255)

                                         ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                                 <  1:2764646/93‑1 (MQ=255)

                                             gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAt                                             <  1:2507225/93‑1 (MQ=255)

                                              aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAtt                                            <  2:831102/93‑1 (MQ=255)

                                              aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAtt                                            <  2:2793341/93‑1 (MQ=255)

                                               aTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTg                                           >  1:2859109/1‑93 (MQ=255)

                                               aTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTg                                           >  1:5835076/1‑93 (MQ=255)

                                                 gtgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGcc                                         >  1:1389054/1‑93 (MQ=255)

                                                    aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                      <  2:5948955/93‑1 (MQ=255)

                                                     ttGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                                     >  2:1913889/1‑93 (MQ=255)

                                                     ttGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                                     >  2:3896327/1‑93 (MQ=255)

                                                     ttGCTGACGTTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                                     <  2:3983905/93‑1 (MQ=255)

                                                      tGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGTGATTGTCATcgc                                    >  2:149772/1‑93 (MQ=255)

                                                         tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCggg                                 <  2:2970780/93‑1 (MQ=255)

                                                         tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCggg                                 <  2:4419336/93‑1 (MQ=255)

                                                            cGGTTCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgt                              <  1:1813545/93‑1 (MQ=255)

                                                             ggTGATTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGtgtg                             <  1:553101/93‑1 (MQ=255)

                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                            >  1:6377968/1‑93 (MQ=255)

                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                            <  2:4302719/93‑1 (MQ=255)

                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                            <  2:3752603/93‑1 (MQ=255)

                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                            >  1:4168491/1‑93 (MQ=255)

                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑AATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                            <  2:4578980/93‑1 (MQ=255)

                                                                gCTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTg                          >  2:1289144/1‑93 (MQ=255)

                                                                 cTTTCTTCTTCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTgg                         >  1:5665842/1‑93 (MQ=255)

                                                                  tttCTTCTGCCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt                        <  2:588161/93‑1 (MQ=255)

                                                                      ttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                    >  2:4334217/1‑93 (MQ=255)

                                                                      ttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                    >  2:6359462/1‑93 (MQ=255)

                                                                      ttcttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                    <  2:2079163/93‑1 (MQ=255)

                                                                        ctgccTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATcc                  <  2:3503920/90‑1 (MQ=255)

                                                                         ttcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCg                 <  1:3057592/93‑1 (MQ=255)

                                                                          tcCTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATTAGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa                <  2:3548047/93‑1 (MQ=255)

                                                                            cTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGt              <  1:2662510/93‑1 (MQ=255)

                                                                            cTGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGt              >  1:6234818/1‑93 (MQ=255)

                                                                             tGCCGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtt             <  1:4566286/93‑1 (MQ=255)

                                                                                cGCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGttttt          >  1:4614813/1‑92 (MQ=255)

                                                                                 gCTGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAt         >  1:5022313/1‑93 (MQ=255)

                                                                                   tGATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTg       >  2:5424567/1‑93 (MQ=255)

                                                                                    gATGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGa      <  1:5875507/93‑1 (MQ=255)

                                                                                     aTGGTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGaa     >  2:4605312/1‑93 (MQ=255)

                                                                                        gTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGCGTTTTATTGAACTg  >  1:3326259/1‑93 (MQ=255)

                                                                                        gTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTg  <  2:5769085/93‑1 (MQ=255)

                                                                                        gTCG‑CGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTg  >  2:2194514/1‑93 (MQ=255)

GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCG‑C‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTG  >  NC_000913/2837951‑2838131

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