BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	671,943	T→C	A688A (GCA→GCG)	leuS ←	leucine tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	671,943	0	T	C	97.3%	467.7 / ‑0.6	150	A688A (GCA→GCG)	leuS	leucine tRNA synthetase
Reads supporting (aligned to +/- strand): ref base T (4/0); new base C (81/65); total (85/65)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 161 total aligned reads displayed.

ACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCC  >  NC_000913/671853‑672030

acgGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAGCGTTCAgcgc                                                                                       <  1:1200695/93‑1 (MQ=255)

acgGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAgcgc                                                                                       >  2:912784/1‑93 (MQ=255)

acgGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAgcgc                                                                                       >  1:2345192/1‑93 (MQ=255)

acgGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAgcgc                                                                                       >  2:4173683/1‑93 (MQ=255)

 cgGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGcc                                                                                      <  1:3142755/93‑1 (MQ=255)

 cgGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGcc                                                                                      <  1:2478717/93‑1 (MQ=255)

  ggCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCg                                                                                     >  1:3319556/1‑93 (MQ=255)

     cGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCaa                                                                                  >  2:4263332/1‑93 (MQ=255)

      gATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAAc                                                                                 >  2:2394005/1‑93 (MQ=255)

      gATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAAc                                                                                 >  1:218319/1‑93 (MQ=255)

      gATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAAc                                                                                 >  2:2298390/1‑93 (MQ=255)

        tatCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAAGGTTCAGCGCCGCAACAt                                                                               >  2:535381/1‑93 (MQ=255)

        tatCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACAt                                                                               >  1:2099585/1‑93 (MQ=255)

         atcatcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATc                                                                              >  1:1763394/1‑93 (MQ=255)

          tcatcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCa                                                                             >  1:1618043/1‑93 (MQ=255)

          tcatcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCa                                                                             <  1:4093317/93‑1 (MQ=255)

            atcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCAcc                                                                           <  1:1523006/93‑1 (MQ=255)

             tcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCt                                                                          >  1:4560638/1‑93 (MQ=255)

               ggTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCttt                                                                        <  2:1643089/93‑1 (MQ=255)

               ggTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCttt                                                                        <  1:1292980/93‑1 (MQ=255)

               ggTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCttt                                                                        <  2:1852559/93‑1 (MQ=255)

                gTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCtttt                                                                       <  1:4063637/93‑1 (MQ=255)

                gTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCtttt                                                                       >  2:1264977/1‑93 (MQ=255)

                gTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCtttt                                                                       <  1:1958744/93‑1 (MQ=255)

                  cACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGc                                                                     >  2:1992281/1‑93 (MQ=255)

                   aCTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCt                                                                    >  2:2642284/1‑93 (MQ=255)

                   aCTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCt                                                                    >  1:3942518/1‑93 (MQ=255)

                    cTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtg                                                                   <  2:2621800/93‑1 (MQ=255)

                      ttAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtgtg                                                                 >  2:446361/1‑93 (MQ=255)

                      ttAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtgtg                                                                 <  2:4002406/93‑1 (MQ=255)

                         gCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGc                                                              >  1:104214/1‑93 (MQ=255)

                         gCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGc                                                              >  2:3320618/1‑93 (MQ=255)

                          cGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCt                                                             >  1:959118/1‑93 (MQ=255)

                          cGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCt                                                             >  1:2966787/1‑93 (MQ=255)

                          cGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCt                                                             <  2:3563208/93‑1 (MQ=255)

                           gATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTc                                                            >  2:1683150/1‑93 (MQ=255)

                               gTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTaa                                                        >  1:2778617/1‑93 (MQ=255)

                                ttttATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTaaa                                                       <  1:3500133/93‑1 (MQ=255)

                                  ttATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAAcc                                                     <  2:3795996/93‑1 (MQ=255)

                                     tGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGt                                                  <  2:2927919/93‑1 (MQ=255)

                                     tGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGt                                                  >  2:2257093/1‑93 (MQ=255)

                                       cacaTCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGttt                                                >  2:3915424/1‑93 (MQ=255)

                                       cacaTCGCGACGCAGCGCTTTCTGATTGTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGttt                                                >  1:2089871/1‑93 (MQ=255)

                                        acaTCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTc                                               <  2:3138819/93‑1 (MQ=255)

                                        acaTCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTc                                               >  1:4684257/1‑93 (MQ=255)

                                        acaTCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGCGCTCGTAAACCAGTTTc                                               >  2:2118803/1‑93 (MQ=255)

                                         caTCGCGACGCAGCGCTTTCTGATTTTCAGTCATCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTcc                                              <  2:2198054/93‑1 (MQ=255)

                                         caTCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTcc                                              <  1:3693841/93‑1 (MQ=255)

                                          aTCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCa                                             >  1:3726515/1‑93 (MQ=255)

                                          aTCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCa                                             <  1:3040623/93‑1 (MQ=255)

                                          aTCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCa                                             >  1:3855524/1‑93 (MQ=255)

                                          aTCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCa                                             >  2:2221830/1‑93 (MQ=255)

                                           tCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAg                                            >  2:2982260/1‑93 (MQ=255)

                                           tCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAg                                            <  2:3610510/93‑1 (MQ=255)

                                             ccgACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGac                                          <  2:1667265/92‑1 (MQ=255)

                                              cgACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGaca                                         <  1:4267782/93‑1 (MQ=255)

                                              cgACGCAGCACTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGaca                                         <  1:449950/93‑1 (MQ=255)

                                               gACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGacac                                        <  1:3780207/93‑1 (MQ=255)

                                               gACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGacac                                        >  2:4714837/1‑93 (MQ=255)

                                                aCGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACg                                       >  2:2365878/1‑93 (MQ=255)

                                                aCGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACg                                       >  2:2695492/1‑93 (MQ=255)

                                                aCGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACg                                       >  2:1703051/1‑93 (MQ=255)

                                                 cGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGt                                      <  2:4435263/93‑1 (MQ=255)

                                                 cGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGt                                      >  2:4447303/1‑93 (MQ=255)

                                                   cAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGttt                                    >  1:1889390/1‑93 (MQ=255)

                                                   cAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGttt                                    <  2:2161793/93‑1 (MQ=255)

                                                   cAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGttt                                    <  1:2563359/93‑1 (MQ=255)

                                                   cAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGttt                                    <  2:2486363/93‑1 (MQ=255)

                                                    aGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTc                                   >  1:1026717/1‑93 (MQ=255)

                                                     gcgcTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTGTGCTGTGTGCTCGTAAACCAGATTCCAGACACGCTTCa                                  >  1:4408399/1‑93 (MQ=255)

                                                     gcgcTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCa                                  >  2:1324643/1‑93 (MQ=255)

                                                     gcgcTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCa                                  >  2:3581774/1‑93 (MQ=255)

                                                     gcgcTTTCTGATTTTCAGTCAGAGCATCACCGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCa                                  <  2:2537286/93‑1 (MQ=255)

                                                       gcTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAgg                                >  2:3951968/1‑93 (MQ=255)

                                                       gcTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAgg                                >  2:3604227/1‑93 (MQ=255)

                                                          ttCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAg                             >  1:3175529/1‑93 (MQ=255)

                                                          ttCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAg                             >  1:3175528/1‑93 (MQ=255)

                                                            cTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCg                           <  1:3017960/93‑1 (MQ=255)

                                                             tGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCgg                          >  1:4337317/1‑93 (MQ=255)

                                                               aTTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGtt                        >  1:3341212/1‑93 (MQ=255)

                                                               aTTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGtt                        >  1:612722/1‑93 (MQ=255)

                                                               aTTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGtt                        <  1:4395584/93‑1 (MQ=255)

                                                                ttttCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTa                       >  2:2937391/1‑93 (MQ=255)

                                                                  ttCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGc                     <  2:724298/93‑1 (MQ=255)

                                                                       tcagCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCtt                >  1:2455525/1‑93 (MQ=255)

                                                                        cagCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTc               >  2:4581740/1‑93 (MQ=255)

                                                                        cagCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTc               <  2:1507178/93‑1 (MQ=255)

                                                                         agCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTcc              >  1:4122967/1‑93 (MQ=255)

                                                                            gcATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCcaca           >  2:3280041/1‑93 (MQ=255)

                                                                             cATCAACGTTCAGCGCCGCAACATCGCCGTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCacac          <  1:2278101/93‑1 (MQ=255)

                                                                              aTCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACAcc         <  1:813800/93‑1 (MQ=255)

                                                                              aTCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACAcc         <  1:532666/93‑1 (MQ=255)

                                                                                  aCGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt     >  1:803478/1‑93 (MQ=255)

                                                                                  aCGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt     <  2:1857578/93‑1 (MQ=255)

                                                                                   cGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAtt    <  1:1859055/93‑1 (MQ=255)

                                                                                    gTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTc   >  2:3786281/1‑93 (MQ=255)

                                                                                    gTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTc   >  1:4080550/1‑93 (MQ=255)

                                                                                     ttCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTcc  <  2:1732797/93‑1 (MQ=255)

                                                                                     ttCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTcc  >  1:835101/1‑93 (MQ=255)

                                                                                     ttCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTcc  >  2:3893856/1‑93 (MQ=255)

ACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCC  >  NC_000913/671853‑672030

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