BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,874,579	A→G	M1T (ATG→ACG) †	yidA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,874,579	0	A	G	98.8%	527.7 / ‑2.7	164	M1T (ATG→ACG) †	yidA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (80/82); total (81/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

Only 100 of 170 total aligned reads displayed.

GCGGCAATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCA  >  NC_000913/3874489‑3874668

  ggCAATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAg                                                                                       <  1:4614263/93‑1 (MQ=255)

   gCAATCGCATTTTTAACGGCGGGGGAAATGGTGGGGTCGGGGAGCAGAAGGGGGCCCTCCATATCGATAGCAAAGAGTTTAATAAACGTGAGt                                                                                      >  2:1167967/1‑93 (MQ=255)

    cAATCGCATTTTTAACGACGGGTGAAATGGTGTGATCGGGCAGCAGATGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGtt                                                                                     <  2:4219240/93‑1 (MQ=255)

     aaTCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTc                                                                                    <  2:265013/93‑1 (MQ=255)

        cGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTcccc                                                                                 <  1:1185218/93‑1 (MQ=255)

        cGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTcccc                                                                                 <  1:4714654/93‑1 (MQ=255)

          cATTTTTAACGGGGGGTGAATTGGTGTGATCGGGCAGAAGAAGGGGCCCATCCATATCGATAGC‑ATGAGTTTAATAGCCGGGAGTTCCCCAgg                                                                              <  2:3029813/93‑1 (MQ=255)

          cATTTTTAACGGCGGGGGGAATGGGGGGCTCGCGGCGCAGAAGGGGGCCATCCCAATCGCTAACAAAGAGGTTTATAGCCGTGAGTTcccccc                                                                               >  2:412061/1‑91 (MQ=255)

           aTTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAgg                                                                              <  1:159198/93‑1 (MQ=255)

           aTTTTTAACGGCGGGGGGAATGGGGTGATCGGGGAGCAGAAGGGGGGCCTCCATATCGATAGCAATGAGGTTTATAGCCGGGAGTTTCCCCgg                                                                              >  2:62280/1‑93 (MQ=255)

            tttttAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGt                                                                             <  2:3350944/93‑1 (MQ=255)

             ttttAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTa                                                                            <  2:106792/93‑1 (MQ=255)

                  agggcggGTGAAA‑GGTGTTAATGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTgg                                                                      <  1:3766716/91‑1 (MQ=255)

                  aCGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTg                                                                       <  2:1140205/93‑1 (MQ=255)

                  aCGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTg                                                                       <  2:1140202/93‑1 (MQ=255)

                     gcggGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGtt                                                                    <  1:4401501/93‑1 (MQ=255)

                      cggGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTc                                                                   >  2:1891218/1‑93 (MQ=255)

                       gggTGAAATGGTGTGATCGGGCAGGAGAAGGGTGCCATGCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTcc                                                                  <  1:2506566/93‑1 (MQ=255)

                       gggTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTcc                                                                  >  1:4601488/1‑93 (MQ=255)

                        ggTGAAATGGTGTGATCGGGCAGAAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCt                                                                 <  2:2601130/93‑1 (MQ=255)

                         gTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTg                                                                >  2:1976121/1‑93 (MQ=255)

                         gTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTg                                                                <  1:4461441/93‑1 (MQ=255)

                           gAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGAc                                                              >  1:44166/1‑93 (MQ=255)

                           gAAATGGTGTGATCGGGCAGAAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCAGTGAGTTCCCCAGGTAGATTGGTTCCTgcc                                                              >  1:3188826/1‑91 (MQ=255)

                              aTGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCa                                                           <  1:2433829/93‑1 (MQ=255)

                              aTGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCa                                                           >  1:2045631/1‑93 (MQ=255)

                              aTGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCa                                                           <  2:664785/93‑1 (MQ=255)

                               tGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCAt                                                          >  2:4730996/1‑93 (MQ=255)

                                ggTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATg                                                         <  1:921526/93‑1 (MQ=255)

                                ggTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATg                                                         <  1:2208161/93‑1 (MQ=255)

                                  ggtgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCt                                                       <  2:31500/92‑1 (MQ=255)

                                    tgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTaa                                                     >  2:1247837/1‑93 (MQ=255)

                                    tgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTaa                                                     <  2:2753625/93‑1 (MQ=255)

                                     gATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAAc                                                    >  1:1822962/1‑93 (MQ=255)

                                          ggCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGAt                                               <  1:1423533/93‑1 (MQ=255)

                                             agcagAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGCTTcc                                            >  2:1222857/1‑93 (MQ=255)

                                              gcagAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCg                                           >  2:4427526/1‑93 (MQ=255)

                                                agAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCt                                         >  1:919493/1‑93 (MQ=255)

                                                agAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCt                                         >  1:1901943/1‑93 (MQ=255)

                                                agAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCt                                         <  2:3464523/93‑1 (MQ=255)

                                                agAAGGGGGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCt                                         >  2:2222819/1‑93 (MQ=255)

                                                 gAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGCTTCCGCTc                                        >  2:1743879/1‑93 (MQ=255)

                                                  aaGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCa                                       <  2:635310/93‑1 (MQ=255)

                                                   aGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCaa                                      >  1:2784673/1‑93 (MQ=255)

                                                    gggTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCaaa                                     >  2:2462055/1‑93 (MQ=255)

                                                        gCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATc                                 >  2:4743630/1‑93 (MQ=255)

                                                         ccatccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCa                                <  2:1532069/93‑1 (MQ=255)

                                                          catccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAg                               >  1:2966654/1‑93 (MQ=255)

                                                          catccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAg                               <  2:1817540/93‑1 (MQ=255)

                                                           atccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGt                              >  2:2737983/1‑93 (MQ=255)

                                                           atccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGt                              >  1:2012412/1‑93 (MQ=255)

                                                           atccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGt                              >  1:774685/1‑93 (MQ=255)

                                                            tccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTa                             <  1:1629310/93‑1 (MQ=255)

                                                            tccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTa                             >  1:1877912/1‑93 (MQ=255)

                                                              catATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCTTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTaca                           >  1:4331611/1‑93 (MQ=255)

                                                              catATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTaca                           >  2:2886005/1‑93 (MQ=255)

                                                              catATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGACTCCGCTCAAAAATCAGTaca                           >  1:3447866/1‑93 (MQ=255)

                                                               atatCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaa                          >  2:106113/1‑93 (MQ=255)

                                                               atatCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaa                          <  1:256121/93‑1 (MQ=255)

                                                               atatCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaa                          >  1:3957665/1‑93 (MQ=255)

                                                                tatCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaac                         >  2:661978/1‑93 (MQ=255)

                                                                tatCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaac                         >  2:4373618/1‑93 (MQ=255)

                                                                tatCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaac                         >  1:214178/1‑93 (MQ=255)

                                                                tatCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaac                         <  1:397186/93‑1 (MQ=255)

                                                                  tCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAacac                       >  1:709660/1‑93 (MQ=255)

                                                                  tCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAacac                       >  1:4109222/1‑93 (MQ=255)

                                                                   cGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACAcc                      <  1:2848752/93‑1 (MQ=255)

                                                                    gATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACAccc                     <  2:1312681/93‑1 (MQ=255)

                                                                    gATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACAccc                     >  2:1811351/1‑93 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCTGa                   >  1:4699579/1‑93 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGa                   <  2:1619016/93‑1 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCACTCAAAAATCAGTACAACACCCGa                   >  1:166441/1‑93 (MQ=255)

                                                                       aGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAg                  <  1:350651/93‑1 (MQ=255)

                                                                       aGCAATGAGTTTAAGAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAg                  <  1:574616/93‑1 (MQ=255)

                                                                        gCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAgg                 <  2:1034901/93‑1 (MQ=255)

                                                                        gCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAgg                 >  1:1948697/1‑93 (MQ=255)

                                                                          aaTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGa               >  2:3958034/1‑93 (MQ=255)

                                                                          aaTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGa               >  1:3332064/1‑93 (MQ=255)

                                                                           aTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaa              >  1:2250554/1‑93 (MQ=255)

                                                                            tGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaaa             >  1:1106394/1‑93 (MQ=255)

                                                                             gAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaaaa            <  1:2398447/93‑1 (MQ=255)

                                                                             gAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaaaa            <  1:3203860/93‑1 (MQ=255)

                                                                              aGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAg           <  2:1937104/93‑1 (MQ=255)

                                                                              aGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAg           <  2:3129427/93‑1 (MQ=255)

                                                                               gTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgg          >  2:967590/1‑93 (MQ=255)

                                                                               gTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgg          >  1:67914/1‑93 (MQ=255)

                                                                                tttAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggg         <  2:3269437/93‑1 (MQ=255)

                                                                                 ttAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgggg        >  1:4670141/1‑93 (MQ=255)

                                                                                 ttAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgggg        <  2:55856/93‑1 (MQ=255)

                                                                                 ttAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgggg        <  2:849954/93‑1 (MQ=255)

                                                                                  tAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggggg       <  2:2518124/93‑1 (MQ=255)

                                                                                  tAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggggg       <  1:969985/93‑1 (MQ=255)

                                                                                   aaTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGa      >  1:496500/1‑93 (MQ=255)

                                                                                   aaTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGa      <  2:1668879/93‑1 (MQ=255)

                                                                                   aaTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGa      <  2:1668878/93‑1 (MQ=255)

                                                                                   aaTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGa      >  2:744738/1‑93 (MQ=255)

                                                                                   aaTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGa      >  2:1412707/1‑93 (MQ=255)

                                                                                    aTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGAt     <  1:3558443/93‑1 (MQ=255)

                                                                                    aTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGAt     >  2:150205/1‑93 (MQ=255)

                                                                                       gCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCa  >  2:1126097/1‑93 (MQ=255)

GCGGCAATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCA  >  NC_000913/3874489‑3874668

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